Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 414,027 C→T 8.7% intergenic (+44/‑159) 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,0270CT 8.7% 715.3 / 42.5 244intergenic (+44/‑159)glnE/waaFglutamate‑ammonia‑ligase adenylyltransferase/ADP‑heptose:LPS heptosyltransferase II
Reads supporting (aligned to +/- strand):  ref base C (71/150);  new base T (5/16);  total (77/167)
Fisher's exact test for biased strand distribution p-value = 6.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.76e-01

Only 100 of 249 total aligned reads displayed.

CCACCGCCAGGCCCTGCAGAAGGAGGCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATAAAGAGGGGAGGCCAGGCTGTACCGCAGCCTGTAGCCTCCCCGAGCGTTTCT  >  NC_002947/413878‑414176
                                                                                                                                                     |                                                                                                                                                      
ccACCGCCAGGCCCTGCAGAAGGAGGCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTAc                                                                                                                                                        <  1:5098037/150‑1 (MQ=255)
ccACCGCCAGGCCCTGCAGAAGGAGGCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTAc                                                                                                                                                        <  1:260468/150‑1 (MQ=255)
 cACCGCCAGGCCCTGCAGAAGGAGGCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACt                                                                                                                                                       <  1:3025805/150‑1 (MQ=255)
 cACCGCCAGGCCCTGCAGAAGGAGGCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACt                                                                                                                                                       >  1:5894396/1‑150 (MQ=255)
 cACCGCCAGGCCCTGCAGAAGGAGGCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACt                                                                                                                                                       <  2:1074551/150‑1 (MQ=255)
 cACCGCCAGGCCCTGCAGAAGGAGGCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACt                                                                                                                                                       <  1:6040003/150‑1 (MQ=255)
  aCCGCCAGGCCCTGCAGAAGGAGGCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTg                                                                                                                                                      >  2:4495299/1‑150 (MQ=255)
   ccGCCAGGCCCTGCAGAAGGAGGCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGc                                                                                                                                                     <  2:7103255/150‑1 (MQ=255)
    cGCCAGGCCCTGCAGAAGGAGGCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGca                                                                                                                                                    <  1:3801376/150‑1 (MQ=255)
    cGCCAGGCCCTGCAGAAGGAGGCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGGATTGca                                                                                                                                                    <  1:5411281/150‑1 (MQ=255)
    cGCCAGGCCCTGCAGAAGGAGGCCGAGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGGATTGca                                                                                                                                                    <  2:3217732/150‑1 (MQ=255)
       cAGGCCCTGCAGAAGGAGGCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAg                                                                                                                                                <  1:3499267/150‑1 (MQ=255)
        aGGCCCTGCAGAAGGAGGCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGACGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAgg                                                                                                                                               >  2:2173980/1‑150 (MQ=255)
         ggCCCTGCAGAAGGAGGCCGGGGTGATCGATGCGGCACCGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGc                                                                                                                                              <  2:3357343/150‑1 (MQ=255)
         ggCCCTGCAGAAGGAGGCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGc                                                                                                                                              >  2:5544756/1‑150 (MQ=255)
          gCCCTGCAGAAGGAGGCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGcc                                                                                                                                             <  2:4593021/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)
            ccTGCAGAAGGAGGCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCt                                                                                                                                           <  1:2791526/150‑1 (MQ=255)
            ccTGCAGAAGGAGGCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCt                                                                                                                                           <  2:5807004/150‑1 (MQ=255)
              tGCAGAAGGAGGCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCTCAGGGCCTgt                                                                                                                                         <  2:2477192/150‑2 (MQ=255)
              tGCAGAAGGAGGCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGa                                                                                                                                         <  1:372514/150‑1 (MQ=255)
              tGCAGAAGGAGGCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGa                                                                                                                                         <  1:1961775/150‑1 (MQ=255)
              tGCAGAAGGAGGCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGa                                                                                                                                         <  1:4706687/150‑1 (MQ=255)
               gCAGAAGGAGGCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAg                                                                                                                                        <  1:2180161/150‑1 (MQ=255)
                cAGAAGGAGGCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGc                                                                                                                                       <  2:1897989/150‑1 (MQ=255)
                 aGAAGGAGGCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCa                                                                                                                                      <  1:5519981/150‑1 (MQ=255)
                   aaGGAGGCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCTCAGGGCCTGTGCAgt                                                                                                                                    <  2:7224012/150‑1 (MQ=255)
                   aaGGAGGCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAgt                                                                                                                                    <  1:3599717/150‑1 (MQ=255)
                      gaggCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTTAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGtgtg                                                                                                                                 <  1:1065523/150‑1 (MQ=255)
                       aggCCGGGGTGATCGATGCGGCACAGTTTGCCGTTGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGc                                                                                                                                <  1:10964/150‑1 (MQ=255)
                         gCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCaa                                                                                                                              <  1:4323835/150‑1 (MQ=255)
                         gCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCaa                                                                                                                              <  2:4785856/150‑1 (MQ=255)
                            aacgctcggggaggctacaggctgcggtacagcctggcctcccctctttatattcaggctatgtggcggacagtgtaccactcacatggccTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATcc                                                                                                                           >  1:7086326‑M1/92‑150 (MQ=255)
                             gggTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATccc                                                                                                                          <  1:5318088/150‑1 (MQ=255)
                              ggTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCt                                                                                                                         <  2:7096754/150‑1 (MQ=255)
                              cgctcggggaggctacaggctgcggtacagcctggcctcccctctttatattcaggctatgtggcggacagtgtaccactcacatggctTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCt                                                                                                                         >  1:6861134‑M1/90‑150 (MQ=255)
                              cgctcggggaggctacaggctgcggtacagcctggcctcccctctttatattcaggctatgtggcggacagtgtaccactcacatggccTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCt                                                                                                                         >  1:5330042‑M1/90‑150 (MQ=255)
                              cgctcggggaggctacaggctgcggtacagcctggcctcccctctttatattcaggctatgtggcggacagtgtaccactcacatggccTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCt                                                                                                                         >  2:4126427‑M1/90‑150 (MQ=255)
                               gctcggggaggctacaggctgcggtacagcctggcctcccctctttatattcaggctatgtggcggacagtgtaccactcacatggccTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCtg                                                                                                                        >  2:3492025‑M1/89‑150 (MQ=255)
                                tGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCtgt                                                                                                                       <  1:5346367/150‑1 (MQ=255)
                                tGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCtgt                                                                                                                       <  1:817667/150‑1 (MQ=255)
                                ctcggggaggctacaggctgcggtacagcctggcctcccctctttatattcaggctatgtggcggacagtgtaccactcacatggccTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCtgt                                                                                                                       >  1:6860930‑M1/88‑150 (MQ=255)
                                 gATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCtgtg                                                                                                                      <  1:2829813/150‑1 (MQ=255)
                                 gATCGATGCGGCACAGATTGCCGATGAACGCCGTGAACTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCATATTCGCTGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCtgtg                                                                                                                      <  2:6686328/150‑1 (MQ=255)
                                  cggggaggctacaggctgcggtacagcctggcctcccctctttatattcaggctatgtggcggacagtgtaccactcacatggccTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTgg                                                                                                                     >  1:6977781‑M1/86‑150 (MQ=255)
                                  aTCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTgg                                                                                                                     <  1:7344942/150‑1 (MQ=255)
                                    gggaggctacaggctgcggtacagcctggcctcccctctttatattcaggctatgtggcggacagtgtaccactcacatggccTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGa                                                                                                                   >  1:3038760‑M1/84‑150 (MQ=255)
                                    cGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGa                                                                                                                   <  1:6522250/150‑1 (MQ=255)
                                      aTGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGc                                                                                                                 >  1:2316890/1‑150 (MQ=255)
                                        gCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCgg                                                                                                               <  2:291032/150‑1 (MQ=255)
                                        gCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCCTATTCGCGGGCAAGCCCGCTCCCACAGGGATTGCACACTGCTCAGGGCCTGTGCAGTACCTGTGGGAGCgg                                                                                                               <  1:4511886/150‑1 (MQ=17)
                                         cGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCggg                                                                                                              <  2:3655742/150‑1 (MQ=255)
                                           tacaggctgcggtacagcctggcctcccctctttatattcaggctatgtggcggacagtgtaccactcacatggctTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTgg                                                                                                                     >  1:1254513‑M1/77‑141 (MQ=255)
                                           tacaggctgcggtacagcctggcctcccctctttatattcaggctatgtggcggacagtgtaccactcacatggccTATTCGCGGGCAAGCCCGCTCCCACAGGGATTGCACAGGCCCTGAGCAGTGTGCAATcc                                                                                                                           <  1:4893651‑M1/59‑1 (MQ=255)
                                            cacaGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCtt                                                                                                           <  2:1135312/150‑1 (MQ=255)
                                             acaGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCGGTGGGAGCGGGCTTg                                                                                                          >  1:5739046/1‑150 (MQ=255)
                                                gTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTAGTCGCGGGCAAGCCCGCTCCCACAGGTACGGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGccc                                                                                                       <  2:601624/150‑1 (MQ=255)
                                                        aTGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCCTATTCGCGGGCAAGCCCGCTCCCACAGGGATTGCACACTGCTCAGGGCCTGTGCAGTACCTGTGGGAGCGGGCTTGCCCGCGAATaa                                                                                               <  1:1225759/150‑2 (MQ=14)
                                                            aCGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAAt                                                                                                                             <  1:3770632/116‑1 (MQ=255)
                                                            aCGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAAt                                                                                                                             >  2:3770632/1‑116 (MQ=255)
                                                            aCGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAAt                                                                                                                             >  2:6885554/1‑116 (MQ=255)
                                                                cGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCGCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCAagtg                                                                                       >  2:5372078/1‑150 (MQ=255)
                                                                 gTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGa                                                                                      <  1:5548156/150‑1 (MQ=255)
                                                                     aGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTgg                                                                                  <  2:1417211/150‑1 (MQ=255)
                                                                                    gggTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCcac                                                                   >  1:3705647/1‑150 (MQ=255)
                                                                                       tatgtggcggacagtgtaccactcacatggccTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCtgt                                                                                                                       <  1:4170310‑M1/63‑1 (MQ=255)
                                                                                        gAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCTGTTTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAg                                                               <  2:833385/150‑1 (MQ=255)
                                                                                         cgtggggTTTAAGCTGATACCGTGTTGGCTTTTCCGCGGGCACGCCCGCGCCCACCGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCGGTGGGCGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGc                                                              <  1:6617665/146‑1 (MQ=255)
                                                                                         aGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAgtgt                                                                                                                                  <  2:5162601/82‑1 (MQ=255)
                                                                                             gggTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAAGCCGTGGGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGa                                                          <  2:4359578/150‑1 (MQ=255)
                                                                                              ggTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGaa                                                         <  2:5050321/150‑1 (MQ=255)
                                                                                              ggTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCCCCCACCGGTACTGCACTGGCCCTGAGCAGTGTGCAATCCCTGTGGGCGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGaa                                                         <  1:1410605/150‑1 (MQ=255)
                                                                                                    aGCTGATACCGTGTGGGCTTATTCGCGGGCAAGCCCGCTCCCACCGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATaaa                                                   <  1:1661495/150‑1 (MQ=255)
                                                                                                             cGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATAAAGAGGGGAgg                                          >  1:4922820/1‑150 (MQ=255)
                                                                                                                  tGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATAAAGAGGGGAGGCCAgg                                     <  2:1706630/150‑1 (MQ=255)
                                                                                                                  tGGCTTATTCGCGGGCAAGCACGCTCCCACAGGTCCTGCAGAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTTCACTGTCCGCCTCATAGCCTGAATATAAAGAGGGGAGGCCAgg                                     <  2:5478180/150‑1 (MQ=39)
                                                                                                                   ggCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATAAAGAGGGGAGGCCAGGc                                    >  2:4870338/1‑150 (MQ=255)
                                                                                                                    gCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCTTGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATAAAGAGGGGAGGCCAGGCt                                   <  2:4940213/150‑1 (MQ=255)
                                                                                                                     cTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATAAAGAGGGGAGGCCAGGCTg                                  <  1:5550992/150‑1 (MQ=255)
                                                                                                                        ATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACTGGCCGTGAGCTGTGTGCAATCCCTGTGGGCGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATAAAGAGGGGAGGCCAGGCTGTAc                               >  1:4413927‑M1/1‑150 (MQ=255)
                                                                                                                                   AGCCGGCTCCCTCAGTTACTGCACTGGCGATGTGCTGTGTTCAATCCGTGTGGTAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATAAAGAGGGGAGGCCAGGCTGTACCGCAGCCTGTa                    >  1:5104249‑M1/1‑150 (MQ=255)
                                                                                                                                     CCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATAAAGAGGGGAGGCCAGGCTGTACCGCAGCCTGTAGc                  >  2:762443‑M1/1‑150 (MQ=255)
                                                                                                                                      CCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATAAagag                                                <  1:1438165‑M1/119‑1 (MQ=255)
                                                                                                                                      CCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATAAAGAGGGGAGGCCAGGCTGTACCGCAGCCTGTAGcc                 <  2:2100557‑M1/150‑1 (MQ=255)
                                                                                                                                       CGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATAAAGAGGGGAGGCCAGGCTGTACCGCAGCCTGTAGCCt                <  1:1689014‑M1/150‑1 (MQ=255)
                                                                                                                                        GCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAg                                                                                               <  2:5601387‑M1/70‑1 (MQ=255)
                                                                                                                                        GCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAg                                                                                               >  1:5601387‑M1/1‑70 (MQ=255)
                                                                                                                                        GCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATAAAGAGGGGAGGCCAGGCTGTACCGCAGCCTGTAGCCTc               >  2:2316890‑M1/1‑150 (MQ=255)
                                                                                                                                          TCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATAAAGAGGGGAGGCCAGGCTGTACCGCAGCCTGTAGCCTccc             >  1:400785‑M1/1‑150 (MQ=255)
                                                                                                                                            CCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATAAAGAGGGGAGGCCAGGCTGTACCGCAGCCTGTAGCCTCCCCg           <  2:627488‑M1/150‑1 (MQ=255)
                                                                                                                                            CCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATAAAGAGGGGAGGCCAGGCTGTACCGCAGCCTGTAGCCTCCCCg           <  2:6097180‑M1/150‑1 (MQ=255)
                                                                                                                                            CCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATAAAGAGGGGAGGCCAGGCTGTACCGCAGCCTGTAGCCTCCCCg           <  2:4365484‑M1/150‑1 (MQ=255)
                                                                                                                                            CCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATAAAGAGGGGAGGCCAGGCTGTACCGCAGCCTGTAGCCTCCCCg           <  2:136976‑M1/150‑1 (MQ=255)
                                                                                                                                            CCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATAAAGAGGGGAGGCCAGGCTGTACCGCAGCCTGTAGCCTCCCCg           <  1:1595792‑M1/150‑1 (MQ=255)
                                                                                                                                            CCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCACGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATAAAGAGGGGAGGCCAGGCTGTACCTCAGCCTGTAGCCTCCCCg           <  2:3952983‑M1/150‑1 (MQ=255)
                                                                                                                                             CACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGATTGGTACACTGTCCGCCACATAGCCTGAATATAAAGAGGGGAGGCCAGGCTGTACCGCAGCCTGTAGCCTCCCCGa          <  1:1603968‑M1/150‑1 (MQ=255)
                                                                                                                                              ACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATAAAGAGGGGAGGCCAGGCTGTACCGCAGCCTGTAGCCTCCCCGAg         <  2:5902453‑M1/150‑1 (MQ=255)
                                                                                                                                                   caCTGCTCAGGGCCTGTGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATAAAGAGGGGAGGCCAGGCTGTACCGCAGCCTGTAGCCTCCCCGAGCGttt    >  1:6886702/2‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
CCACCGCCAGGCCCTGCAGAAGGAGGCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATAAAGAGGGGAGGCCAGGCTGTACCGCAGCCTGTAGCCTCCCCGAGCGTTTCT  >  NC_002947/413878‑414176

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: