Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007793 725,905 T→C intergenic (+31/+557) SAUSA300_RS03485 → / ← SAUSA300_RS03490 inorganic phosphate transporter/peptidase M23

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_007793725,9050TC86.2% 50.8 / 2.9 29intergenic (+31/+557)SAUSA300_RS03485/SAUSA300_RS03490inorganic phosphate transporter/peptidase M23
Reads supporting (aligned to +/- strand):  ref base T (0/4);  new base C (0/25);  total (0/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.09e-02

TTACATTACCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAGTATCACAAACATATGGTGGTATTGAAGTTCAGTTTTTTATGTAAGTAGATAACATATTTTATAGAATTTACGTAGAATGACCAAG  >  NC_007793/725810‑725987
                                                                                               |                                                                                  
ttacattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATACcagc                                                                                <  1:395854/100‑3 (MQ=255)
ttacattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATACcacc                                                                                <  1:117429/100‑2 (MQ=255)
ttacattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATACcacc                                                                                <  1:4699445/100‑2 (MQ=255)
ttacattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATACcacc                                                                                <  1:4068494/100‑2 (MQ=255)
ttacattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATACcacc                                                                                <  1:6005890/100‑2 (MQ=255)
ttacattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATACcacc                                                                                <  1:2141070/100‑2 (MQ=255)
ttacattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATACcaca                                                                                <  1:2001830/100‑1 (MQ=255)
 tacattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATACcacca                                                                               <  1:3448189/100‑3 (MQ=255)
 tacattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATACcacca                                                                               <  1:1856260/100‑3 (MQ=255)
 tacattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATACcacca                                                                               <  1:1507294/100‑3 (MQ=255)
  acattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAGCAAAACTGAACTTCAATACcaccag                                                                              <  1:4590578/100‑4 (MQ=255)
  acattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATACcaccat                                                                              <  1:2485474/100‑4 (MQ=255)
  acattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATACcaccat                                                                              <  1:2657701/100‑4 (MQ=255)
  acattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATACcaccat                                                                              <  1:683137/100‑4 (MQ=255)
  acattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATACcaccat                                                                              <  1:59093/100‑4 (MQ=255)
  acattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATACcaccat                                                                              <  1:3957353/100‑4 (MQ=255)
  acattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATACCACcaa                                                                              <  1:672557/100‑1 (MQ=255)
  acattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATACCACcaa                                                                              <  1:5263904/100‑1 (MQ=255)
  acattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATACCACcaa                                                                              <  1:1496830/100‑1 (MQ=255)
  acattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATACCACcaa                                                                              <  1:3280640/100‑1 (MQ=255)
  acattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATACCACcaa                                                                              <  1:1864937/100‑1 (MQ=255)
  acattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACAAAACTGAACTTCAATACcaccat                                                                              <  1:4862790/100‑4 (MQ=255)
   cattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATACcaccatg                                                                             <  1:5863033/100‑5 (MQ=255)
   cattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATACcaccatg                                                                             <  1:208261/100‑5 (MQ=255)
     ttacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATACcaccatagc                                                                           <  1:4286998/100‑7 (MQ=255)
              tCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAGTATCACAAACATATGGTGGTa                                                                  <  1:859968/100‑1 (MQ=255)
                                                                  tGAAAATAAAACTAAACTGAACTTCAGTATCACAAACATATGGTGGTATTGAAGTTCAGTTTTTTATGTAAGTAGATAACATATTTTATAGAATTTACGt              <  1:366879/100‑1 (MQ=255)
                                                                            actaaactGAACTTCAGTATCACAAACATATGGTGGTATTGAAGTTCAGTTTTTTATGTAAGTAGATAACATATTTTATAGAATTTACGTAGAATGACCa    <  1:672513/100‑1 (MQ=255)
                                                                              taaactGAACTTCAGTATCACAAACATATGGTGGTATTGAAGTTCAGTTTTTTATGTAAGTAGATAACATATTTTATAGAATTTACGTAGAATGACCAAg  <  1:5700097/100‑1 (MQ=255)
                                                                                               |                                                                                  
TTACATTACCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAGTATCACAAACATATGGTGGTATTGAAGTTCAGTTTTTTATGTAAGTAGATAACATATTTTATAGAATTTACGTAGAATGACCAAG  >  NC_007793/725810‑725987

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑