Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,712,338 A→T 8.7% intergenic (+180/‑431) nth → / → dtpA DNA glycosylase and apyrimidinic (AP) lyase (endonuclease III)/dipeptide and tripeptide permease A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,712,3380AT 8.7% 157.5 / 6.4 46intergenic (+180/‑431)nth/dtpADNA glycosylase and apyrimidinic (AP) lyase (endonuclease III)/dipeptide and tripeptide permease A
Reads supporting (aligned to +/- strand):  ref base A (21/21);  new base T (2/2);  total (23/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.35e-01

CCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAA  >  NC_000913/1712193‑1712482
                                                                                                                                                 |                                                                                                                                                
ccTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGa                                                                                                                                                <  3:824885/148‑1 (MQ=255)
                     aaCATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTa                                                                                                                           >  1:119534/1‑148 (MQ=255)
                         tCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACg                                                                                                                                           >  4:82309/1‑128 (MQ=255)
                         tCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGCGATCAAGATCACg                                                                                                                                           <  3:82309/128‑1 (MQ=255)
                             tttAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACAtt                                                                                                                   >  2:489767/1‑148 (MQ=255)
                                     cAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAAt                                                                                                           <  2:119534/148‑1 (MQ=255)
                                             ccATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACAt                                                                                                                    <  4:476102/131‑1 (MQ=255)
                                             ccATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACAt                                                                                                                    >  3:476102/1‑131 (MQ=255)
                                               aTTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAGATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACAt                                                                                                                    <  3:1077644/129‑1 (MQ=255)
                                               aTTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACAt                                                                                                                    >  4:1077644/1‑129 (MQ=255)
                                               aTTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTaaaa                                                                                                 <  1:489767/148‑1 (MQ=255)
                                                ttGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACg                                                                                                                                           <  1:344902/105‑1 (MQ=255)
                                                ttGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACg                                                                                                                                           >  2:344902/1‑105 (MQ=255)
                                                                     tAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATc                                                                                                          <  4:835897/117‑1 (MQ=255)
                                                                     tAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATc                                                                                                          >  3:835897/1‑117 (MQ=255)
                                                                        tCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTaaaa                                                                                                                        <  1:716994/100‑1 (MQ=255)
                                                                        tCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTaaaa                                                                                                                        >  2:716994/1‑100 (MQ=255)
                                                                                                   gTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGt                                                                                                      <  1:1038732/91‑1 (MQ=255)
                                                                                                   gTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGt                                                                                                      >  2:1038732/1‑91 (MQ=255)
                                                                                                    ttAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCa                                                                                              >  4:736229/1‑98 (MQ=255)
                                                                                                    ttAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCa                                                                                              <  3:736229/98‑1 (MQ=255)
                                                                                                         tGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTaaaa                                                                                                                        <  2:799961/67‑1 (MQ=255)
                                                                                                         tGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTaaaa                                                                                                                        >  1:799961/1‑67 (MQ=255)
                                                                                                         tGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCa                                                                                                         >  2:71731/1‑82 (MQ=255)
                                                                                                         tGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCa                                                                                                         <  1:71731/82‑1 (MQ=255)
                                                                                                                       cGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACAtt                                                    <  2:563677/121‑1 (MQ=255)
                                                                                                                       cGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACAtt                                                    >  1:563677/1‑121 (MQ=255)
                                                                                                                                    tcattcGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTac                                                                                   >  4:276924/1‑77 (MQ=255)
                                                                                                                                    tcattcGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTac                                                                                   <  3:276924/77‑1 (MQ=255)
                                                                                                                                       ttcGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACAt                                                                                                                    >  3:449172/1‑41 (MQ=255)
                                                                                                                                       ttcGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACAt                                                                                                                    <  4:449172/41‑1 (MQ=255)
                                                                                                                                       ttcGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGcaacaa                                                                        >  3:661217/1‑85 (MQ=255)
                                                                                                                                       ttcGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTg                                                                  >  1:1024160/1‑91 (MQ=255)
                                                                                                                                       ttcGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTg                                                                  <  2:1024160/91‑1 (MQ=255)
                                                                                                                                       ttcGTGATCAAGATCACGATCCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGcagcaa                                                                        <  4:661217/85‑1 (MQ=255)
                                                                                                                                        tcGTGATCTTGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAatat                                                                    <  1:458326/88‑1 (MQ=255)
                                                                                                                                        tcGTGATCTTGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAatat                                                                    >  2:458326/1‑88 (MQ=255)
                                                                                                                                         cGTGATCTTGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAAc                                                       <  4:961728/100‑1 (MQ=255)
                                                                                                                                         cGTGATCTTGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAAc                                                       >  3:961728/1‑100 (MQ=255)
                                                                                                                                           tGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTAtctc                   >  4:161267/1‑134 (MQ=255)
                                                                                                                                           tGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTAtctc                   <  3:161267/134‑1 (MQ=255)
                                                                                                                                            gATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGa                          >  2:407876/1‑126 (MQ=255)
                                                                                                                                            gATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGa                          <  1:407876/126‑1 (MQ=255)
                                                                                                                                              tCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTaa  >  2:137581/1‑148 (MQ=255)
                                                                                                                                               cAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTAtt                                        >  4:814373/1‑109 (MQ=255)
                                                                                                                                               cAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTAtt                                        <  3:814373/109‑1 (MQ=255)
                                                                                                                                                aaGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGTTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCt                                                              <  3:908160/86‑1 (MQ=255)
                                                                                                                                                aaGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCt                                                              >  4:908160/1‑86 (MQ=255)
                                                                                                                                                 aGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTgg                           <  4:389362/120‑1 (MQ=255)
                                                                                                                                                 aGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTgg                           >  3:389362/1‑120 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
CCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAA  >  NC_000913/1712193‑1712482

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: