Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 414,025 T→G 9.0% intergenic (+42/‑161) glnE → / → waaF glutamate‑ammonia‑ligase adenylyltransferase/ADP‑heptose:LPS heptosyltransferase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947414,0250TG 9.0% 787.5 / 50.5 249intergenic (+42/‑161)glnE/waaFglutamate‑ammonia‑ligase adenylyltransferase/ADP‑heptose:LPS heptosyltransferase II
Reads supporting (aligned to +/- strand):  ref base T (70/156);  new base G (6/16);  total (77/172)
Fisher's exact test for biased strand distribution p-value = 8.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.09e-01

Only 100 of 254 total aligned reads displayed.

TCCCACCGCCAGGCCCTGCAGAAGGAGGCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATAAAGAGGGGAGGCCAGGCTGTACCGCAGCCTGTAGCCTCCCCGAGCG  >  NC_002947/413876‑414171
                                                                                                                                                     |                                                                                                                                                   
tCCCACCGCCAGGCCCTGCAGAAGGAGGCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGt                                                                                                                                                     >  1:4546378/1‑150 (MQ=255)
tCCCACCGCCAGGCCCTGCAGAAGGAGGCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACACGt                                                                                                                                                     >  1:833385/1‑150 (MQ=255)
tCCCACCGCCAGGCCCTGCAGAAGGAGGCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCACCACCAGGt                                                                                                                                                     >  1:4940213/1‑150 (MQ=255)
tCCCACCGCCAGGCCCTGCAGAAGGAGGCAGGGGTGAACGATGAGGCACAGTTTGCCGATGAAAGCCGTGAAGTGCGGCGGATCTGGGCTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCACACAGGt                                                                                                                                                     >  2:5104249/1‑150 (MQ=255)
 cccACCGCCAGGCCCTGCAGAAGGAGGCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTa                                                                                                                                                    <  2:6646431/150‑1 (MQ=255)
  ccACCGCCAGGCCCTGCAGAAGGAGGCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTAc                                                                                                                                                   <  1:2464652/150‑1 (MQ=255)
   cACCGCCAGGCCCTGCAGAAGGAGGCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACt                                                                                                                                                  >  1:5894396/1‑150 (MQ=255)
   cACCGCCAGGCCCTGCAGAAGGAGGCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACt                                                                                                                                                  <  1:6040003/150‑1 (MQ=255)
    aCCGCCAGGCCCTGCAGAAGGAGGCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTg                                                                                                                                                 >  2:1672927/1‑150 (MQ=255)
     ccGCCAGGCCCTGCAGAAGGAGGCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGc                                                                                                                                                <  2:7103255/150‑1 (MQ=255)
      cGCCAGGCCCTGCAGAAGGAGGCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGca                                                                                                                                               <  2:7148834/150‑1 (MQ=255)
      cGCCAGGCCCTGCAGAAGGAGGCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGGATTGca                                                                                                                                               <  1:5411281/150‑1 (MQ=255)
      cGCCAGGCCCTGCAGAAGGAGGCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGGATTGca                                                                                                                                               <  1:2916050/150‑1 (MQ=255)
      cGCCAGGCCCTGCAGAAGGAGGCCGAGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGGATTGca                                                                                                                                               <  2:3217835/150‑1 (MQ=255)
        ccAGGCCCTGCAGAAGGAGGCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGcaca                                                                                                                                             <  1:4284006/150‑1 (MQ=255)
        ccAGGCCCTGCAGAAGGAGGCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGcaca                                                                                                                                             <  1:3896125/150‑1 (MQ=255)
        ccAGGCCCTGCAGAAGGAGGCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGcaca                                                                                                                                             <  1:2803910/150‑1 (MQ=255)
          aGGCCCTGCAGAAGGAGGCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGACGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAgg                                                                                                                                          >  2:2173980/1‑150 (MQ=255)
          aGGCCCTGCAGAAGGAGGCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAgg                                                                                                                                          <  1:3037721/150‑1 (MQ=255)
           ggCCCTGCAGAAGGAGGCCGGGGTGATCGATGCGGCACCGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGc                                                                                                                                         <  2:3357343/150‑1 (MQ=255)
            gCCCTGCAGAAGGAGGCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGcc                                                                                                                                        <  2:877890/150‑1 (MQ=255)
            gCCCTGCAGAAGGAGGCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGGATTGcacactgc                                                                                                                                         <  1:5920895/150‑5 (MQ=255)
              ccTGCAGAAGGAGGCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCt                                                                                                                                      <  1:1014975/150‑1 (MQ=255)
              ccTGCAGAAGGAGGCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGACAAGCCCGCTCCCACAGGTa                                                                                                                                                    <  2:6133064/137‑1 (MQ=255)
               cTGCAGAAGGAGGCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTg                                                                                                                                     <  1:3712050/150‑1 (MQ=255)
                tGCAGAAGGAGGCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGa                                                                                                                                    <  1:372515/150‑1 (MQ=255)
                tGCAGAAGGAGGCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGa                                                                                                                                    <  1:1961775/150‑1 (MQ=255)
                atgcttagtccagaaacgctcggggaggctacaggctgcggtacagcctggcctcccctctttatattcaggctatgtggcggacagtgtaccactcacatggccTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGa                                                                                                                                    >  2:802534‑M1/106‑150 (MQ=255)
                 tgcttagtccagaaacgctcggggaggctacaggctgcggtacagcctggcctcccctctttatattcaggctatgtggcggacagtgtaccactcacatggccTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAg                                                                                                                                   >  2:2048418‑M1/105‑150 (MQ=255)
                 gCAGAAGGAGGCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAg                                                                                                                                   <  1:4346707/150‑1 (MQ=255)
                   aGAAGGAGGCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCa                                                                                                                                 <  1:5519981/150‑1 (MQ=255)
                     tagtccagaaacgctcggggaggctacaggctgcggtacagcctggcctcccctctttatattcaggctatgtggcggacagtgtaccactcacatggccTATTCGCGGGCAAGCCCGCTCCCACAGGGATTGCACAGGCCCTGAGCAgt                                                                                                                               >  2:4483890‑M1/101‑150 (MQ=255)
                     aaGGTGGCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGACGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCCAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAgt                                                                                                                               <  1:3600438/150‑1 (MQ=255)
                     aaGGAGGCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAgt                                                                                                                               <  1:7357187/150‑1 (MQ=255)
                       ggaggCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAgtgt                                                                                                                             <  2:6000766/150‑1 (MQ=255)
                       ggaggCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAgtgt                                                                                                                             <  1:6921413/150‑1 (MQ=255)
                        gaggCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTTAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGtgtg                                                                                                                            <  1:1065523/150‑1 (MQ=255)
                        gaggCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGtgtg                                                                                                                            <  2:5889101/150‑1 (MQ=255)
                         aggCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCTCAGGGCCTGTGCAgtacct                                                                                                                           <  1:2279130/150‑5 (MQ=255)
                           gCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCaa                                                                                                                         >  1:1134974/1‑150 (MQ=255)
                            ccGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTAGTCGCGGGCAAGCCCGCTCCCACAGGTACTGCTCAGGGCCTGTGCAgtacctgtg                                                                                                                        <  2:6557565/150‑8 (MQ=255)
                              aacgctcggggaggctacaggctgcggtacagcctggcctcccctctttatattcaggctatgtggcggacagtgtaccactcacatggccTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATcc                                                                                                                      >  1:5339566‑M1/92‑150 (MQ=255)
                              aacgctcggggaggctacaggctgcggtacagcctggcctcccctctttatattcaggctatgtggcggacagtgtaccactcacatggccTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATcc                                                                                                                      >  1:7086351‑M1/92‑150 (MQ=255)
                                ggTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTCAGCTGAGACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCt                                                                                                                    <  1:4973910/150‑1 (MQ=255)
                                cgctcggggaggctacaggctgcggtacagcctggcctcccctctttatattcaggctatgtggcggacagtgtaccactcacatggctTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCt                                                                                                                    >  1:6861134‑M1/90‑150 (MQ=255)
                                cgctcggggaggctacaggctgcggtacagcctggcctcccctctttatattcaggctatgtggcggacagtgtaccactcacatggctTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCt                                                                                                                    >  1:3317587‑M1/90‑150 (MQ=255)
                                cgctcggggaggctacaggctgcggtacagcctggcctcccctctttatattcaggctatgtggcggacagtgtaccactcacatggccTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCt                                                                                                                    >  1:5330042‑M1/90‑150 (MQ=255)
                                 ctGATCGATGCCGCACACTTTGCCGATGAACGCCGTGAACAGCGGCGCATCTGGGGTGAGTTGGGATTAAGCTCATACCGAGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACACGCCCTGAGCAGTGTGCAATCCCtg                                                                                                                   <  2:7090219/149‑1 (MQ=255)
                                  tGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCtgt                                                                                                                  <  1:3597142/150‑1 (MQ=255)
                                  ctcggggaggctacaggctgcggtacagcctggcctcccctctttatattcaggctatgtggcggacagtgtaccactcacatggctTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCtgt                                                                                                                  >  2:1729109‑M1/88‑150 (MQ=255)
                                  ctcggggaggctacaggctgcggtacagcctggcctcccctctttatattcaggctatgtggcggacagtgtaccactcacatggccTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCtgt                                                                                                                  >  2:377651‑M1/88‑150 (MQ=255)
                                    aTCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTgg                                                                                                                <  2:1394748/150‑1 (MQ=255)
                                    aTCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTgg                                                                                                                <  1:4495299/150‑1 (MQ=255)
                                    aTCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTgg                                                                                                                <  1:1604509/150‑1 (MQ=255)
                                    aTCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTgg                                                                                                                <  1:1524930/150‑1 (MQ=255)
                                        gaggctacaggctgcggtacagcctggcctcccctctttatattcaggctatgtggcggacagtgtaccactcacatggccTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGc                                                                                                            >  1:6782429‑M1/82‑150 (MQ=255)
                                          gCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCgg                                                                                                          <  2:291032/150‑1 (MQ=255)
                                           cGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCggg                                                                                                         <  2:3655742/150‑1 (MQ=255)
                                           cGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCggg                                                                                                         <  1:4128589/150‑1 (MQ=255)
                                            ggCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGc                                                                                                        <  2:2708561/150‑1 (MQ=255)
                                             tacaggctgcggtacagcctggcctcccctctttatattcaggctatgtggcggacagtgtaccactcacatggctTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTgg                                                                                                                <  2:1254513‑M1/65‑1 (MQ=255)
                                             tacaggctgcggtacagcctggcctcccctctttatattcaggctatgtggcggacagtgtaccactcacatggccTATTCGCGGGCAAGCCCGCTCCCACAGGGATTGCACAGGCCCTGAGCa                                                                                                                                 >  2:5537545‑M1/77‑124 (MQ=255)
                                             tacaggctgcggtacagcctggcctcccctctttatattcaggctatgtggcggacagtgtaccactcacatggccTATTCGCGGGCAAGCCCGCTCCCACAGGGATTGCACAGGCCCTGAGCAGTGTGCAATcc                                                                                                                      <  1:4893651‑M1/59‑1 (MQ=255)
                                             tacaggctgcggtacagcctggcctcccctctttatattcaggctatgtggcggacagtgtaccactcacatggccTATTCGCGGGCAAGCCCGCTCCCACAGGGATTGCACAGGCCCTGAGCAGTGTGCAATcc                                                                                                                      >  2:4893713‑M1/77‑135 (MQ=255)
                                              cacaGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCtt                                                                                                      <  2:1134974/150‑1 (MQ=255)
                                               acaGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCGGTGGGAGCGGGCTTg                                                                                                     >  1:5739046/1‑150 (MQ=255)
                                                          aTGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCCTATTCGCGGGCAAGCCCGCTCCCACAGGGATTGCACACTGCTCAGGGCCTGTGCAGTACCTGTGGGAGCGGGCTTGCCCGCGAATaa                                                                                          <  1:1225759/150‑2 (MQ=14)
                                                              aCGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAAt                                                                                                                        >  2:6885554/1‑116 (MQ=255)
                                                                   gTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGa                                                                                 <  1:5548156/150‑1 (MQ=255)
                                                                   gTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCGGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGa                                                                                 <  2:3584042/150‑1 (MQ=255)
                                                                                         tatgtggcggacagtgtaccactcacatggccTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCtgt                                                                                                                  >  2:4170310‑M1/33‑95 (MQ=255)
                                                                                          gAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAg                                                          <  2:4546378/150‑1 (MQ=255)
                                                                                               gggTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAAGCCGTGGGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGa                                                     <  2:4359578/150‑1 (MQ=255)
                                                                                                ggTTTAAGCTGATCCCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGCGCAGTGTGCAAGCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGaa                                                    <  1:6549078/150‑1 (MQ=255)
                                                                                                ggTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGaa                                                    <  2:5050321/150‑1 (MQ=255)
                                                                                                ggTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCCCCCACCGGTACTGCACTGGCCCTGAGCAGTGTGCAATCCCTGTGGGCGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGaa                                                    <  1:1410605/150‑1 (MQ=255)
                                                                                                  tttAAGCTGTTACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGGGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAata                                                  <  1:4855142/150‑1 (MQ=255)
                                                                                                  tttAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAata                                                  >  2:5457979/1‑150 (MQ=255)
                                                                                                      cgCTGTTACCGTGTGGGGTTATTCGGGGGCAAGCCCGCTCCCACAGGTAGTGCACTGGCCCTGAGCAGTGTGCAATCCGTGTGGGCGGGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATaaa                                              <  1:3792197/149‑1 (MQ=255)
                                                                                                      aGCTGATACCGTGTGGGCTTATTCGCGGGCAAGCCCGCTCCCACCGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATaaa                                              <  1:1661495/150‑1 (MQ=255)
                                                                                                       gCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTAATGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCAGGTGAGTGGTACACAGTCCGCCACATAGCCGGAATATAAag                                             >  1:6332490/1‑150 (MQ=255)
                                                                                                             aCCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATAAAGAGGGGa                                       <  2:5573595/150‑1 (MQ=255)
                                                                                                                     ggCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATAAAGAGGGGAGGCCAGGc                               >  2:4870338/1‑150 (MQ=255)
                                                                                                                       cTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATAAAGAGGGGAGGCCAGGCTg                             <  1:5550992/150‑1 (MQ=255)
                                                                                                                                     AGCCGGCTCCCTCAGTTACTGCACTGGCGATGTGCTGTGTTCAATCCGTGTGGTAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATAAAGAGGGGAGGCCAGGCTGTACCGCAGCCTGTa               >  1:5104249‑M1/1‑150 (MQ=255)
                                                                                                                                     AGCCCGCTCCCACCGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATAAAGAGGGGAGGCCAGGCTGTACCGCAGCCTGTa               >  2:5753720‑M1/1‑150 (MQ=255)
                                                                                                                                       CCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATAAAGAGGGGAGGCCAGGCTGTACCGCAGCCTGTAGc             >  2:762443‑M1/1‑150 (MQ=255)
                                                                                                                                        CCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATAAAGAGGGGAGGCCAGGCTGTACCGCAGCCTGTAGcc            <  2:47635‑M1/150‑1 (MQ=255)
                                                                                                                                        CCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATAAAGAGGGGAGGCCAGGCTGTACCGCAGCCTGTAGcc            <  1:2649684‑M1/150‑1 (MQ=255)
                                                                                                                                         CGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATAAAGAGGGGAGGCCAGGCTGTACCGCAGCCTGTAGCCt           <  2:5291145‑M1/150‑1 (MQ=255)
                                                                                                                                          GCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAg                                                                                          >  1:5601387‑M1/1‑70 (MQ=255)
                                                                                                                                          GCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATAAAGAGGGGAGGCCAGGCTGTACCGCAGCCTGTAGCCTc          <  1:5304985‑M1/150‑1 (MQ=255)
                                                                                                                                          GCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATAAAGAGGGGAGGCCAGGCTGTACCGCAGCCTGTAGCCTc          >  2:2316890‑M1/1‑150 (MQ=255)
                                                                                                                                           CTCCCACAGGCACTGCCCAGGCCCTGAGCAGTGTGCAAGCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATAAAGAGGGGAGGCCAGGCTGTACCGCAGCCTGTAGCCTcc         >  1:707928‑M1/1‑150 (MQ=255)
                                                                                                                                              CCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATAAAGAGGGGAGGCCAGGCTGTACCGCAGCCTGTAGCCTCCCCg      <  2:6084754‑M1/150‑1 (MQ=255)
                                                                                                                                              CCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATAAAGAGGGGAGGCCAGGCTGTACCGCAGCCTGTAGCCTCCCCg      <  1:1595792‑M1/150‑1 (MQ=255)
                                                                                                                                              CCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATAAAGAGGGGAGGCCAGGCTGTACCGCAGCCTGTAGCCTCCCCg      <  1:1263957‑M1/150‑1 (MQ=255)
                                                                                                                                              CCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCACGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATAAAGAGGGGAGGCCAGGCTGTACCTCAGCCTGTAGCCTCCCCg      <  2:3952983‑M1/150‑1 (MQ=255)
                                                                                                                                                ACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATAAAGAGGGGAGGCCAGGCTGTACCGCAGCCTGTAGCCTCCCCGAg    <  2:1068203‑M1/150‑1 (MQ=255)
                                                                                                                                                  AGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATAAAGAGGGGAGGCCAGGCTGTACCGCAGCCTGTAGCCTCCCCGAGCg  <  1:1722189‑M1/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TCCCACCGCCAGGCCCTGCAGAAGGAGGCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATAAAGAGGGGAGGCCAGGCTGTACCGCAGCCTGTAGCCTCCCCGAGCG  >  NC_002947/413876‑414171

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