Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,712,518 (C)6→7 intergenic (+360/‑251) 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,5121.C97.4% 128.4 / ‑3.2 38intergenic (+354/‑257)nth/dtpADNA glycosylase and apyrimidinic (AP) lyase (endonuclease III)/dipeptide and tripeptide permease A
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (22/15);  total (22/16)
Fisher's exact test for biased strand distribution p-value = 4.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

TCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACAAAGGACTT  >  NC_000913/1712377‑1712639
                                                                                                                                         |                                                                                                                               
tCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACATACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAATCCATCAccc                                                                                                                              >  2:759851/1‑139 (MQ=255)
 cAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCAcccc                                                                                                                             <  2:273590/139‑1 (MQ=255)
  aGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCAccccc                                                                                                                            <  1:78289/139‑1 (MQ=255)
                                         gACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCAcc                                                                                                                               >  1:946569/1‑97 (MQ=255)
                                         gACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCAcc                                                                                                                               <  2:946569/97‑1 (MQ=255)
                                         gACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGCACTATCTCACTAGTGAAAACGCTAAATAGCAGAACATATCGCCGATAAGCCATCACCCCCCCTCTCAGTGCAGTGCAAAAATCTGCCGTTACGtttttt                                                                                      >  2:707663/1‑139 (MQ=255)
                                           ccTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATAACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGa                                                                                    >  2:1058150/1‑139 (MQ=255)
                                               cAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGaaaaa                                                                                >  2:627853/1‑139 (MQ=255)
                                                         cacTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTg                                                                      >  2:381336/1‑139 (MQ=255)
                                                            tGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGAt                                                                   >  2:398989/1‑139 (MQ=255)
                                                              gCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGa                                                                    >  2:49053/1‑136 (MQ=255)
                                                              gCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGa                                                                    <  1:49053/136‑1 (MQ=255)
                                                              gCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATaa                                                                 >  1:1034095/1‑139 (MQ=255)
                                                                                    aTCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAAtt                                           >  1:940174/1‑139 (MQ=255)
                                                                                     tctcACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTc                                          >  2:907367/1‑139 (MQ=255)
                                                                                       tcACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTccc                                        <  2:940174/139‑1 (MQ=255)
                                                                                          cTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGc                                     <  1:627853/139‑1 (MQ=255)
                                                                                          cTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCAcc                                                                                                                               >  2:262545/1‑48 (MQ=255)
                                                                                          cTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCAcc                                                                                                                               <  1:262545/48‑1 (MQ=255)
                                                                                           tAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGcc                                    <  1:907367/139‑1 (MQ=255)
                                                                                            aGTAGAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCt                                   <  1:1026676/139‑1 (MQ=255)
                                                                                            aGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCt                                   <  1:866/139‑1 (MQ=255)
                                                                                            aGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCGCTGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCt                                   <  2:1034095/139‑1 (MQ=255)
                                                                                             gTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGc                                     <  1:278975/136‑1 (MQ=255)
                                                                                             gTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGc                                     >  2:278975/1‑136 (MQ=255)
                                                                                               aaaaaCGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCt                                   <  2:323736/136‑1 (MQ=255)
                                                                                               aaaaaCGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCt                                   >  1:204633/1‑136 (MQ=255)
                                                                                               aaaaaCGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCt                                   >  1:323736/1‑136 (MQ=255)
                                                                                               aaaaaCGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCt                                   <  2:204633/136‑1 (MQ=255)
                                                                                               aaaaaCGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAg                                >  1:990795/1‑139 (MQ=255)
                                                                                               aaaaaCGCTAAATAGCAGAACATATCGCCGTTAAACCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCGCGCCTCAg                                >  1:638921/1‑139 (MQ=255)
                                                                                                 aaaCGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGtt                              <  1:1058150/139‑1 (MQ=255)
                                                                                                   aCGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTtat                            >  1:423121/1‑139 (MQ=255)
                                                                                                    cGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTtata                           <  1:398989/139‑1 (MQ=255)
                                                                                                     gCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTcag                                                                                                                     >  1:595201/1‑48 (MQ=255)
                                                                                                     gCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTcag                                                                                                                     <  2:595201/48‑1 (MQ=39)
                                                                                                       tAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGaa                                                     <  1:800035/110‑1 (MQ=255)
                                                                                                       tAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGaa                                                     >  2:800035/1‑110 (MQ=255)
                                                                                                       tAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGt                        >  1:245744/1‑139 (MQ=255)
                                                                                                        aaaTAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTa                       >  1:159810/1‑139 (MQ=255)
                                                                                                         aaTATCAGAACCTTTCCCCGTTAATCCATCTCCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGTTAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTaa                      <  2:245744/139‑1 (MQ=255)
                                                                                                                        cGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACAAAg       >  1:537005/1‑139 (MQ=255)
                                                                                                                           cGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACAAAGGAc    <  2:802177/139‑1 (MQ=255)
                                                                                                                             ttAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACAAAGGACtt  >  2:197623/1‑139 (MQ=255)
                                                                                                                             ttAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACAAAGGACtt  >  1:769976/1‑139 (MQ=255)
                                                                                                                                         |                                                                                                                               
TCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACAAAGGACTT  >  NC_000913/1712377‑1712639

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