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.C93.2% 142.5 / 3.0 44intergenic (+354/‑257)nth/dtpADNA glycosylase and apyrimidinic (AP) lyase (endonuclease III)/dipeptide and tripeptide permease A
Reads supporting (aligned to +/- strand):  ref base . (1/2);  new base C (20/21);  total (21/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 = 9.87e-01

AATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACAAAGGACTTGTCTGA  >  NC_000913/1712386‑1712645
                                                                                                                               |                                                                                                                                     
aaTCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCAccccc                                                                                                                                  <  1:1149947/132‑1 (MQ=255)
aaTCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCAccccc                                                                                                                                  >  2:1149947/1‑132 (MQ=255)
   cATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTcagtg                                                                                                                         <  1:760807/139‑1 (MQ=255)
     tttAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCAcc                                                                                                                                     >  2:198372/1‑124 (MQ=255)
     tttAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCAcc                                                                                                                                     <  1:198372/124‑1 (MQ=255)
          aTTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGaa                                                                                                                  <  1:509292/139‑1 (MQ=255)
                                             ttACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCa                                                                                                                                        <  1:398067/82‑1 (MQ=255)
                                             ttACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCa                                                                                                                                        >  2:398067/1‑82 (MQ=255)
                                              tACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGc                                                                              >  2:294863/1‑139 (MQ=255)
                                              tACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCAccccccctct                                                                                                                             >  2:1656997/1‑87 (MQ=255)
                                              tACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCAccccccctct                                                                                                                             <  1:1656997/91‑5 (MQ=255)
                                                   tGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTcagtg                                                                                                                         <  2:1804020/91‑1 (MQ=255)
                                                   tGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTcagtg                                                                                                                         >  1:1804020/1‑91 (MQ=255)
                                                     gCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATaa                                                                       >  2:1804258/1‑139 (MQ=255)
                                                     gCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATaa                                                                       >  1:1520037/1‑139 (MQ=255)
                                                     gCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATaa                                                                       >  2:1494865/1‑139 (MQ=255)
                                                     gCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATaa                                                                       >  2:1275330/1‑139 (MQ=255)
                                                     gCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAAc                                                                      >  1:1161052/1‑139 (MQ=255)
                                                                 aaaCTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGaa                                                           <  2:1520037/139‑1 (MQ=255)
                                                                             ctcACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTcc                                               >  2:1031147/1‑139 (MQ=255)
                                                                                 cTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGc                                           >  1:1127037/1‑139 (MQ=255)
                                                                                 cTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGc                                           >  1:572613/1‑139 (MQ=255)
                                                                                 cTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGcc                                          <  1:55882/139‑1 (MQ=255)
                                                                                  tAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGcc                                          <  1:1275330/139‑1 (MQ=255)
                                                                                  tAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGcc                                          <  1:294863/139‑1 (MQ=255)
                                                                                   aGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCt                                         <  1:1804258/139‑1 (MQ=255)
                                                                                   aGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCt                                         <  1:1494865/139‑1 (MQ=255)
                                                                                   aGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCt                                         <  2:572613/139‑1 (MQ=255)
                                                                                      aaaaaCGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGa                                                                          >  2:876939/1‑103 (MQ=255)
                                                                                      aaaaaCGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGa                                                                          <  1:876939/103‑1 (MQ=255)
                                                                                      aaaaaCGCTAAATAGCAGAACATATCGCCGTTAAGCCATCAccccccctct                                                                                                                             >  2:222319/1‑47 (MQ=255)
                                                                                      aaaaaCGCTAAATAGCAGAACATATCGCCGTTAAGCCATCAccccccctct                                                                                                                             <  1:222319/51‑5 (MQ=255)
                                                                                      aaaaaCGCTAAATAGCAGAACATATCGCCGTTAAGCCATCAcc                                                                                                                                     <  1:953147/43‑1 (MQ=255)
                                                                                      aaaaaCGCTAAATAGCAGAACATATCGCCGTTAAGCCATCAcc                                                                                                                                     >  2:953147/1‑43 (MQ=255)
                                                                                          aCGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTg                                                     <  1:641303/120‑1 (MQ=255)
                                                                                          aCGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTg                                                     >  2:641303/1‑120 (MQ=255)
                                                                                           cGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTtata                                 <  1:1031147/139‑1 (MQ=255)
                                                                                              tAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGATAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAAtt                                                 <  2:1119642/120‑1 (MQ=255)
                                                                                              tAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAAtt                                                 >  1:1119642/1‑120 (MQ=255)
                                                                                              tAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCt                                         >  2:1262948/1‑127 (MQ=255)
                                                                                              tAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCt                                         <  1:1262948/127‑1 (MQ=255)
                                                                                                     cAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTccc                                              <  2:1804764/116‑1 (MQ=255)
                                                                                                     cAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTccc                                              >  1:1804764/1‑116 (MQ=255)
                                                                                                     cAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGat                       <  1:342192/139‑1 (MQ=255)
                                                                                                      aGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAAtt                                                 >  1:1024010/1‑112 (MQ=255)
                                                                                                      aGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAAtt                                                 <  2:1024010/112‑1 (MQ=255)
                                                                                                        aaCATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGc                                           <  1:293634/116‑1 (MQ=255)
                                                                                                        aaCATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGc                                           >  2:293634/1‑116 (MQ=255)
                                                                                                          cATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGAttatta                  >  2:982651/1‑139 (MQ=255)
                                                                                                             atCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCgg                                                             <  1:1575852/93‑1 (MQ=255)
                                                                                                             atCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCgg                                                             >  2:1575852/1‑93 (MQ=255)
                                                                                                               cGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGtt                                    >  1:1699850/1‑116 (MQ=255)
                                                                                                               cGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGtt                                    <  2:1699850/116‑1 (MQ=255)
                                                                                                                   gTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGc                                           >  1:953902/1‑105 (MQ=255)
                                                                                                                   gTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGc                                           <  2:953902/105‑1 (MQ=255)
                                                                                                                          cATCACCCCCCCTCTCAATGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGCATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACAAAGGACTTGTCTgg  <  2:1127037/139‑2 (MQ=255)
                                                                                                                               |                                                                                                                                     
AATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACAAAGGACTTGTCTGA  >  NC_000913/1712386‑1712645

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: