Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 972,232 G→C 9.7% G175A (GGC→GCC)  secD → protein translocase subunit SecD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947972,2320GC 9.7% 178.0 / 9.7 62G175A (GGC→GCC) secDprotein translocase subunit SecD
Reads supporting (aligned to +/- strand):  ref base G (31/25);  new base C (3/3);  total (34/28)
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.61e-01

GACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGATTTCGACCTGACCACCACCGATCGCAACGAGCTTGCGGTGCTGCGTCTGGCGCTGACCCAGGCGA  >  NC_002947/972090‑972372
                                                                                                                                              |                                                                                                                                            
gACCTGTCCGGAGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAAGATGGCGGcat                                                                                                                                        >  1:93102/1‑149 (MQ=255)
      tCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGc                                                                                                                                   <  8:29877/148‑1 (MQ=255)
       ccGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCATGATGCCGCCATCctgctg                                                                                                                                 <  4:73608/149‑1 (MQ=255)
        cGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTg                                                                                                                                 <  5:129036/148‑1 (MQ=255)
            ggtgTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGTCATCATGc                                                                                                                                   <  2:147538/142‑1 (MQ=255)
            ggtgTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGc                                                                                                                                   >  1:147538/1‑142 (MQ=255)
              tgtgCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGAGCGTCCGCTACAGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCg                                                                                                                          >  5:105707/1‑149 (MQ=255)
               gtgCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGc                                                                                                                         >  8:239346/1‑149 (MQ=255)
               gtgCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGc                                                                                                                         >  2:152861/1‑149 (MQ=255)
                tgCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCt                                                                                                                        >  3:302369/1‑149 (MQ=255)
                   gcTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATggcg                                                                                                                                            <  8:232838/125‑1 (MQ=255)
                   aCTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGAc                                                                                                                     <  1:152861/149‑1 (MQ=255)
                   aCTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGAACTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGAc                                                                                                                     >  7:314952/1‑149 (MQ=255)
                   aCTTCCATCTGTAAGTGGACATTGACAAGGACATGACAGCCCGCATGAAAGTCTACTAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGAc                                                                                                                     <  2:264133/149‑1 (MQ=255)
                     ttCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTc                                                                                                                           <  1:163197/141‑1 (MQ=255)
                     ttCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTc                                                                                                                           >  2:163197/1‑141 (MQ=255)
                       ccTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAAGATGGCGGCATCATGCTGGGCTTCGCTGACGAt                                                                                                                  <  2:93102/148‑1 (MQ=255)
                           ctgGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCa                                                                                                              <  6:275027/148‑1 (MQ=255)
                                   ggacatggacaAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGaa                                                                                                      >  4:155878/1‑148 (MQ=255)
                                    gacatggacaAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGaa                                                                                                      >  4:322892/1‑147 (MQ=255)
                                     acatggacaAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAg                                                                                                   >  7:198827/1‑149 (MQ=255)
                                      catggacaAGTCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCATGATGCCGCCATCctgctga                                                                                                                                <  6:322478/119‑2 (MQ=255)
                                      catggacaAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCATGATGCCGCCATCctgctga                                                                                                                                >  5:322478/1‑118 (MQ=255)
                                         ggacaAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCAc                                                                                               >  8:116656/1‑149 (MQ=255)
                                               ggCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGccc                                                                                         >  5:169431/1‑149 (MQ=255)
                                                 ccATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTg                                                                                                                                 <  6:194139/107‑1 (MQ=255)
                                                 ccATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTg                                                                                                                                 >  5:194139/1‑107 (MQ=255)
                                                        tGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGATGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCACCCAGCGAACAGGCACATGCCCTGATCCGc                                                                                 >  4:253388/1‑148 (MQ=255)
                                                           ccGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAg                                                                              >  6:113413/1‑148 (MQ=255)
                                                           ccGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAg                                                                              >  6:353346/1‑148 (MQ=255)
                                                                tGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTc                                                                        <  7:116656/149‑1 (MQ=255)
                                                                tGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTc                                                                        <  5:113413/149‑1 (MQ=255)
                                                                   aaGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCaa                                                                      <  3:47820/148‑1 (MQ=255)
                                                                   aaGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAAt                                                                     <  3:139177/149‑1 (MQ=255)
                                                                     gTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGa                                                                   >  4:192667/1‑149 (MQ=255)
                                                                         aCGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCCGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGATTTc                                                               >  5:91983/1‑149 (MQ=255)
                                                                          cGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGATTTc                                                               <  6:169431/148‑1 (MQ=255)
                                                                            aaGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGATTTCGAc                                                            <  7:239346/149‑1 (MQ=255)
                                                                            aaGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGATTTCGAc                                                            <  3:86335/149‑1 (MQ=255)
                                                                              ggCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGATTTCGACCt                                                          <  5:353346/149‑1 (MQ=255)
                                                                                cGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTAGCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCTCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGATTTCGACCTGa                                                        <  6:105707/149‑1 (MQ=255)
                                                                                 gAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGATTTCGACCTGac                                                       <  3:322892/149‑1 (MQ=255)
                                                                                   ggTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGATTTCGACCTGacca                                                     >  4:103967/1‑149 (MQ=255)
                                                                                      cAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGAGCCCTGATCCGCAAGAATTTCAATGATTTCGACCTCaccacca                                                  >  5:98414/1‑149 (MQ=255)
                                                                                         gACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACGGGCACGTGCCCTGATCCGCAAGAATTTCAATGATTTCGACCTGACCACCACCg                                               >  6:128824/1‑149 (MQ=255)
                                                                                         gACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGATTTCGACCTGaccaccacc                                                >  4:299414/1‑148 (MQ=255)
                                                                                           ccTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGATTTCGACCTGACCACCACCGa                                              <  3:155878/148‑1 (MQ=255)
                                                                                              tgctgcGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGATTTCGACCTGACCACCACCGATCg                                           >  7:111752/1‑148 (MQ=255)
                                                                                                    gcAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGATTTCGACCTGACCACCACCGATCGCAACGAg                                    >  8:101923/1‑149 (MQ=255)
                                                                                                     cAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGATTTCGACCTGACCACCACCGATCGCAAGAGCt                                  >  2:69379/1‑149 (MQ=255)
                                                                                                                  cGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGATTTCGACCTGACCACCACCGATCGCAACGAGCTTGCGGtgctgc                       <  7:101923/148‑1 (MQ=255)
                                                                                                                  cGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGATTTCGACCTGACCACCACCGATCGCAACGAGCTTGCGGtgctgc                       <  3:192667/148‑1 (MQ=255)
                                                                                                                          cAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGATTTCGACCTGACCACCACCGATCGCAACGAGCTTGCGGTGCTGCGTCTGgcgc              >  6:123846/1‑149 (MQ=255)
                                                                                                                                gcctCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGATTTCGACCTGACCACCACCGATCGCAACGAGCTTGCGGTGCTGCGTCTGGCGCTGAcc         >  6:262348/1‑148 (MQ=255)
                                                                                                                                  ctCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGATTTCGACCTGACCACCACCGATCGCAACGAGCTTGCGGTGCTGCGTCTGGCGCTGACCCAg      <  8:198827/149‑1 (MQ=255)
                                                                                                                                  cccagcatGATGCCGCCATCCTGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGATTTCGACCTGACCACCACCGATCGCAACGAGCTTGCGGTGCTGCGTCt                   >  2:198004/3‑136 (MQ=255)
                                                                                                                                  cccagcatGATGCCGCCATCCTGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGATTTCGACCTGACCACCACCGATCGCAACGAGCTTGCGGTGCTGCGTCt                   <  1:198004/134‑1 (MQ=255)
                                                                                                                                  cccagcatGATGCCGCCATCCTGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGATTTCGACCTGACCACCACCGATCGCAACGAGCTTGCGGTGCTGCGTCTGGCGCTGACCCAg      >  8:214861/3‑149 (MQ=255)
                                                                                                                                   tCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGATTTCGACCTGACCACCACCGATCGCAACGAGCTTGCGGtgct                         >  5:163058/1‑129 (MQ=255)
                                                                                                                                   tCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGATTTCGACCTGACCACCACCGATCGCAACGAGCTTGCGGtgct                         <  6:163058/129‑1 (MQ=255)
                                                                                                                                    cagcagGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGATTTCGACCTGACCACCACCGATCGCAACGAGCTTGCGGTGCTGCGTCTGGCGCTGACCCAGGc    <  8:111752/149‑1 (MQ=255)
                                                                                                                                      gcagGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGATTTCGACCTGACCACCACCGATCGCAACGAGCTTGCGGTGCTGCGTCTGGCGCTGACCCAGGCg   >  7:217950/1‑148 (MQ=255)
                                                                                                                                        agGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGATTTCGACCTGACCACCACCGATCGCAACGAGCTTGCGGTGCTGCGTCTGGCGCTGACCCAGGCGa  >  3:72458/1‑147 (MQ=255)
                                                                                                                                              |                                                                                                                                            
GACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGATTTCGACCTGACCACCACCGATCGCAACGAGCTTGCGGTGCTGCGTCTGGCGCTGACCCAGGCGA  >  NC_002947/972090‑972372

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: