Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,541,711 T→G 7.7% I330L (ATC→CTC)  dnaB ← replicative DNA helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,541,7110TG 7.7% 198.0 / 5.6 64I330L (ATC→CTC) dnaBreplicative DNA helicase
Reads supporting (aligned to +/- strand):  ref base T (31/28);  new base G (3/2);  total (34/30)
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.89e-01

AGCGGTTCAGCTGCGAAAGGGCAATGACCGGGCAGTTGAATTCTTTGGCCAGGGCC‑TTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAGTCGTCGTCGTC  >  NC_002947/5541564‑5541858
                                                                                                                                                    |                                                                                                                                                   
agcgGTTCAGCTGCGAAAGGGCAATGACCGGGCAGTTGAATTCTTTGGCCAGGGCC‑TTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGAt                                                                                                                                                     >  2:277498/1‑148 (MQ=255)
agcgGTTCAGCTGCGAAAGGGCAATGACCGGGCAGTTGAATTCTTTGGCCAGGGCC‑TTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGAt                                                                                                                                                     >  6:161534/1‑148 (MQ=255)
agcgGTTCAGCTGCGAAAGGGCAATGACCGGGCAGTTGAATTCTTTGGCCAGGGCC‑TTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATc                                                                                                                                                    <  1:325572/149‑1 (MQ=255)
           tGCGAAAGGGCAATGACCGGGCAGTTGAATTCTTTGGCCAGGGCC‑TTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGAGCATGGCGaaatc                                                                                                                                         >  1:234359/1‑145 (MQ=255)
               caacgggaaTGTCCGGGCACTTGAATTATTTTGCCATGGCCTTTCTGCATGCGGTAAAGCCCGTAAATTTGGTTGGT‑CGGGTGGCAACGGGCGCACCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCa                                                                                                                                     <  3:345717/142‑1 (MQ=255)
                aaGGGCAATGACCGGGCAGTTTAATTCTTTGGCCAGGGCA‑GTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGAAATCGCCAt                                                                                                                                    <  4:531/149‑1 (MQ=255)
                  gggCAATGACCGGGCAGTGGAATTCTTTGGCCAGGGCC‑TTGAGCGAGCGGGAAATCTCGGAAATTTCGGTGGTACGGTTGTAACCTGCAGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGc                                                                                                                                  <  8:161761/149‑1 (MQ=255)
                           ccGGGCAGTTGAATTCTTTGGCCAGGGCC‑TTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGAGCATGGCGAAATCGCCATGATCATggtc                                                                                                                          <  2:234359/148‑3 (MQ=255)
                                  gTTGAATTCTTTGGCCAGGGCC‑TTGAGCGAGCGGGAAATCTCTGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGCAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAggcggc                                                                                                                  <  2:115086/149‑1 (MQ=255)
                                    tGAATTCTTTGGCCAGGGCC‑TTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACTGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGcggcgg                                                                                                                >  4:234266/1‑149 (MQ=255)
                                     gAATTCTTTGGCCAGGGCC‑TTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCggg                                                                                                               >  1:354076/1‑149 (MQ=255)
                                     gAATTCTTTGGCCAGGGCC‑TTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCggg                                                                                                               >  3:318192/1‑149 (MQ=255)
                                     gAATTCTTTGGCCAGGGCC‑TTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCggg                                                                                                               >  4:11774/1‑149 (MQ=255)
                                     gAATTCTTTGGCCAGGGCC‑TTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGGACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTTCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGcggcgg                                                                                                                <  1:99617/148‑1 (MQ=255)
                                           tttGGCCAGGGCC‑TTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTAcgcg                                                                                                         >  8:176889/1‑149 (MQ=255)
                                                ccAGGGCC‑TTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGcggcgg                                                                                                                <  8:92847/137‑1 (MQ=255)
                                                ccAGGGCC‑TTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGcggcgg                                                                                                                >  7:92847/1‑137 (MQ=255)
                                                       c‑TTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCTCCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTcc                                                                                              >  1:259874/1‑148 (MQ=255)
                                                         ttGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTcc                                                                                              >  3:263543/1‑147 (MQ=255)
                                                           gagcgagcGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAg                                                                                           <  6:47976/148‑1 (MQ=255)
                                                           gagcgagcGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCACGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGa                                                                                            <  6:348325/147‑1 (MQ=255)
                                                                agcgGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTg                                                                                     <  3:11774/149‑1 (MQ=255)
                                                                 gcgGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGa                                                                                    <  2:354076/149‑1 (MQ=255)
                                                                      aaaTCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATAccc                                                                               <  2:140699/149‑1 (MQ=255)
                                                                            cGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTa                                                                         >  2:258153/1‑149 (MQ=255)
                                                                             ggAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACccgcc                                                                            <  2:259874/145‑1 (MQ=255)
                                                                             ggAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGt                                                                          <  4:263543/147‑1 (MQ=255)
                                                                               aaaTTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTAtcg                                                                      >  7:354451/1‑149 (MQ=255)
                                                                                    tCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGa                                                                  <  1:277498/148‑1 (MQ=255)
                                                                                    tCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGa                                                                  <  5:161534/148‑1 (MQ=255)
                                                                                    tCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGAt                                                                 >  7:267849/1‑149 (MQ=255)
                                                                                      gTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATg                                                                >  2:348574/1‑148 (MQ=255)
                                                                                      gTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATg                                                                >  1:46319/1‑148 (MQ=255)
                                                                                       ttGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGTTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGaa                                                              <  3:234266/149‑1 (MQ=255)
                                                                                            aCGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAAtt                                                         >  8:70839/1‑149 (MQ=255)
                                                                                               tttGTCACCGGCGGACCCCGGAATCTGCACCATCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTc                                                      <  5:330728/148‑1 (MQ=255)
                                                                                               gTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTc                                                      >  8:75380/1‑149 (MQ=255)
                                                                                                ttGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCg                                                     <  1:258153/149‑1 (MQ=255)
                                                                                                ttGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCg                                                     <  7:75380/149‑1 (MQ=255)
                                                                                                     aCCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGt                                                >  3:68551/1‑149 (MQ=255)
                                                                                                       cGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTg                                              <  8:74180/149‑1 (MQ=255)
                                                                                                        ggCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAgg                                                                                          >  6:168815/1‑104 (MQ=255)
                                                                                                         gCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATTGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAg                                            <  8:354451/149‑1 (MQ=255)
                                                                                                         cccGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAg                                            >  8:17335/2‑149 (MQ=255)
                                                                                                            gACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAg                                         >  3:159598/1‑149 (MQ=255)
                                                                                                            gACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTAGGTCGTTGAGCAg                                         >  8:282006/1‑149 (MQ=255)
                                                                                                              ccccGGAATCTGCATCAGCTGCAGGTAGTCGACGATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGt                                       <  1:346240/149‑1 (MQ=255)
                                                                                                                 cGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTAGGTCGTTGAGCAGGTTCa                                    >  1:103249/1‑149 (MQ=255)
                                                                                                                         gCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAg                            >  8:198691/1‑149 (MQ=255)
                                                                                                                            tCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCAAGCACACCGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGAAGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAgct                          >  4:20072/1‑146 (MQ=255)
                                                                                                                                cTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGgc                     <  7:70839/149‑1 (MQ=255)
                                                                                                                                   cAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGcgcg                  <  7:282006/149‑1 (MQ=255)
                                                                                                                                   cAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGcgcg                  <  8:267849/149‑1 (MQ=255)
                                                                                                                                       tAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCa              <  2:46319/149‑1 (MQ=255)
                                                                                                                                          tCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAgtc           >  5:298085/1‑149 (MQ=255)
                                                                                                                                          tCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAgtc           >  4:81000/1‑149 (MQ=255)
                                                                                                                                          tCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAgtc           >  2:164804/1‑149 (MQ=255)
                                                                                                                                            gACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAgtcgt         >  3:183940/1‑149 (MQ=255)
                                                                                                                                            gACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAgtcgt         >  4:162159/1‑149 (MQ=255)
                                                                                                                                             aCCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAgtcgtc        >  3:245753/1‑149 (MQ=255)
                                                                                                                                               cgtGAGCATGGCGAAATCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAAc                                                             >  2:137947/3‑94 (MQ=39)
                                                                                                                                               cgtGAGCATGGCGAAATCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAAc                                                             <  1:137947/92‑1 (MQ=39)
                                                                                                                                               cgtGAGCATGGCGAAATCGCCATGATCATGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAgtcgtcgt      >  3:79955/3‑149 (MQ=255)
                                                                                                                                                aTGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAgtcgtcgtc     >  6:52589/1‑149 (MQ=255)
                                                                                                                                                  gATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAgtc           >  3:23758/1‑141 (MQ=255)
                                                                                                                                                  gATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAgtc           <  4:23758/141‑1 (MQ=255)
                                                                                                                                                    tCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAgtcgtcgtcgtc  >  7:165594/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AGCGGTTCAGCTGCGAAAGGGCAATGACCGGGCAGTTGAATTCTTTGGCCAGGGCC‑TTGAGCGAGCGGGAAATCTCGGAAATTTCGTTGGTACGGTTGTCACCGGCCGACCCCGGAATCTGCATCAGCTGCAGGTAGTCGACCATGATCATGGCGATTTCGCCATGCTCACGGGCCAGGCGGCGGGTACGCGCGCGCATTTCCGAGGGGCTGATACCCGCCGTATCGTCGATGAACAATTTTCGGTCGTTGAGCAGGTTCACTGCCGAGGTCAGGCGCGGCCAGTCGTCGTCGTC  >  NC_002947/5541564‑5541858

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑