Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,950,510 C→G 5.1% G129A (GGG→GCG)  fliL ← flagellar protein FliL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,950,5100CG 5.1% 285.9 / 4.8 91G129A (GGG→GCG) fliLflagellar protein FliL
Reads supporting (aligned to +/- strand):  ref base C (39/45);  new base G (3/2);  total (43/48)
Fisher's exact test for biased strand distribution p-value = 6.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGGGACAGCAGGTCCTGTACGGCCATTGCGGACTCCTACTGCAATACGAAATTGGTGAACAGCAACTGGTCGATGACCGGCTTGCCGACTTCCTTCTGCGCCACTTCCTGGACCACCGCGGTGGCCTTCTGGCGCAGCATCTCCTGCCCTACCGGGCTGCTGGCGAGGGTGTCGAACCCTTGCCCGGAGAACATCATCACCAGGTTGTTGCGGATCACCGGCATGTGCACCTTGAGGGCATCCAGGTCGGCCTGGTTACGCGCCTGTAGGGTGATACTCAGCTGCATGTAGCGCTGG  >  NC_002947/4950362‑4950659
                                                                                                                                                    |                                                                                                                                                      
cTGGGACAGCAGGTCCTGTACGGCCATTGCGGACTCCTACTGCAATACGAAATTGGTGAACAGCAACTGGTCGATGACCGGCTTGCCGACTTCCTTCTGCGCCACTTCCTGGACCACCGCGGTGGCCTTCTGGCGCAGCATCTCCTGccc                                                                                                                                                       >  1:914297/1‑150 (MQ=255)
   ggACAGCAGGTCCTGTACGGCCATTGCGGACTCCTACTGCAATACGAAATTGGTGAACAGCAACTGGTCGATGACCGGCTTGCCGACTTCCTTCTGCGCCACTTCCTGGACCACCGCGGTGGCCTTCTGGCGCAGCATCTCCTGCCCTAc                                                                                                                                                    <  1:924168/150‑1 (MQ=255)
   ggACAGCAGGTCCTGTACGGCCATTGCGGACTCCTACTGCAATACGAAATTGGTGAACAGCAACTGGTCGATGACCGGCTTGCCGACTTCCTTCTGCGCCACTTCCTGGACCACCGCGGTGGCCTTCTGGCGCAGCATCTCCTGCCCTAc                                                                                                                                                    <  1:1495235/150‑1 (MQ=255)
     aCAGCAGGTCCTGTACGGCCATTGCGGACTCCTACTGCAATACGAAATTGGTGAACAGCAACTGGTCGATGACCGGCTTGCCGACTTCCTTCTGCGCCACTTCCTGGACCACCGCGGTGGCCTTCTGGCGCAGCATCTCCTGCCCTACCg                                                                                                                                                 <  2:563584/150‑1 (MQ=255)
        gcagGTCCTGTACGGCCATTGCGGACTCCTACTGCAATACGAAATTGGTGAACAGCAACTGGTCGATGACCGGCTTGCCGACTTCCTTCTGCGCCACTTCCTGGACCACCGCGGTGGCCTTCTGGCGCAGCATCTCCTGCCCTACCGGgc                                                                                                                                              >  1:1612629/1‑150 (MQ=255)
           ggTCCTGTACGGCCATTGCGGACTCCTACTGCAATACGAAATTGGTGAACAGCAACTGGTCGATGACCGGCTTGCCGACTTCCTTCTGCGCCACTTCCTGGACCACCGCGGTGGCCTTCTGGCGCAGCATCTCCTGCCCTACCGGgctgc                                                                                                                                           >  1:1660690/1‑150 (MQ=255)
             tCCTGTACGGCCATTGCGGACTCCTACTGCAATACGAAATTGGTGAACAGCAACTGATCGATGACCGGCTTGCCGACTTCCTTCTGCGCCACTTCCTGGACCACCGCGGTGGCCTTCTGGCGCAGCATCTCCTGCCCTACCGGGctgctg                                                                                                                                         >  2:1260533/1‑150 (MQ=255)
                 gTACGGCCATTGCGGACTCCTACTGCAATACGAAATTGGTGAACAGCAACTGGTCGATGACCGGCTTGCCGACTTCCTTCTGCGCCACTTCCTGGACCACCGCGGTGGCCTTCTGGCGCAGCATCTCCTGCCCTACCGGGCTGCTGGCGa                                                                                                                                     <  2:489148/150‑1 (MQ=255)
                     ggCCATTGCGGACTCCTACTGCAATACGAAATTGGTGAACAGCAACTGGTCGATGACCGGCTTGCCGACTTCCTTCTGCGCCACTTCCTGGACCACCGCGGTGGCCTTCTGGCGCAGCATCTCCTGCCCTACCGGGCTGCTGGCGAGGgt                                                                                                                                 <  2:1895445/150‑1 (MQ=255)
                     ggCCATTGCGGACTCCTACTGCAATACGAAATTGGTGAACAGCAACTGGTCGATGACCGGCTTGCCGACTTCCTTCTGCGCCACTTCCTGGACCACCGCGGTGGCCTTCTGGCGCAGCATCTCCTGCCCTACCGGGCTGCTGGCGAGGgt                                                                                                                                 >  2:468610/1‑150 (MQ=255)
                        cATTGCGGACTCCTACTGCAATACGAAATTGGTGAACAGCAACTGGTCGATGACCGGCTTGCCGACTTCCTTCTGCGCCACTTCCTGGACCACCGCGGTGGCCTTCTGGCGCAGCATCTCCTGCCCTACCGGgct                                                                                                                                             >  1:685998/1‑135 (MQ=255)
                        cATTGCGGACTCCTACTGCAATACGAAATTGGTGAACAGCAACTGGTCGATGACCGGCTTGCCGACTTCCTTCTGCGCCACTTCCTGGACCACCGCGGTGGCCTTCTGGCGCAGCATCTCCTGCCCTACCGGgct                                                                                                                                             <  2:685998/135‑1 (MQ=255)
                        cATTGCGGACTCCTACTGCAATACGAAATTGGTGAACAGCAACTGGTCGATGACCGGCTTGCCGACTTCCTTCTGCGCCACTTCCTGGACCACCGCGGTGGCCTTCTGGCGCAGCATCTCCTGCCCTACCGGGCTGCTGGCGAGGGTGTc                                                                                                                              <  2:450935/150‑1 (MQ=255)
                        cATTGCGGACTCCTACTGCAATACGAAATTGGTGAACAGCAACTGGTCGATGACCGGCTTGCCGACTTCCTTCTGCGCCACTTCCTGGACCACCGCGGTGGCCTTCTGGCGCAGCATCTCCTGCCCTACCGGGCTGCTGGCGAGGGTGTc                                                                                                                              <  1:583375/150‑1 (MQ=255)
                            gCGGACTCCTACTGCAATACGAAATTGGTGAACAGCAACTGGTCGATGCCTTGCTTGCCGACTTCCTTCTGCGCCACTTCCTGGACCACCGCGGTGGCCTTCTGGCGCAGCATCTCCTGCCCTACCGGGCTGTTGGCGAGGGTGTCGAAc                                                                                                                          <  1:1087635/150‑1 (MQ=255)
                            gCGGACTCCTACTGCAATACGAAATTGGTGAACAGCAACTGGTCGATGACCGGCTTGCCGACTTCCTTCTGCGCCACTTCCTGGACCACCGCGGTGGCCTTCTGGCGCAGCATCTCCTGCCCTACCGGGCTGCTGGCGAGGGTGTCGAAc                                                                                                                          <  1:468610/150‑1 (MQ=255)
                             cGGACTCCTACTGCAATACGAAATTGGTGAACAGCAACTGGTCGATGACCGGCTTGCCGACTTCCTTCTGCGCCACTTCCTGGACCACCGCGGTGGCCTTCTGGCGCAGCATCTCCTGCCCTACCGGGCTGCTGGCGAGGGTGTCGAAcc                                                                                                                         <  2:212672/150‑1 (MQ=255)
                             cGGACTCCTACTGCAATACGAAATTGGTGAACAGCAACTGGTCGATGACCGGCTTGCCGACTTCCTTCTGCGCCACTTCCTGGACCACCGCGGTGGCCTTCTGGCGCAGCATCTCCTGCCCTACCGGGCTGCTGGCGAGGGTGTCGAAcc                                                                                                                         <  1:351358/150‑1 (MQ=255)
                             cGGACTCCTACTGCAATACGAAATGGGAGAACAGCAACTGGTCGATGCCCGGCTTGCCGACTTCCTTCTGCGCCACTTCCTGGACCACCGCGGTGGCCTTCTGGCGCAGCATCTCCTGCCCTAACGGGGTGCTGCCGAGGGTGTGGAAcc                                                                                                                         >  1:651770/1‑150 (MQ=255)
                                  tCCTACTGCAATACGAAATTGGTGACCGGCAACTGGTCGATGACCGGCTTGCCGACTTACTTATGCGCCTCTTCCTAGACGACCGGGGTGGCCTTCTGGCGCAGCATCGCCTTCGCTACCGGGCTGCTGGCGAGCGTGTCGAACCCTTGc                                                                                                                    >  2:255877/1‑150 (MQ=255)
                                  tCCTACTGCAATACGAAATTGGTGAACAGCAACTGGTCGATGACCGGCTTGCCGACTTCCTTCTGCGCCACTTCCTGGACCACCGCGGTGGCCTTCTGGCGCAGCATCTCCTGCCCTACCGGGCTGCTGGCGAGGGTGTCGAACCCtt                                                                                                                      <  1:1715459/148‑1 (MQ=255)
                                  tCCTACTGCAATACGAAATTGGTGAACAGCAACTGGTCGATGACCGGCTTGCCGACTTCCTTCTGCGCCACTTCCTGGACCACCGCGGTGGCCTTCTGGCGCAGCATCTCCTGCCCTACCGGGCTGCTGGCGAGGGTGTCGAACCCtt                                                                                                                      >  2:1715459/1‑148 (MQ=255)
                                    cTACTGCAATACGAAATTGGTGAACAGCAACTGGTCGATGACCGGCTTGCCGACTTCCTTCTGCGCCACTTCCTGGACCACCGCGGTGGCCTTCTGGCGCAGCATCTCCTGCGCTACCGGGCTGCTGGCGAGGGTGTCGAACCCTTGccc                                                                                                                  >  1:1759133/1‑150 (MQ=255)
                                    cTACTGCAATACGAAATTGGTGAACAGCAACTGGTCGATGACCGGCTTGCCGACTTCCTTCTGCGCCACTTCCTGGACCACCGCGGTGGCCTTCTGGCGCAGCATCTCCTGCGCTACCGGGCTGCTGGCGAGGGTGTCGAACCCTTGccc                                                                                                                  >  1:1760542/1‑150 (MQ=255)
                                     tACTGCAATACGAAATTGGTGAACAGCAACTGGTCGATGACCGGCTTGCCGACTTCCTTCTGCGCCACTTCCTGGACCACCGCGGTGGCCTTCTGGCGCAGCATCTCCTGCCCTACCGGGCTGCTGGCGAGGGTGTCGAACCCTTGCCCg                                                                                                                 >  1:1805571/1‑150 (MQ=255)
                                     tACTGCAATACGAAATTGGTGAACAGCAACTGGTCGATGACCGGCTTGCCGACTTCCTTCTGCGCCACTTCCTGGACCACCGCGGTGGCCTTCTGGCGCAGCATCTCCTGCCCTACCGGGCTGCTGGCGAGGGTGTCGAACCCTTGCCCg                                                                                                                 >  2:1619330/1‑150 (MQ=255)
                                      aCTGCAATACGAAATTGGTGAACAGCAACTGGTCGATGACCGGCTTGCCGACTTCCTTCTGCGCCACTTCCTGGACCACCGCGGTGGCCTTCTGGCGCAGCATCTCCTGCCCTACCGGGCTGCTGGCGAGGGTGTCGAACCCTTGCCCgg                                                                                                                >  2:256661/1‑150 (MQ=255)
                                      aCTGCAATACGAAATTGGTGAACAGCAACTGGTCGATGACCGGCTTGCCGACTTCCTTCTGCGCCACTTCCTGGACCACCGCGGTGGCCTTCTGGCGCAGCATCTCCTGCCCTACCGGGCTGCTGGCGAGGGTGTCGAACCCTTGCCCgg                                                                                                                >  1:1712994/1‑150 (MQ=255)
                                       cTGCAATACGAAATTGGTGAACAGCAACTGGTCGATGACCGGCTTGCCGACTTCCTTCTGCGCCACTTCCTGGACCACCGCGGTGGCCTTCTGGCGCAGCATCTCCAGCACTACCGGGCTGCTGGCGAGGGTGTCGCACCCTTGCCCGga                                                                                                               <  1:1214184/150‑1 (MQ=255)
                                             tACGAAATTGGTGAACAGCAACTGGTCGATGACCGGCTTGCCGACTTCCTTCTGCGCCACTTCCTGGACCACCGCGGTGGCCTTCTGGCGCAGCGACGCCGTCCCTACCGGGCTGCTGGCGAGGGTGTCGAACCCTTGCCCGGAGAAcat                                                                                                         <  2:1627350/150‑1 (MQ=255)
                                                gAAATTGGTGAACAGCAACTGGTCGATGACCGGCTTGCCGACTTCCTTCTGCGCCACTTCCTGGACCACCGCGGTGGCCTTCTGGCGCAGCATCTCCTGCCCTACCggg                                                                                                                                               <  1:620000/109‑1 (MQ=255)
                                                gAAATTGGTGAACAGCAACTGGTCGATGACCGGCTTGCCGACTTCCTTCTGCGCCACTTCCTGGACCACCGCGGTGGCCTTCTGGCGCAGCATCTCCTGCCCTACCggg                                                                                                                                               >  2:620000/1‑109 (MQ=255)
                                                 aaaTTGGTGAACAGCAACTGGTCGATGACCGGCTTGCCGACTTCCTTCTGCGCCACTTCCTGGACCACCGCGGTGGCCTTCTGGCGCAGCATCTCCTGCCCTACCGGGCTGCTGGCGAGGGTGTCGAACCCTTGCCCGGAGAACatcatc                                                                                                     <  2:821416/150‑1 (MQ=255)
                                                  aaTTGGTGAACAGCAACTGGTCGATGACCGGCTTGCCGACTTCCTTCTGCGCCACTTCCTGGACCACCGCGGTGGCCTTCTGGCGCAGCATCTCCTGCCCTACCGGGCTGCTGGCGAGGGTGTCGAACCCTTGCCCGGAGAACAtcatca                                                                                                    >  2:3850/1‑150 (MQ=255)
                                                   aTTGGTGAACAGCAACTGGTCGATGACCGGCTTGCCGACTTCCTTCTGCGCCACTTCCTGGACCACCGCGGTGGCCTTCTGGCGCAGCATCTCCTGCCCTACCGGGCTGCTGGCGAGGGTGTCGAACCCTTGCCCGGAGAACATCATCAc                                                                                                   <  1:256661/150‑1 (MQ=255)
                                                       gTGAACAGCAACTGGTCGATGACCGGCTTGCCGACTTCCTTCTGCGCCACTTCCTGGACCACCGCGGTGGCCTTCTGGCGCAGCATCTCCTGCCCTACCGGGCTGCTGGCGAGGGTGTCGAACCCTTGCCCGGAGAACATCATCACCAgg                                                                                               >  2:647821/1‑150 (MQ=255)
                                                             aGCAACTGGTCGATGACCGGCTTGCCGACTTCCTTCTGCGCCACTTCCTGGACCACCGCGGTGGCCTTCTGGCGCAGCATCTCCTGCCCTACCGGGCTGCTGGCGAGGGTGTCGAACCCTTGCCCGGAGAACATCATCACCa                                                                                                 >  2:30633/1‑142 (MQ=255)
                                                             aGCAACTGGTCGATGACCGGCTTGCCGACTTCCTTCTGCGCCACTTCCTGGACCACCGCGGTGGCCTTCTGGCGCAGCATCTCCTGCCCTACCGGGCTGCTGGCGAGGGTGTCGAACCCTTGCCCGGAGAACATCATCACCa                                                                                                 <  1:30633/142‑1 (MQ=255)
                                                              gCAACTGGTCGATGACCGGCTTGCCGACTTCCTTCTGCGCCACTTCCTGGACCACCGCGGTGGCCTTCTGGCGCAGCATCTCCTGCCCTACCGGGCTGCTGGCGAGGGTGTCGAACCCTTGCCCGGAGAACATCATCACCAGGTTGTTGc                                                                                        >  2:265587/1‑150 (MQ=255)
                                                                aaCTGGTCGATGACCGGCTTGCCGACTTCCTTCTGCGCCACTTCCTGGACCACCGCGGTGGCCTTCTGGCGCAGCATCTCCTGCCCTACCGGGCTGCTGGCGAGGGTGTCGAACCCTTGCCCg                                                                                                                 <  1:700034/123‑1 (MQ=255)
                                                                aaCTGGTCGATGACCGGCTTGCCGACTTCCTTCTGCGCCACTTCCTGGACCACCGCGGTGGCCTTCTGGCGCAGCATCTCCTGCCCTACCGGGCTGCTGGCGAGGGTGTCGAACCCTTGCCCg                                                                                                                 <  1:699463/123‑1 (MQ=255)
                                                                aaCTGGTCGATGACCGGCTTGCCGACTTCCTTCTGCGCCACTTCCTGGACCACCGCGGTGGCCTTCTGGCGCAGCATCTCCTGCCCTACCGGGCTGCTGGCGAGGGTGTCGAACCCTTGCCCg                                                                                                                 >  2:700034/1‑123 (MQ=255)
                                                                aaCTGGTCGATGACCGGCTTGCCGACTTCCTTCTGCGCCACTTCCTGGACCACCGCGGTGGCCTTCTGGCGCAGCATCTCCTGCCCTACCGGGCTGCTGGCGAGGGTGTCGAACCCTTGCCCg                                                                                                                 >  2:699463/1‑123 (MQ=255)
                                                                   tGGTCGATGACCGGCTTGCCGACTTCCTTCTGCGCCACTTCCTGGACCACCGCGGTGGCCTTCTGGCGCAGCATCTCCTGCCCTACCGGGCTGCTGGCGAgg                                                                                                                                   >  2:495636/1‑102 (MQ=255)
                                                                   tGGTCGATGACCGGCTTGCCGACTTCCTTCTGCGCCACTTCCTGGACCACCGCGGTGGCCTTCTGGCGCAGCATCTCCTGCCCTACCGGGCTGCTGGCGAgg                                                                                                                                   <  1:495636/102‑1 (MQ=255)
                                                                   tGGTCGATGACCGGCTTGCCGACTTCCTTCTGCGCCACTTCCTGGACCACCGCGGTGGCCTTCTGGCGCAGCATCTCCTGCCCTACCGGGCTGCTGGCGAGGGTGTCGAACCCTTGCCCGGAGAACATCATCACCAGGTTGTTGCGGATc                                                                                   >  1:1050646/1‑150 (MQ=255)
                                                                     gTCGATGACCGGCTTGCCGACTTCCTTCTGCGCCACTTCCTGGACCACCGCGGTGGCCTTCTGGCGCAGCATCTCCTGCCCTACCGGGCTGCTGGCGAGGGTGTCGAACCCTTGCCCGGAGAACATCATCACCAGGTTGTTGCGGATCAc                                                                                 >  2:93548/1‑150 (MQ=255)
                                                                     gTCGATGACCGGCTTGCCGACTTCCTTCTGCGCCACTTCCTGGACCACCGCGGTGGCCTTCTGGCGCAGCATCTCCTGCCCTACCGGGCTGCTGGCGAGGGTGTCGAACCCTTGCCCGGAGAACATCATCACCAGGTTGTTGCGGATCAc                                                                                 >  1:1469441/1‑150 (MQ=255)
                                                                     gTCGATGACCGGCTTGCCGACTTCCTTCTGCGCCACTCAATGGACCACCGCGGTGGCCTTTTGGCGCAGCATCTACTGATGCACTCGTAGCGCTTGGAGGGTGTCGAACCCTTGCCCGGAGAACATCATCACCAGGTTGTTGCGGATCAc                                                                                 >  1:925015/1‑150 (MQ=255)
                                                                        gATGACCGGCTTGCCGACTTCCTTCTGCGCCACTTCCTGGACCACCGCGGTGGCCTTCTGGCGCAGCATCTCCTGCCCTACCGGGCTGCTGGCGa                                                                                                                                     <  2:1868724/95‑1 (MQ=255)
                                                                        gATGACCGGCTTGCCGACTTCCTTCTGCGCCACTTCCTGGACCACCGCGGTGGCCTTCTGGCGCAGCATCTCCTGCCCTACCGGGCTGCTGGCGa                                                                                                                                     >  1:1868724/1‑95 (MQ=255)
                                                                        gATGACCGGCTTGCCGACTTCCTTCTGCGCCACTTCCTGGACCACCGCGGTGGCCTTCTGGCGCAGCATCTCCTGCCCTACCGGGCTGCTGGCGAGGGTGTCGAACCCTTGCCCGGAGAACATCATCACCAGgttgtt                                                                                          >  2:1267135/1‑138 (MQ=255)
                                                                        gATGACCGGCTTGCCGACTTCCTTCTGCGCCACTTCCTGGACCACCGCGGTGGCCTTCTGGCGCAGCATCTCCTGCCCTACCGGGCTGCTGGCGAGGGTGTCGAACCCTTGCCCGGAGAACATCATCACCAGgttgtt                                                                                          <  1:1267135/138‑1 (MQ=255)
                                                                              cGGCTTGCCGACTTCCTTCTGCGCCACTTCCTGGACCACCGCGGTGGCCTTCTGGCGCAGCATCTCCTGCCCTACCGGGCTGCTGGCGAgg                                                                                                                                   <  2:1243302/91‑1 (MQ=255)
                                                                              cGGCTTGCCGACTTCCTTCTGCGCCACTTCCTGGACCACCGCGGTGGCCTTCTGGCGCAGCATCTCCTGCCCTACCGGGCTGCTGGCGAgg                                                                                                                                   >  1:1243302/1‑91 (MQ=255)
                                                                                 cTTGCCGACTTCCTTCTGCGCCACTTCCTGGACCACCGCGGTGGCCTTCTGGCGCAGCATCTCCTGCCCTACCGGGCTGCTGGCGAGGGTGTCGAACCCTTGCCCGGAGAACATCATCACCAGGTTGTTGCGGATCACCGGCATGTGCAc                                                                     <  2:1660690/150‑1 (MQ=255)
                                                                                     ccGACTTCCTTCTGCGCCACTTCCTGGACCACCGCGGTGGCCTTCTGGCGCAGCATCTCCTGCCCTACCGGGCTGCTGGCGAGGGTGTCGAACCCTTGCCCGGAGAACATCATCACCAGGTTGTTGCGGATCACCGGCATGTGCACCTTg                                                                 >  1:609479/1‑150 (MQ=255)
                                                                                      cGACTTCCTTCTGCGCCACTTCCTGGACCACCGCGGTGGCCTTCTGGCGCAGCATCTCCTGCGCTACCGGGCTGCTGGCGAGGGTGTCGAACCCTTGCCCGGAGAACATCATCACCAGGTTGTTGCGGATCACCGGCATGTGCACCTTGa                                                                <  2:1759133/150‑1 (MQ=255)
                                                                                      cGACTTCCTTCTGCGCCACTTCCTGGACCACCGCGGTGGCCTTCTGGCGCAGCATCTCCTGCGCTACCGGGCTGCTGGCGAGGGTGTCGAACCCTTGCCCGGAGAACATCATCACCAGGTTGTTGCGGATCACCGGCATGTGCACCTTGa                                                                <  2:1760542/150‑1 (MQ=255)
                                                                                      cGACTTCCTTCTGCGCCACTTCCTGGACCACCGCGGTGGCCTTCTGGCGCAGCATCTCCTGCCCTACCGGGCTGCTGGCGAGGGTGTCGAACCCTTGCCCGGAGAACATCATCACCAGGTTGTTGCGGATCACCGGCATGTGCACCTTGa                                                                >  2:632516/1‑150 (MQ=255)
                                                                                       gACTTCCTTCTGCGCCACTTCCTGGACCACCGCGGTGGCCTTCTGGCGCAGCATCTCCTGCCCTACCGGGCTGCTGGCGAGGGTGTCGAACCCTTGCCCGGAGAACATCATCACCAGGTTGTTGCGGATCACCGGCATGTGCACCTTGAg                                                               <  2:609479/150‑1 (MQ=255)
                                                                                       gACTTCCTTCTGCGCCACTTCCTGGACCACCGCGGTGGCCTTCTGGCGCAGCATCTCCTGCCCTACCGGGCTGCTGGCGAGGGTGTCGAACCCTTGCCCGGAGAACATCATCACCAGGTTGTTGCGGATACCCCGTACGTGCACCtgtag                                                               <  2:1050646/150‑5 (MQ=255)
                                                                                              ttcTGCGCCACTTCCTGGACCACCGCGGTGGCCTTCTGGCGCAGCATCTCCTGCCCTACCGGGCTGCTGGCGAGGGTGTCGAACCCTTGCCCGGAGAACATCATCACCa                                                                                                 >  1:133630/1‑109 (MQ=255)
                                                                                              ttcTGCGCCACTTCCTGGACCACCGCGGTGGCCTTCTGGCGCAGCATCTCCTGCCCTACCGGGCTGCTGGCGAGGGTGTCGAACCCTTGCCCGGAGAACATCATCACCa                                                                                                 <  2:133630/109‑1 (MQ=255)
                                                                                                 tGCGCCACTTCCTGGACCACCGCGGTGGCCTTCTGGCGCAGCATCTCCTGCCCTACCGGGCTGCTGGCGAGGGTGTCGAACCCTTGCCCGGAGAACATCATCACCAGGTTGTTGCGGATCACCGGCATGTGCACCTTGAGGGCATCCAgg                                                     <  2:1712994/150‑1 (MQ=255)
                                                                                                   cgcCACCTCCCGCACCACCGCGGTGGCCGTCTGGCGCAGCATCTCCTGCCCTACCGGGCTGCCGGCGTGGGTGTCGAACCCTTGCCCGGAGAACCGCATCACCAGGTTGTTGCGGCTCTCCGGCATGTGCACCTTGAGGGCATCCAGGTc                                                   <  1:255877/150‑1 (MQ=255)
                                                                                                    gcCACTTCCTGGACCACCGCGGTGGCCTTCTGGCGCAGCATCTCCTGCCCTACCGGGCTGCTGGCGAGGGTGTCGAACCCTTGCCCGGAGAACATCATCACCAGGTTGTTGCGGATCACCGGCATGTGCACCTTGAGGGCATCCAGGTCg                                                  >  2:1416342/1‑150 (MQ=255)
                                                                                                    gcCACTTCCTGGACCACCGCGGTGGCCTTCTGGCGCAGCATCTCCTGCCCTACCGGGCTGCTGGCGAGGGTGTCGAACCCTTGCCCGGAGAACATCATCACCAGGTTGTTGCGGATCACCGGCATGTGCACCTTGAGGGCATCCAGGTCg                                                  >  2:1893080/1‑150 (MQ=255)
                                                                                                      cACTTCCTGGACCACCGCGGTGGCCTTCTGGCGCAGCATCTCCTGCCCTACCGGGCTGCTGGCGAGGGTGTCGAACCCTTGCCCGGAGAACATCATCACCAGGTTGTTGCGGATCACCGGCATGTGCACCTTGAGGGCATCCAGGTCGGc                                                >  1:1023994/1‑150 (MQ=255)
                                                                                                         ttCCTGGACCACCGCGTTGGCCTTCTGGCGCAGCATCTCCTGCCCTACCGGGCTGCTGGCGAGGGTGTCGAACCCTTGCCCGGAGAACATCATCACCAGGTTGTTGCGGATCACCGGCATGTGCACCTTGAGGGCATCCAGGTCGGCCTg                                             <  2:1805571/150‑1 (MQ=255)
                                                                                                                  caccGCGGTGGCCTTCTGGCGCAGCATCTCCTGCCCTACCGGGCTGCTGGCGAGGGTGTCGAACCCTTGCCCGGAGAACATCATCACCAGGTTGTTGCGGATCACCGGCATGTGCACCTTGAGGGCATCCAGGTCGGcc                                               >  2:1549745/1‑139 (MQ=255)
                                                                                                                  caccGCGGTGGCCTTCTGGCGCAGCATCTCCTGCCCTACCGGGCTGCTGGCGAGGGTGTCGAACCCTTGCCCGGAGAACATCATCACCAGGTTGTTGCGGATCACCGGCATGTGCACCTTGAGGGCATCCAGGTCGGcc                                               <  1:1549745/139‑1 (MQ=255)
                                                                                                                  caccGCGGTGGCCTTCTGGCGCAGCATCTCCTGCCCTACCGGGCTGCTGGCGAGGGTGTCGAACCCTTGCCCGGAGAACATCATCACCAGGTTGTTGCGGATCACCGGCATGTGCACCTTGAGGGCATCCAGGTCGGCCTGGTTACgcgc                                    <  2:1612629/150‑1 (MQ=255)
                                                                                                                  caccGCGGAGGCCTTCTGGCGCAGCATCTCCTGCCCTACCGGGCTGCTGGCGAGGGTGTCGAACCCTTGCCCGGAGAACATCATCACCAGGTTGTTGCGGATAACCGGCATGTGCACCTTGAGGACATCCAGGTAGGCCTGGTTACgcgc                                    >  1:1386489/1‑150 (MQ=255)
                                                                                                                        ggTGGCCTTCTGGCGCAGCATCTCCTGCCCTACCGGGCTGCTGGCGAGGGTGTCGAACCCTTGCCCGGAGAACATCATCACCAGGTTGTTGCGGATCACCGGCATGTGCACCTTGAGGGCATCCAGGTCGGCCTGGTTACGCGCCTGTAg                              <  1:1619330/150‑1 (MQ=255)
                                                                                                                          tGGCCTTCTGGCGCAGCATCTCCTGCCCTACCGGGCTGCTGGCGAGGGTGTCGAACCCTTGCCCGGAGAACATCATCACCAGGTTGTTGCGGATCACCGGCATGTGCACCTTGAGGGCATCCAGGTCGGCCTGGTTACGCGCCTGTAggg                            <  1:647821/150‑1 (MQ=255)
                                                                                                                          tGGCCTTCTGGCGCAGCATCTCCTGCCCTACCGGGCTGCTGGCGAGGGTGTCGAACCCTTGCCCGGAGAACATCATCACCAGGTTGTTGCGGATCACCGGCATGTGCACCTTGAGGGCATCCAGGTCGGCCTGGTTACGCGCCTGTAggg                            <  2:1173622/150‑1 (MQ=255)
                                                                                                                            gCCTTCTGGCGCAGCATCTCCTGCCCTACCGGGCTGCTGGCGAGGGTGTCGAACCCTTGCCCGGAGAACATCATCACCAGGTTGTTGCGGATCACCGGCATGTGCACCTTGAGGGCATCCAGGTCGGCCTGGTTACGCGCCTGTAGGGTg                          <  2:746150/150‑1 (MQ=255)
                                                                                                                             ccTTCTGGCGCAGCATCTCCTGCCCTACCGGGCTGCTGGCGAGGGTGTCGAACCCTTGCCCGGAGAACATCATCACCAGGTTGTTGCGGATCACCGGCATGTGCACCTTGAGGGCATCCAGGTCGGCCTGGTTACGCGCCt                                  >  2:431673/1‑141 (MQ=255)
                                                                                                                             ccTTCTGGCGCAGCATCTCCTGCCCTACCGGGCTGCTGGCGAGGGTGTCGAACCCTTGCCCGGAGAACATCATCACCAGGTTGTTGCGGATCACCGGCATGTGCACCTTGAGGGCATCCAGGTCGGCCTGGTTACGCGCCt                                  <  1:431673/141‑1 (MQ=255)
                                                                                                                               ttCTGGCGCAGCATCTCCTGCCCTACCGGGCTGCTGGCGAGGGTGTCGAACCCTTGCCCGGAGAACATCATCACCAGGTTGTTGCGGATCACCGGCATGTGCACCTTGAGGGCATCCAGGTCGGCCTGGTTACGCGCCTGTAGGGTGATa                       <  2:925015/150‑1 (MQ=255)
                                                                                                                               ttCTGGCGCAGCATCTCCTGCCCTACCGGGCTGCTGGCGAGGGTGTCGAACCCTTGCCCGGAGAACATCATCACCAGGTTGTTGCGGATCACCGGCATGTGCACCTTGAGGGCATCCAGGTCGGCCTGGTTACGCGCCTGTAGGGTGATa                       <  1:632516/150‑1 (MQ=255)
                                                                                                                                   ggCGCAGCATCTCCTGCCCTACCGGGCTGCTGGCGAGGGTGTCGAACCCTTGCCCGGAGAACATCATCACCAGGTTGTTGCGGATCACCGGCATGTGCACCTTGAGGGCATCCAGGTCGGCCTGGTTACGCGCCTGTAGGGTGATACTCa                   <  2:1469441/150‑1 (MQ=255)
                                                                                                                                          caTCTCCTGCCCTACCGGGCTGCTGGCGAGGGTGTCGAACCCTTGCCCGGAGAACATCATCACCAGGTTGTTGCGGATCACCGGCATGTGCACCTTGAGGGCATCCAGGTCGGCCTgg                                            <  2:885730/118‑1 (MQ=255)
                                                                                                                                          caTCTCCTGCCCTACCGGGCTGCTGGCGAGGGTGTCGAACCCTTGCCCGGAGAACATCATCACCAGGTTGTTGCGGATCACCGGCATGTGCACCTTGAGGGCATCCAGGTCGGCCTGGTTACGCGCCTGTAGGGTGATACTCAGCTGCAt            <  1:3850/150‑1 (MQ=255)
                                                                                                                                          caTCTCCTGCCCTACCGGGCTGCTGGCGAGGGTGTCGAACCCTTGCCCGGAGAACATCATCACCAGCTTGTTGCGGATCACCGGCATGTGCACCTTGAGGGCATCCAGGTCGGCCTgg                                            >  1:885730/1‑118 (MQ=255)
                                                                                                                                            tctcCTGCCCTACCGGGCTGCTGGCGAGGGTGTCGAACCCTTGCCCGGAGAACATCATCACCAGGTTGTTGCGGATCACCGGCATGTGCACCTTGAGGGCATCCAGGTCGGCCTGGTTACGCGCCTGTAGGGTGATACTCAGCTGCATGt          >  1:522845/1‑150 (MQ=255)
                                                                                                                                             ctcCTGCCCTACCGGGCTGCTGGCGAGGGTGTCGAACCCTTGCCCGGAGAACATCATCACCAGGTTGTTGCGGATCACCGGCATGTGCACCTTGAGGGCATCCAGGTCGGCCTGGTTACGCGCCTGTAGGGTGATACTCAGCTGCATGTa         >  1:1267466/1‑150 (MQ=255)
                                                                                                                                             ctcCTGCCCTACCGGGCTGCTGGCGAGGGTGTCGAACCCTTGCCCGGAGAACATCATCACCAGGTTGTTGCGGATCACCGGCATGTGCACCTTGAGGGCATCCAGGTCGGCCTGGTTACGCGCCTGTAGGGTGATACTCAGCTGCATGTa         <  2:1267466/150‑1 (MQ=255)
                                                                                                                                                cTGCCCTACCGGGCTGCTGGCGAGGGTGTCGAACCCTTGCCCGGAGAACATCATCACCAGGTTGTTGCGGATCACCGGCATGTGCACCTTGAGGGCATCCAGGTCGGCCTGGTTACGCGCCTGTAGGGTGATACTCAGCTGCATGTAgcg      >  1:87218/1‑150 (MQ=255)
                                                                                                                                                  gCCCTACCGGGCTGCTGGCGAGGGTGTCGAACCCTTGCCCGGAGAACATCATCACCAGGTTGTTGCg                                                                                       <  1:274332/67‑1 (MQ=255)
                                                                                                                                                  gCCCTACCGGGCTGCTGGCGAGGGTGTCGAACCCTTGCCCGGAGAACATCATCACCAGGTTGTTGCg                                                                                       >  2:274332/1‑67 (MQ=255)
                                                                                                                                                    ccTACCGGGCTGCTGGCGAGGGTGTCGAACCCTTGCCCGGAGAACATCATCACCAGGTTGTTGCGGATCACCGGCATGTGCACCTTGAGGGCATCCAGGTCGGCCTGGTTACGCGCCTGTAGGGTGATACTCAGCTGCATGTAGCGCTgg  >  2:1360702/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
CTGGGACAGCAGGTCCTGTACGGCCATTGCGGACTCCTACTGCAATACGAAATTGGTGAACAGCAACTGGTCGATGACCGGCTTGCCGACTTCCTTCTGCGCCACTTCCTGGACCACCGCGGTGGCCTTCTGGCGCAGCATCTCCTGCCCTACCGGGCTGCTGGCGAGGGTGTCGAACCCTTGCCCGGAGAACATCATCACCAGGTTGTTGCGGATCACCGGCATGTGCACCTTGAGGGCATCCAGGTCGGCCTGGTTACGCGCCTGTAGGGTGATACTCAGCTGCATGTAGCGCTGG  >  NC_002947/4950362‑4950659

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