Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 417,618 G→C 7.2% K69N (AAG→AAC)  waaP → lipopolysaccharide core heptose (I) kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947417,6180GC 7.2% 175.6 / 4.6 56K69N (AAG→AAC) waaPlipopolysaccharide core heptose (I) kinase
Reads supporting (aligned to +/- strand):  ref base G (21/31);  new base C (2/2);  total (23/33)
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.06e-01

ACGCCGTGGAGGCGCTGCAAGGCGAGGTCTATCGCGAGCTGGAAGGGCGCCGCACGCTGCGCACCGAAGTGGCTGGCGAGGGCTTTTTCGTCAAGATCCACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGTGCCGACCATGACCGCAGTCGCCTACGGCGAGCGTGGCAGCAACCCGGCCGCGCAGCACTCGTTCATCATCACCGAAGAGCTGGCGC  >  NC_002947/417470‑417766
                                                                                                                                                    |                                                                                                                                                    
aCGCCGTGGAGGCGCTGCAAGGCGAGGTCTATCGCGAGCTGGAAGGGCGCCGCACGCTGCGCACCGAAGTGGCTGGCGAGGGCTTTTTCGTCAAGATCCACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAg                                                                                                                                                      >  1:310524/1‑149 (MQ=255)
  gCCGTGGAGGCGCTGCAAGGCGAGGTCTATCGCGAGCTGGAAGGGCGCCGCACGCTGCGCACCGAAGTGGCTGGCGAGGGCTTTTTCGTCAAGATCCACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGtt                                                                                                                                                    <  3:123855/149‑1 (MQ=255)
    cGTGGAGGCGCTGCAAGGCGAGGTCTATCGCGAGCTGGAAGGGCGCCGCACGCTGCGCACCGAAGTGGCTGGCGAGGGCTTTTTCGTCAAGATCCACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGGCAACTTgg                                                                                                                                                  >  4:69851/1‑148 (MQ=255)
    cGTGGAGGCGCTGCAAGGCGAGGTCTATCGCGAGCTGGAAGGGCGCCGCACGCTGCGCACCGAAGTGGCTGGCGAGGGCTTTTTCGTCAAGATCCACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGc                                                                                                                                                  <  4:209630/149‑1 (MQ=255)
     gTGGAGGCGCTGCAAGGCGAGGTCTATCGCGAGCTGGAAGGGCGCCGCACGCTGCGCACCGAAGTGGCTGGCGAGGGCTTTTTCGTCAAGATCCACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGcc                                                                                                                                                 <  3:200127/149‑1 (MQ=255)
      tGGAGGCGCTGCAAGGCGAGGTCTATCGCGAGCTGGAAGGGCGCCGCACGCTGCGCACCGAAGTGGCTGGGGAGGGCTTTTTCGTCAAGATCCACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCg                                                                                                                                                <  6:295545/149‑1 (MQ=255)
       ggAGGCGCTGCAAGGCGAGGTCTATCGCGAGCTGGAAGGGCGCCGCACGCTGCGCACCGAAGTGGCTGGCGAGGGCTTTTTCGTCAAGATCCACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCgg                                                                                                                                               >  7:228013/1‑149 (MQ=255)
             gcTGCAAGGCGAGGTCTATCGCGAGCTGGAAGGGCGCCGCACGCTGCGCACCGAAGTGGCTGGCGAGGGCTTTTTCGTCAAGATCCACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGGCAACTTGGCGGTgaaca                                                                                                                                         <  3:69851/149‑5 (MQ=255)
              cTGCAAGGCGAGGTCTATCGCGAGCTGGAAGGGCGCCGCACGCTGCGCACCGAAGTGGCTGGCGAGGGCTTTTTCGTCAAGATCCACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCgg                                                                                                                                        >  5:31937/1‑149 (MQ=255)
              cTGCAAGGCGAGGTCTATCGCGAGCTGGAAGGGCGCCGCACGCTGCGCACCGAAGTGGCTGGCGAGGGCTTTTTCGTCAAGATCCACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCgg                                                                                                                                        <  4:326450/149‑1 (MQ=255)
              cTGCAAGGCGAGGTCTATCGCGAGCTGGAAGGGCGCCGCACGCTGCGCACCGAAGTGGCTGGCGAGGGCTTTTTCGTCAAGATCCACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCgg                                                                                                                                        <  2:61475/149‑1 (MQ=255)
              cTGCAAGGCGAGGGCTATCGCGAGCTGGAAGGGCGCCGCACGCTGCGCACCGAAGGGGCTGGCGAGGGCTTTTTCGTCAAGATCCACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCgg                                                                                                                                        <  4:126951/149‑1 (MQ=255)
               tGCAAGGCGAGGTGTATCTCGAGCTGGAAGGGCGCCTCCCGCTGCGCACCGAAGTGGCTGGCGAGGGCTTTTTCGTCAAGATCCACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCggg                                                                                                                                       <  8:266013/149‑1 (MQ=255)
               tGCAAGGCGAGGTCTATCGCGAGCTGGAAGGGCGCCGCACGCTGCGCACCGAAGTGGCTGGCGAGGGCTTTTTCGTCAAGATCCACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCggg                                                                                                                                       >  3:280821/1‑149 (MQ=255)
                 cAAGGCGAGGTCTATCGCGAGCTGGAAGGGCGCCGCACGCTGCGCACCGAAGTGGCTGGCGAGGGCTTTTTCGTCAAGATCCACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGc                                                                                                                                     >  8:277530/1‑149 (MQ=255)
                 cAAGGCGAGGTCTATCGCGAGCTGGAAGGGCACCGCACGCTGCGCACCGAAGTGGCTGGCGAGGGCTTTTTCGTCAAGATCCACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCg                                                                                                                                                >  3:333783/1‑138 (MQ=255)
                 cAAGGCGAGGTCTATCGCGAGCTGGAAGGGCACCGCACGCTGCGCACCGAAGTGGCTGGCGAGGGCTTTTTCGTCAAGATCCACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCg                                                                                                                                                <  4:333783/138‑1 (MQ=255)
                  aaGGCGAGGGCTATCGCGAGCTGGAAGGGCGCCGCACGCTGCGCACCGAAGTGGCTGGCGAGGGCTTTTTCGTCAAGATCCACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCGCGGGGCt                                                                                                                                    <  2:140425/149‑1 (MQ=255)
                        atgTCTATAGCGAGCTGTAAGGGCGCCGCACGCTGCGCACCGAAGTGGCTGGCGAGGGCTTTTTCGGCAAGATCCACCGTGGCATCGGCTGGGGCGAGATCTTCATGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGcctg                                                                                                                              <  4:32583/147‑3 (MQ=255)
                           tCTATCGCGAGCTGGAAGGGCGCCGCACGCTGCGCACCGAAGTGGCTGGCGAGGGCTTTTTCGTCAAGATCCACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGaa                                                                                                                           <  7:229372/149‑1 (MQ=255)
                               tCGCGAGCTGGAAGGGCGCCGCACGCTGCGCACCGAAGTGGCTGGCGAGGGCTTTTTCGTCAAGATCCACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATggc                                                                                                                       <  5:286483/149‑1 (MQ=255)
                               tCGCGAGCTGGAAGGGCGCCGCACGCTGCGCACCGAAGTGGCTGGCGAGGGCTTTTTCGTCAAGATCCACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATgg                                                                                                                        <  6:40174/148‑1 (MQ=255)
                                 gcgAGCTGGAAGGGCGCCGCACGCTGCGCACCGAAGTGGCTGGCGAGGGCTTTTTCGTCAAGATCCACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATggcg                                                                                                                      <  6:31937/148‑1 (MQ=255)
                                  cgAGCTGGAAGGGCGCCGCACGCTGCGCACCGAAGTGGCTGGCGAGGGCTTTTTCGTCAAGATCCACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATggcggg                                                                                                                    <  3:318727/149‑1 (MQ=255)
                                  cgAGCTGGAAGGGCGCCGCACGCTGCGCACCGAAGTGGCTGGCGAGGGCTTTTTCGTCAAGATCCACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATggcggg                                                                                                                    >  1:231932/1‑149 (MQ=255)
                                   gAGCTGGAAGGGCGCCGCACGCTGCGCACCGAAGTGGCTGGCGAGGGCTTTTTCGTCAAGATCCACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATggcgggc                                                                                                                   <  4:229514/149‑1 (MQ=255)
                                       tGGAAGGGCGCCGCACGCTGCGCACCGAAGTGGCTGGCGAGGGCTTTTTCGTCAAGATCCACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATc                                                                                                               >  8:124433/1‑149 (MQ=255)
                                       tGGAAGGGCGCCGCACGCTGCGCACCGAAGTGGCTGGCGAGGGCTTTTTCATCAAGATCCACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATc                                                                                                               >  1:92538/1‑149 (MQ=255)
                                         gAAGTGCGCCGCACGCTGCGCACCGAAGTGGCTGGCGAGGGCTTTTTCGTCAAGATCCAGCGTGGCATCGGGTGGGGCGAGATCTTCAAGAACCTGTTCTCCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATc                                                                                                               <  8:248728/147‑1 (MQ=255)
                                         gAAGGGCGCCGCACGCTGCGCACCGAAGTGGCTGGCGAGGGCTTTTTCGTCAAGATCCACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATc                                                                                                               <  4:280821/147‑1 (MQ=255)
                                         gAAGGGCGCCGCACGCTGCGCACCGAAGTGGCTGGCGAGGGCTTTTTCGTCAAGATCCACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATc                                                                                                               <  2:231932/147‑1 (MQ=255)
                                         gAAGGGCGCCGCACGCTGCGCACCGAAGTGGCTGGCGAGGGCTTTTTCGTCAAGATCCACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCGCGGGGCTGGCCAGGAATGGCGGGCGATc                                                                                                               <  2:92538/147‑1 (MQ=255)
                                                              aCCGAAGTGGCTGGCGAGGGCTTTTTCGTCAAGATCCACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGt                                                                                        >  3:45841/1‑149 (MQ=255)
                                                                 gAAGTGGCTGGCGAGGGCTTTTTCGTCAAGATCCACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCa                                                                                                    <  7:31405/134‑1 (MQ=255)
                                                                 gAAGTGGCTGGCGAGGGCTTTTTCGTCAAGATCCACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGTGcc                                                                                     <  7:124433/149‑1 (MQ=255)
                                                                         tggcGAGGGCTTTTTCGTCAAGATCCACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGTGCCGACCATGa                                                                             >  5:165678/1‑149 (MQ=255)
                                                                              aGGGCTTTTTCGTCAAGATCCACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGTGCCGACCATGACCGCa                                                                        <  7:277530/149‑1 (MQ=255)
                                                                              aGGGCTTTTTCGTCAAGATCCACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGTACCGACCATGACCg                                                                          >  1:102218/1‑147 (MQ=255)
                                                                                         gTCAAGATCCACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGTGCCGACCATGACCGCAGTCGCCTACgg                                                             >  2:150754/1‑149 (MQ=255)
                                                                                          tCAAGATCCACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGTGCCGACCATGACCGCAGTCGCCTACGGc                                                            >  6:161656/1‑149 (MQ=255)
                                                                                            aaGATCCACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGTGCCGACCATGACCGCAGTCGCCTACGGCGa                                                          <  2:310524/149‑1 (MQ=255)
                                                                                                 ccACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGTGCCGACCATGACCGCAGTCGCCTACGGCGAGCGTg                                                     >  8:198868/1‑149 (MQ=255)
                                                                                                  cACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGTGCCGACCATGACCGCAGTCGCCTACGGCGAGCGTgg                                                    <  6:165678/149‑1 (MQ=255)
                                                                                                      gTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGTGCCGACCATGACCGCAGTCGCCTACGGCGAGCGTGgcag                                                 >  6:176133/1‑148 (MQ=255)
                                                                                                                   gggCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGTGCCGACCATGACCGCAGTCGCCTACGGCGAGCGTGGCAGCAACCCGGCcgcg                                    <  4:45841/148‑1 (MQ=255)
                                                                                                                             ttCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGTGCCGACCATGACCGCAGTCGCCTACGGCGAGCGTGGCAGCAACCCGGCCGCGCAGCACTCGtt                         >  6:237855/1‑149 (MQ=255)
                                                                                                                             ttCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGTGCCGACCATGACCGCAGTCGCCTACGGCGAGCGTGGCAGCAACCCGGCCGCGCAGCACTCGt                          >  4:142432/1‑148 (MQ=255)
                                                                                                                               cAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGTGCCGACCATGACCGCAGTCGCCTACGGCGAGCGTGGCAGCAACCCGGCCGCGCAGCACTc                            <  4:132720/144‑1 (MQ=255)
                                                                                                                               cAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGTGCCGACCATGACCGCAGTCGCCTACGGCGAGCGTGGCAGCAACCCGGCCGCGCAGCACTCGTtc                        <  5:161656/148‑1 (MQ=255)
                                                                                                                                aaGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGTGCCGACCATGACCGCAGTCGCCTACGGCGAGCGTGGCAGCAACCCGGCCGCGCAGCACTCGTtcat                      >  3:145055/1‑149 (MQ=255)
                                                                                                                                  gAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGTGCCGACCATGACCGCAGTCGCCTACGGCGAGCGTGGCAGCAACCCGGCCGCGCAGCACTCGTtcatca                    >  1:227607/1‑149 (MQ=255)
                                                                                                                                      cTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCTGGCGATCCAACGCCTGCATGAAGTGGGCGTGCCGACCATGACCGCAGTCGCCTACGGCGAGCGTGGCAGCAACCCGGCCGCGCAGCACTCGTTCATCATCAc                <  4:298443/149‑1 (MQ=255)
                                                                                                                                      cTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGTACCGACCATGACCGCAGTCGCCTACGGCGAGCGTGGCAGCAACCCGGCCGCGCAGCACTCGTTCATCATCAc                <  2:102218/149‑1 (MQ=255)
                                                                                                                                       cgagcACCGGCAACTTGGCGGTGCTCGGGGCAGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGTGCCGACCATGACCGCAGTCGCCTACGGCGAGCGTGGCAGCAACCCGGCCGCGCAGCACTCGTTCATCATCAcc               >  3:334144/5‑149 (MQ=255)
                                                                                                                                            aCCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGTGCCGACCATGACCGCAGTCGCCTACGGCGAGCGTGGCAGCAACCCGGCCGCGCAGCACTCGTTCATCATCACCGaa            >  8:152208/1‑147 (MQ=255)
                                                                                                                                                gcAACTTGGCGGTGCTCGGGGCAGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGTGCCGACCATGACCGCAGTCGCCTACGGCGAGCGTGGCAGCAACCCGGCCGCGCAGCACTCGTTCATCATCACCGAAGAGCTg      <  4:334144/148‑1 (MQ=255)
                                                                                                                                                 cAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGTGCCGACCATGACCGCAGTCGCCTACGGCGAGCGTGGCAGCAACCCGGCCGCGCAGCACTCGTTCATCATCACCGAAGAGCTgg     >  7:264732/1‑149 (MQ=255)
                                                                                                                                                    gTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGTGCCGACCATGACCGCAGTCGCCTACGGCGAGCGTGGCAGCAACCCGGCCGCGCAGCACTCGTTCATCATCACCGAAGAGCTGgcgc  <  7:198868/149‑1 (MQ=255)
                                                                                                                                                    gTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGTGCCGACCATGACCGCAGTCGCCTACGGCGAGCGTGGCAGCAACCCGGCCGCGCAGCACTCGTTCATCATCACCGAAGAGCTGgcgc  >  8:177993/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ACGCCGTGGAGGCGCTGCAAGGCGAGGTCTATCGCGAGCTGGAAGGGCGCCGCACGCTGCGCACCGAAGTGGCTGGCGAGGGCTTTTTCGTCAAGATCCACCGTGGCATCGGCTGGGGCGAGATCTTCAAGAACCTGTTCACCGCCAAGTTGCCGGTGCTCGGGGCTGGCCAGGAATGGCGGGCGATCCAACGCCTGCATGAAGTGGGCGTGCCGACCATGACCGCAGTCGCCTACGGCGAGCGTGGCAGCAACCCGGCCGCGCAGCACTCGTTCATCATCACCGAAGAGCTGGCGC  >  NC_002947/417470‑417766

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: ‑