Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 972,235 G→C 9.7% G176A (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,2350GC 9.7% 169.6 / 9.7 62G176A (GGC→GCC) secDprotein translocase subunit SecD
Reads supporting (aligned to +/- strand):  ref base G (31/24);  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 = 8.91e-02

GACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGATTTCGACCTGACCACCACCGATCGCAACGAGCTTGCGGTGCTGCGTCTGGCGCTGACCCAGGCGAAAGTCGCCGAG  >  NC_002947/972090‑972383
                                                                                                                                                 |                                                                                                                                                    
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                                                                                                                                    >  2:152861/1‑149 (MQ=255)
               gtgCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGc                                                                                                                                    >  8:239346/1‑149 (MQ=255)
                tgCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCt                                                                                                                                   >  3:302369/1‑149 (MQ=255)
                   aCTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGAc                                                                                                                                <  1:152861/149‑1 (MQ=255)
                   aCTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGAACTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGAc                                                                                                                                >  7:314952/1‑149 (MQ=255)
                   aCTTCCATCTGTAAGTGGACATTGACAAGGACATGACAGCCCGCATGAAAGTCTACTAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGAc                                                                                                                                <  2:264133/149‑1 (MQ=255)
                     ttCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTc                                                                                                                                      >  2:163197/1‑141 (MQ=255)
                     ttCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTc                                                                                                                                      <  1:163197/141‑1 (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:353346/1‑148 (MQ=255)
                                                           ccGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAg                                                                                         >  6:113413/1‑148 (MQ=255)
                                                                tGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTc                                                                                   <  5:113413/149‑1 (MQ=255)
                                                                tGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTc                                                                                   <  7:116656/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                                  <  3:192667/148‑1 (MQ=255)
                                                                                                                  cGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGATTTCGACCTGACCACCACCGATCGCAACGAGCTTGCGGtgctgc                                  <  7:101923/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                              <  1:198004/134‑1 (MQ=255)
                                                                                                                                  cccagcatGATGCCGCCATCCTGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGATTTCGACCTGACCACCACCGATCGCAACGAGCTTGCGGTGCTGCGTCt                              >  2:198004/3‑136 (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)
                                                                                                                                                ggcATGCTGGTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCTCGTGCCCTGATCCGCAAGAATATCAATGATTTCGACCTGACCACCACCGATCGCACCGAGCTTGCGGTGCTCCGTCTGGCGCTGACCCAGGCGAAAGGCGCCga   <  8:314952/149‑1 (MQ=255)
                                                                                                                                                 gcAACATGCTTGTCTTCGCTGACGAGGCCACTCGCTACTAGGCCCGCGCCCTGATCCGCAAGAATTTCAATGATTTCGACCTGACCACCACCGATCGCAACGAGCTTGCGGTGCTGCGTCTGGCGCTGACCCAGGCGAAAGTCGCCgag  <  6:292365/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCCGCGAACAGGCACGTGCCCTGATCCGCAAGAATTTCAATGATTTCGACCTGACCACCACCGATCGCAACGAGCTTGCGGTGCTGCGTCTGGCGCTGACCCAGGCGAAAGTCGCCGAG  >  NC_002947/972090‑972383

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