Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 416,121 G→C 5.9% D297H (GAC→CAC)  waaC → ADP‑heptose‑‑LPS heptosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947416,1210GC 5.9% 319.2 / 6.7 98D297H (GAC→CAC) waaCADP‑heptose‑‑LPS heptosyltransferase
Reads supporting (aligned to +/- strand):  ref base G (46/46);  new base C (2/4);  total (48/50)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.21e-01

GCGGCGGCAAAGGCCTGCGTGGCGGTCGATACCGGCCTTGGCCACCTGGCGGCGGCACTCGATGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCCTACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAACGCGAGTGGCCGCTGTGCTTCACTCGCCTGAATCCCGAGCATGTGGCGAGCCGCTTGAGCGCGT  >  NC_002947/415977‑416265
                                                                                                                                                 |                                                                                                                                                
gcggcggcAAAGGCCTGCGTGGCGGTCGATACCGGCCTTGGGCACCTGGCGGCGGCACTCGATGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCCTACGGACGGACCCAGATTCACCAGGCCAGTGACTg                                                                                                                                              <  5:23691/149‑1 (MQ=255)
  ggcggcAAAGGCCTGCGTGGCGGTCGATACCGGCCTTGGCCACCTGGCGGCGGCACTCGATGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCCTACGGACGGACCCAGATTCACCAGGCCAGTGACTGGc                                                                                                                                            <  2:363289/149‑1 (MQ=255)
   gcggcAAAGGCCTGCGTGGCGGTCGATACCGGCCTTGGCCACCTGGCGGCGGCACTCGATGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCCTACGGACGGACCCAGATTCACCAGGCCAGTGACTGGcc                                                                                                                                           <  3:292449/149‑1 (MQ=255)
   gcggcAAAGGCCTGCGTGGCGGTCGATACCGGCCTTGGCCACCTGGCGGCGGCACTCGATGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCCTACGGACGGACCCAGATTCACCAGGCCAGTGACTGGcc                                                                                                                                           >  5:8995/1‑149 (MQ=255)
   gcggcAAAGGCCTGCGTGGCGGTCGATACCGGCCTTGGCCACCTGGCGGCGGCACTCGATGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCCTACGGACGGACCCAGATTCACCAGGCCAGTCACTGGcc                                                                                                                                           >  6:193767/1‑149 (MQ=255)
     ggcAAAGGCCTGCGTGGCGGTCGATACCGGCCTTGGCCACCTGGCGGCGGCACTCGATGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCCTACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCtt                                                                                                                                         <  2:48008/149‑1 (MQ=255)
     ggcAAAGGCCTGCGTGGCGGTCGATACCGGCCTTGGCCACCTGGCGGCGGCACTCGATGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCCTACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCtt                                                                                                                                         >  5:212247/1‑149 (MQ=255)
      gcAAAGGCCTGCGTGGCGGTCGATACCGGCCTTGGCCACCTGGCGGCGGCACTCGATGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCCTACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTg                                                                                                                                        <  1:285312/149‑1 (MQ=255)
      gcAAAGGCCTGCGTGGCGGTCGATACCGGCCTTGGCCACCTGGCGGCGGCACTCGATGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCCTACGGACGGACCCAGATTCACCAGGCCAGTCACTGGCCtgg                                                                                                                                        <  6:201893/149‑3 (MQ=255)
       cAAAGGCCTGCGTGGCGGTCGATACCGGCCTTGGCCACCTGGCGGCGGCACTCGATGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCCTACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTgc                                                                                                                                       >  1:391953/1‑149 (MQ=255)
       cAAAGGCCTGCGTGGCGGTCGATACCGGCCTTGGCCACCTGGCGGCGGCACTCGATGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCCTACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTgc                                                                                                                                       >  3:57597/1‑149 (MQ=255)
       cAAAGGCCTGCGTGGCGGTCGATACCGGCCTTGGCCACCTGGCGGCGGCACTCGATGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCCTACGGACGGACCCAGATTCACCAGGCCAGTCACTGGCCtggt                                                                                                                                       <  5:193767/149‑4 (MQ=255)
       cAAAGGCCTGCGTGGCGGTCGATACCGGCCTTGGCCACCTGGCGGCGGCACTCGATGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCCTACGGACGGACCCAGATTCACCAGGCCAGTCACTGGCCtggt                                                                                                                                       <  6:361563/149‑4 (MQ=255)
         aaGGCCTGCGTGGCGGTCGATACCGGCCTTGGCCACCTGGCGGCGGCACTCGATGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCCTACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTgcgc                                                                                                                                     >  5:221899/1‑149 (MQ=255)
            gCCTGCGTGGCGGTCGATACCGGCCTTGGCCACCTGGCGGCGGCACTCGATGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCCTACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTcc                                                                                                                                  <  6:164690/149‑1 (MQ=255)
            gCCTGCGTGGCGGTCGATACCGGCCTTGGCCACCTGGCGGCGGCACTCGATGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCCTACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTcc                                                                                                                                  <  4:333381/149‑1 (MQ=255)
              cTGCGTGGCGGTCGATACCGGCCTTGGCCACCTGGCGGCGGCACTCGATGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCCTACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCcct                                                                                                                                >  4:116144/1‑149 (MQ=255)
                 cGTGGCGGTCGATACCGGCCTTGGCCACCTGGCGGCGGCACTCGATGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCCTACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCcctgcc                                                                                                                             >  5:419386/1‑149 (MQ=255)
                 cGTGGCGGTCGATACCGGCCTTGGCCACCTGGCGGCGGCACTCGATGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCCTACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCcctgc                                                                                                                              <  6:212247/148‑1 (MQ=255)
                     gCGGTCGATACCGGCCTTGGCCACCTGGCGGCGGCACTCGATGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCCTACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCctgcctgc                                                                                                                          <  6:434076/148‑1 (MQ=255)
                      cGGTCGATACCGGCCTTGGCCACCTGGCGGCGGCACTCGATGTGCCCACCCTTTTGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCCTACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCcctgcctg                                                                                                                           >  1:176193/1‑146 (MQ=255)
                      cGGTCGATACCGGCCTTGGCCACCTGGCGGCGGCACTCGATGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCCTACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCcctgcctg                                                                                                                           >  4:119520/1‑146 (MQ=255)
                        gTCGATACCGGCCTTGGCCACCTGGCGGCGGCACTCGATGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCCTACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCagaa                                                                                                                      <  4:57597/149‑1 (MQ=255)
                         tCGATACCGGCCTTGGCCACCTGGCGGCGGCACTCGATGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCCTACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCagaag                                                                                                                     <  5:400080/149‑1 (MQ=255)
                           gATACCGGCCTTGGCCACCTGGCGGCGGCACTCGATGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCCTACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCagaag                                                                                                                     >  2:413687/1‑147 (MQ=255)
                            aTACCGGCCTTGGCCACCTGGCGGCGGCACTCGATGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCCTACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGaagaag                                                                                                                  <  2:408295/149‑1 (MQ=255)
                            aTACCGGCCTTGGCCACCTGGCGGCGGCACTCGATGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCCTACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGaagaag                                                                                                                  <  3:281141/149‑1 (MQ=255)
                              aCCGGCCTTGGCCACCTGGCGGCGGCACGTCGATGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCCTACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGt                                                                                                                 >  1:247733/1‑149 (MQ=255)
                              aCCGGCCTTGGCCACCTGGCGGCGGCACTCGATGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCCTACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGt                                                                                                                 <  6:8995/148‑1 (MQ=255)
                              aCCGGCCTTGGCCACCTGGCGGCGGCACTCGATGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCCTACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTg                                                                                                                >  3:197192/1‑149 (MQ=255)
                                cGGCCTTGGCCACCTGGTGGCGGCACTCGATGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCCTACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGt                                                                                                                 >  5:35132/1‑146 (MQ=255)
                                cGGCCTTGGCCACCTGGTGGCGGCACTCGATGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCCTACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGt                                                                                                                 <  6:35132/146‑1 (MQ=255)
                                cGGCCTTGGCCACCTGGCGGCGGCACTCGATGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCCTACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCa                                                                                                              <  2:335034/149‑1 (MQ=255)
                                 ggCCTTGGCCACCTGGCGGCGGCACTCGATGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCCTACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCAc                                                                                                             >  4:26444/1‑149 (MQ=255)
                                       ggCCACCTGGCGGCGGCACTCGATGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCCTACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGCAGTGCACCTACa                                                                                                        >  4:261228/1‑148 (MQ=255)
                                        gCCACCTGGCGGCGGCACTCGATGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCCTACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACaaa                                                                                                      >  7:35209/1‑149 (MQ=255)
                                        gCCACCTGGCGGCGGCACTCGATGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCCTACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCACCCTGCCTGCAGAAGAAGTGCACCTACaaa                                                                                                      <  1:35419/149‑1 (MQ=255)
                                            gcTGGCGGCGGCACTCGATGTGCGCACCCTTTAGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCCTACGGACGGAGCCAGATTCACGAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACTAACCGa                                                                                                  <  2:354535/148‑1 (MQ=255)
                                             cTGGCGGCGGCACTCGATGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCCTACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAg                                                                                                 >  6:263302/1‑149 (MQ=255)
                                                gcggcggcACTCGATGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCCTACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAgc                                                                                                <  1:335916/147‑1 (MQ=255)
                                                    cgtcACTCGATGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCCTACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCTAt                                                                                          <  2:314230/146‑1 (MQ=255)
                                                      gcACTCGATGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCCTACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGa                                                                                        <  3:26444/149‑1 (MQ=255)
                                                        aCTCGATGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCCTACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGAcc                                                                                      >  3:259214/1‑149 (MQ=255)
                                                         cTCGATGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCCTACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCt                                                                                     >  1:2905/1‑149 (MQ=255)
                                                         cTCGATGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCCTACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCt                                                                                     <  2:391953/149‑1 (MQ=255)
                                                              aTGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCCTACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGcc                                                                                 <  6:419386/149‑1 (MQ=255)
                                                                   cccACCCTTTCGCTGTTTGGCCCTCCCAACCCGGGCCTGACTGGCGCCTACGGACGGACCCAGATTCACCAGGCCAATGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGt                                                                             <  3:247987/148‑1 (MQ=255)
                                                                   cccACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCCTACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGtt                                                                            >  1:383198/1‑149 (MQ=255)
                                                                      aCCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCCTACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGa                                                                         >  4:190202/1‑149 (MQ=255)
                                                                       cccTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCCTACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGAt                                                                        <  1:413687/149‑1 (MQ=255)
                                                                          tttCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCACCTACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCgcact                                                                      >  6:274796/1‑144 (MQ=255)
                                                                           ttCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCCTACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTg                                                                     >  3:423980/1‑148 (MQ=255)
                                                                            tCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCCTACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCAATCTGa                                                                    <  1:300129/148‑1 (MQ=255)
                                                                                  tttGGCCCTACCAACCCGGGCCTGACTGGCGCCTACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAAcgc                                                               >  3:68114/1‑147 (MQ=255)
                                                                                  tttGGCCCTACCAACCCGGGCCTGACTGGCGCCTACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAAcgc                                                               >  1:262713/1‑147 (MQ=255)
                                                                                  tttGGCCCTACCAACCCGGGCCTGACTGGCGCCTACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAAcgc                                                               >  1:24958/1‑147 (MQ=255)
                                                                                      gCCCTACCAACCCGGGCCTGACTGGCGCCTACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAACGCGAGTgg                                                         >  6:177998/1‑149 (MQ=255)
                                                                                        ccTACCAACCCGGGCCTGACTGGCGCGTACGGGGGGACGCGGATTCACCAGGCCAGTGAGGGGCGTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGCGCGCCGATGGGCTGGCGCGGTGCGATCTGAAGCGCGAGt                                                           >  8:87529/1‑145 (MQ=255)
                                                                                        ccTACCAACCCGGGCCTGACTGGCGCCTACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAACGCGAGt                                                           <  7:87529/145‑1 (MQ=255)
                                                                                        ccTACCAACCCGGGCCTGACTGGCGCCTACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAACGCGAGTGGcc                                                       >  6:259615/1‑149 (MQ=255)
                                                                                          tACCAACCCGGGCCTGACTGGCGCCTACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAACGCGAGTGGCCGc                                                     >  5:126941/1‑149 (MQ=255)
                                                                                          tACCAACCCGGGCCTGACTGGCGCCTACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAACGCGAGTGGCCGc                                                     >  7:311530/1‑149 (MQ=255)
                                                                                           aCCAACCCGGGCCTGACTGGCGCCTACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAACGCGAGTGGCCGCt                                                    >  5:59142/1‑149 (MQ=255)
                                                                                            ccAACCCGGGCCTGACTGGCGCCTACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAACGCGAGTGGCCGCtg                                                   <  3:211870/149‑1 (MQ=255)
                                                                                             cAACCCGGGCCTGACTGGCGCCTACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGCTCTGAAACGCGAGTGGCCGCtgt                                                  >  7:125211/1‑149 (MQ=255)
                                                                                             cAACCCGGGCCTGACTGGCGCCTACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAACGCGAGTGGCCGCtgt                                                  <  2:247733/149‑1 (MQ=255)
                                                                                             cAACCCGGGCCTGACTGGCGCCTACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAACGCGAGTGGCCGCtgt                                                  <  3:116144/149‑1 (MQ=255)
                                                                                               aCCCGGGCCTGACTGGCGCCTACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAACGCGAGTGGCCGCTGTGc                                                <  2:176193/149‑1 (MQ=255)
                                                                                                       cTGACTGGCGCCTACGGCCGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCGCTGCCTGCAGAAGAAGTGCGCCTACAAACCGAGCGCCGATGCCCTGGCCCGGTTCGCTCTGACACGCGAGTGGCCGCTGTGCTTCACtag                                        >  8:78611/1‑147 (MQ=255)
                                                                                                          aCTGGCGCCTACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCg                                                                                            <  8:42488/94‑1 (MQ=255)
                                                                                                          aCTGGCGCCTACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCg                                                                                            >  7:42488/1‑94 (MQ=255)
                                                                                                          aCTGGCGCCTACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAACGCGAGTGGCCGCTGTGCTTCACTCGc                                       <  4:197192/147‑1 (MQ=255)
                                                                                                          aCTGGCGCCTACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAACGCGAGTGGCCGCTGTGCTTCACTCGc                                       <  2:262713/147‑1 (MQ=255)
                                                                                                          aCTGGCGCCTACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAACGCGAGTGGCCGCTGTGCTTCACTCGCCt                                     <  3:190202/149‑1 (MQ=255)
                                                                                                                gcCTACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAACGCGAGTGGCCGCTGTGCTTCACTCGCCTGAATcc                               >  8:358146/1‑149 (MQ=255)
                                                                                                                gcCTACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAACGCGAGTGGCCGCTGTGCTTCACTCGCCTGAATcc                               >  8:160313/1‑149 (MQ=255)
                                                                                                                gcCTACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTACACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAACGCGAGTGGCCGCTGTGCTTCACTCGCCTGAATcc                               >  4:38534/1‑149 (MQ=255)
                                                                                                                 ccTACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAACGCGAGTGGCCGCTGTGCTTCACTCGCCTGAATccc                              <  2:238726/149‑1 (MQ=255)
                                                                                                                  cTACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTg                                                                                                                <  3:20061/66‑1 (MQ=255)
                                                                                                                  cTACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTg                                                                                                                >  4:20061/1‑66 (MQ=255)
                                                                                                                    acggacggACCCAGATTCGCCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAACGCGAGTGGCCGc                                                     >  7:129672/1‑123 (MQ=255)
                                                                                                                    acggacggACCCAGATTCGCCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAACGCGAGTGGCCGc                                                     <  8:129672/123‑1 (MQ=255)
                                                                                                                     cggacggaCCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAACGCGAGTGGCCGCTGTGCTTCACTCGCCTGAATCCCGAg                           <  4:259214/148‑1 (MQ=255)
                                                                                                                     cggacggaCCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAACGCGAGTGGCCGCTGTGCTTCACTCGCCTGAATCCCGAg                           <  6:221899/148‑1 (MQ=255)
                                                                                                                     cggacggaCCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAACGCGAGTGGCCGCTGTGCTTCACTCGCCTGAATCCCGAGc                          <  3:119520/149‑1 (MQ=255)
                                                                                                                      ggacggacCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAACGCGAGTGGCCGCTGTGCTTCACTCGCCTGAATCCCGAGCa                         <  7:358146/149‑1 (MQ=255)
                                                                                                                        aaggacCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAACGCGAGTGGCCGCTGTGCTTCACTCGCCTGAATCCCGAGCAtg                       >  7:278629/3‑149 (MQ=255)
                                                                                                                              ccAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAACGCGAGTGGCCGCTGTGCTTCACTCGCCTGAATCCCGa                            <  6:260669/138‑1 (MQ=255)
                                                                                                                              ccAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAACGCGAGTGGCCGCTGTGCTTCACTCGCCTGAATCCCGa                            >  5:260669/1‑138 (MQ=255)
                                                                                                                              ccAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAACGCGAGTGGCAGCTGTGCTTCACTCGCCTGAATCCCGAGCATGTGGc                   >  6:19184/1‑147 (MQ=255)
                                                                                                                               cAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAACGCGAGTGGCCGCTGTGCTTCACTCGCCTGAATCCCGAGCATGTGGCGAg                <  6:59142/149‑1 (MQ=255)
                                                                                                                                 gATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAACGCGAGTGGCCGCTGTGCTTCACTCGCCTGAATCCCGAGCATGTGGCGAGcc              >  3:228507/1‑149 (MQ=255)
                                                                                                                                      aCCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAACGCGAGTGGCCGCTGTGCTTCACTCGCCTGAATCCCGAGCATGTGGCGAGCCGCtt          <  2:2905/148‑1 (MQ=255)
                                                                                                                                      aCCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGGTCGATCTGAAACGCGAGTGGCCGCTGTGCTTCACTCGCCTGAATCCCGAGCATGTGGCCAGCCGCtt          >  6:46303/1‑148 (MQ=255)
                                                                                                                                        aaGGCCAGTCACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAACGCGAGTGGCCGCTGTGCTTCACTCGCCTGAATCCCGAGCATGTGGCGAg                >  7:216058/2‑140 (MQ=255)
                                                                                                                                        aaGGCCAGTCACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAACGCGAGTGGCCGCTGTGCTTCACTCGCCTGAATCCCGAGCATGTGGCGAg                <  8:216058/139‑1 (MQ=255)
                                                                                                                                            ccAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAACGCGAGTGGCCGCTGTGCTTCACTCGCCTGAATCCCGAGCATGTGGCGAGCCGCTTGAGcgcg   <  4:68114/149‑1 (MQ=255)
                                                                                                                                             cAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGCGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAACGCGAGTGGCCGCTGTGCTTCACTCGCCTGAATCCCGAGCATGTGGCGAGCCGCTTGAGCGCGt  >  5:11938/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GCGGCGGCAAAGGCCTGCGTGGCGGTCGATACCGGCCTTGGCCACCTGGCGGCGGCACTCGATGTGCCCACCCTTTCGCTGTTTGGCCCTACCAACCCGGGCCTGACTGGCGCCTACGGACGGACCCAGATTCACCAGGCCAGTGACTGGCCTTGCGCTCCCTGCCTGCAGAAGAAGTGCACCTACAAACCGAGCGCCGATGACCTGGCCCGGTTCGATCTGAAACGCGAGTGGCCGCTGTGCTTCACTCGCCTGAATCCCGAGCATGTGGCGAGCCGCTTGAGCGCGT  >  NC_002947/415977‑416265

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: