Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007793 725,902 G→A intergenic (+28/+560) 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,9020GA84.9% 67.6 / 3.8 33intergenic (+28/+560)SAUSA300_RS03485/SAUSA300_RS03490inorganic phosphate transporter/peptidase M23
Reads supporting (aligned to +/- strand):  ref base G (1/4);  new base A (0/28);  total (1/32)
Fisher's exact test for biased strand distribution p-value = 1.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.77e-01

GGGTGATTACATTACCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAGTATCACAAACATATGGTGGTATTGAAGTTCAGTTTTTTATGTAAGTAGATAACATATTTTATAGAATTTACGTAGAATGACCAAG  >  NC_007793/725804‑725987
                                                                                                  |                                                                                     
gggTGATTACATTACCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTCCAga                                                                                      <  1:4959164/100‑2 (MQ=255)
 ggTGATTACATTACCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAGTa                                                                                     >  1:6121327/1‑100 (MQ=255)
   tGATTACATTACCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATacc                                                                                   <  1:1311523/100‑3 (MQ=255)
    gATTACATTACCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAA‑AAAGCAAAACTGAACTTCAATacgaa                                                                                 <  1:4589493/100‑5 (MQ=255)
      ttacattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATACcagc                                                                                <  1:395854/100‑3 (MQ=255)
      ttacattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATACcacc                                                                                <  1:117429/100‑2 (MQ=255)
      ttacattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATACcacc                                                                                <  1:6005890/100‑2 (MQ=255)
      ttacattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATACcacc                                                                                <  1:4699445/100‑2 (MQ=255)
      ttacattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATACcacc                                                                                <  1:2141070/100‑2 (MQ=255)
      ttacattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATACcacc                                                                                <  1:4068494/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:3957353/100‑4 (MQ=255)
        acattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATACcaccat                                                                              <  1:2657701/100‑4 (MQ=255)
        acattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATACcaccat                                                                              <  1:2485474/100‑4 (MQ=255)
        acattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATACcaccat                                                                              <  1:683137/100‑4 (MQ=255)
        acattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATACcaccat                                                                              <  1:59093/100‑4 (MQ=255)
        acattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATACCACcaa                                                                              <  1:3280640/100‑1 (MQ=255)
        acattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATACCACcaa                                                                              <  1:672557/100‑1 (MQ=255)
        acattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATACCACcaa                                                                              <  1:1864937/100‑1 (MQ=255)
        acattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATACCACcaa                                                                              <  1:5263904/100‑1 (MQ=255)
        acattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATACCACcaa                                                                              <  1:1496830/100‑1 (MQ=255)
        acattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACAAAACTGAACTTCAATACcaccat                                                                              <  1:4862790/100‑4 (MQ=255)
         cattacCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAATacctagaca                                                                             <  1:1175324/100‑9 (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)
                                                                                                  |                                                                                     
GGGTGATTACATTACCTATTTCAGCATTGTTAGCAGGTTTACTATTCTATATACTTAACTTATTTTTCTAATTGAAAATAAAACTAAACTGAACTTCAGTATCACAAACATATGGTGGTATTGAAGTTCAGTTTTTTATGTAAGTAGATAACATATTTTATAGAATTTACGTAGAATGACCAAG  >  NC_007793/725804‑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: ‑