Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,926,782 G→A *221* (TAG→TAA
D686D (GAC→GAT
exoX →
ptrB ←
exodeoxyribonuclease 10, DNA exonuclease X
protease II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,926,7820GA100.0% 115.4 / NA 37*221* (TAG→TAA
D686D (GAC→GAT
exoX
ptrB
exodeoxyribonuclease 10, DNA exonuclease X
protease II
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (18/19);  total (18/19)

GGCAAATACCGTGGCAAAGCGGTTTCCGACGTTGCCGAACGCGATCCGGGCTATCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAG  >  NC_000913/1926648‑1926908
                                                                                                                                      |                                                                                                                              
ggCAAATACCGTGGCAAAGCGGTTTCCGACGTTGCCGAACGCGATCCGGGCTATCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCg                                                                                                                            <  1:246305/139‑1 (MQ=255)
                                         cgATCCGGGCTATCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGa                                                                                   <  1:89992/139‑1 (MQ=255)
                                          gATCCGGGCTATCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGaa                                                                                  <  2:1487872/139‑1 (MQ=255)
                                           aTCCGGGCTATCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGaaa                                                                                 >  1:507616/1‑139 (MQ=255)
                                                            tGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTa                                                                <  2:507616/139‑1 (MQ=255)
                                                               ttATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGc                                                             >  2:1256229/1‑139 (MQ=255)
                                                                        aaCCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACg                                                    >  2:9708/1‑139 (MQ=255)
                                                                         aCCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGa                                                   >  2:1660453/1‑139 (MQ=255)
                                                                             ggACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTa                                               >  1:791499/1‑139 (MQ=255)
                                                                                  gCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCg                                          <  1:1660453/139‑1 (MQ=255)
                                                                                   cATGTGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGa                                         <  1:1551588/139‑1 (MQ=255)
                                                                                   cATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTa                                                                                                                    >  2:279070/1‑64 (MQ=255)
                                                                                   cATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTa                                                                                                                    <  1:279070/64‑1 (MQ=255)
                                                                                   cATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGa                                                                                   >  2:864730/1‑97 (MQ=255)
                                                                                   cATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGa                                                                                   <  1:864730/97‑1 (MQ=255)
                                                                                     tGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGAcc                                       >  1:765646/1‑139 (MQ=255)
                                                                                        gCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGa                                    >  2:1359844/1‑139 (MQ=255)
                                                                                        gCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGa                                    >  2:512399/1‑139 (MQ=255)
                                                                                             gAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGAc                                                                                       >  2:1055505/1‑83 (MQ=255)
                                                                                             gAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGAc                                                                                       <  1:1055505/83‑1 (MQ=255)
                                                                                             gAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCt                                                                 >  2:1502365/1‑105 (MQ=255)
                                                                                             gAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCt                                                                 <  1:1502365/105‑1 (MQ=255)
                                                                                             gAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTAc                               >  1:668941/1‑139 (MQ=255)
                                                                                                  gCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCa                          <  1:9708/139‑1 (MQ=255)
                                                                                                   cGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGAc                                                                                       <  1:1549321/77‑1 (MQ=255)
                                                                                                   cGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGAc                                                                                       >  2:1549321/1‑77 (MQ=255)
                                                                                                      ttAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGc                      <  1:1284560/139‑1 (MQ=255)
                                                                                                             tGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAg                                            >  1:927645/1‑110 (MQ=255)
                                                                                                             tGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAg                                            <  2:927645/110‑1 (MQ=255)
                                                                                                             tGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTc               <  1:1367331/139‑1 (MQ=255)
                                                                                                                     tATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCgg       >  2:1401676/1‑139 (MQ=255)
                                                                                                                       tCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTa     <  2:1529422/139‑1 (MQ=255)
                                                                                                                          ggAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACag  <  1:216321/139‑1 (MQ=255)
                                                                                                                             aaaTACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGa                                                                                   >  1:605023/1‑55 (MQ=255)
                                                                                                                             aaaTACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGa                                                                                   <  2:605023/55‑1 (MQ=255)
                                                                                                                                 aCTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGAc                                                                                       <  2:1594230/47‑1 (MQ=255)
                                                                                                                                 aCTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGAc                                                                                       >  1:1594230/1‑47 (MQ=255)
                                                                                                                                      |                                                                                                                              
GGCAAATACCGTGGCAAAGCGGTTTCCGACGTTGCCGAACGCGATCCGGGCTATCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAG  >  NC_000913/1926648‑1926908

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