Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 414,039 T→A 5.2% intergenic (+56/‑147) 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,0390TA 5.2% 608.9 / 11.3 200intergenic (+56/‑147)glnE/waaFglutamate‑ammonia‑ligase adenylyltransferase/ADP‑heptose:LPS heptosyltransferase II
Reads supporting (aligned to +/- strand):  ref base T (61/118);  new base A (6/6);  total (74/126)
Fisher's exact test for biased strand distribution p-value = 3.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 205 total aligned reads displayed.

CCTGCAGAAGGAGGCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATAAAGAGGGGAGGCCAGGCTGTACCGCAGCCTGTAGCCTCCCCGAGCGTTTCTGGACTAAGCATG  >  NC_002947/413890‑414188
                                                                                                                                                      |                                                                                                                                                     
 cTGCAGAAGGAGGCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTg                                                                                                                                                      >  2:4950732/1‑150 (MQ=255)
 cTGCAGAAGGAGGCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTGGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCGCCCACAGGTACTGCACAGGCCCTg                                                                                                                                                      >  2:8721265/1‑150 (MQ=255)
  tGCAGAAGGAGGCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGa                                                                                                                                                     <  2:3071503/150‑1 (MQ=255)
  tGCAGAAGGAGGCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGa                                                                                                                                                     <  1:2143180/150‑1 (MQ=255)
    cAGAAGGAGGCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGc                                                                                                                                                   >  1:8033567/1‑150 (MQ=255)
     aGAAGGAGGCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCa                                                                                                                                                  <  1:8435675/150‑1 (MQ=255)
       aaGGAGGCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAgt                                                                                                                                                <  1:1472925/150‑1 (MQ=255)
       aaGGAGGCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAgt                                                                                                                                                >  1:5067590/1‑150 (MQ=255)
         ggaggCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAgtgt                                                                                                                                              <  2:4136691/150‑1 (MQ=255)
           aggCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGc                                                                                                                                            <  1:4792715/150‑1 (MQ=255)
             gCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCaa                                                                                                                                          >  2:4265300/1‑150 (MQ=255)
             gCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCaa                                                                                                                                          <  1:5753578/150‑1 (MQ=255)
             gCCGGGGTGAACGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATAACGTGTTGGCTTAATCGCGGGCAAGCCCGCTCCAACAGGAACTGCACAGGCCCTGAGAAGTGTGCaa                                                                                                                                          >  2:906759/1‑150 (MQ=255)
                ggggTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATcc                                                                                                                                       <  2:1489101/150‑1 (MQ=255)
                 gggTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATccc                                                                                                                                      <  2:4855401/150‑1 (MQ=255)
                 gggTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATccc                                                                                                                                      <  2:1411169/150‑1 (MQ=255)
                 gggTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATccc                                                                                                                                      <  1:1880431/150‑1 (MQ=255)
                  ggTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCt                                                                                                                                     <  2:1395197/150‑1 (MQ=255)
                  ggTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCt                                                                                                                                     <  1:1033606/150‑1 (MQ=255)
                  ggTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCt                                                                                                                                     <  1:3864751/150‑1 (MQ=255)
                    tGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCtgt                                                                                                                                   <  2:6240145/150‑1 (MQ=255)
                      aTCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTgg                                                                                                                                 >  2:5872446/1‑150 (MQ=255)
                        cGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGa                                                                                                                               <  2:1432289/150‑1 (MQ=255)
                          aTGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGa                                                                                                                                                     <  1:7629425/126‑1 (MQ=255)
                          aTGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGc                                                                                                                             <  1:2513327/150‑1 (MQ=255)
                          aTGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATACCTGTGGGAGc                                                                                                                             <  1:1263080/150‑1 (MQ=255)
                            gCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCgg                                                                                                                           <  2:35800/150‑1 (MQ=255)
                            gCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGGATTGCACACTGCTCAGGGCCTGTGCAGTACCTGTGGGAGCgg                                                                                                                           <  2:453622/150‑1 (MQ=18)
                                cacaGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCtt                                                                                                                       <  2:3994060/150‑1 (MQ=255)
                                     tttGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCcg                                                                                                                  <  2:2539004/150‑1 (MQ=255)
                                      ttGCCGATGACCGCCGTGAAGTGCGGCGGTTCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCTGTGTGCAATCCCTGTGGTAGCGGGCTTGCCcgc                                                                                                                 <  2:1701496/150‑1 (MQ=255)
                                        gCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACCGGGCCGGGGCGGGGGGCCATCCCTGTGGGAGCGGGCTTGCCcgcgt                                                                                                               >  1:6881157/1‑149 (MQ=38)
                                        gCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGa                                                                                                               >  2:4416973/1‑150 (MQ=255)
                                           gATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATa                                                                                                            <  2:8033567/150‑1 (MQ=255)
                                            aTGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTTGCGGGCAAGCCCGCTCCCACAGGGATTGCACACTGCTCAGGGCCTGTGCAGTACCTGTGGGAGCGGGCTTGCCCGCGAATaa                                                                                                           <  2:5893106/150‑2 (MQ=14)
                                            aTGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAg                                                                                                           <  2:2683737/150‑1 (MQ=255)
                                                   ccGTGAAGTGCGGCGGATGTGGGGTGCGTGGGGTTTAAGCTGATCCCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACCGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCAtgt                                                                                                    <  1:3563709/150‑1 (MQ=38)
                                                   ccGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCATTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCAtgt                                                                                                    <  2:294245/150‑1 (MQ=255)
                                                      tGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAg                                                                                                 <  2:9145500/150‑1 (MQ=255)
                                                        caGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTg                                                                                               <  2:5085486/149‑1 (MQ=255)
                                                        aaGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTg                                                                                               <  2:8947715/150‑1 (MQ=255)
                                                           tGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTAATGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCTAATATGCCATGTGAGTGGTa                                                                                            <  2:507846/150‑1 (MQ=255)
                                                           tGCGGCGGATCGGGGGGGCGTTGGGTTTAAGCTGATACCGTGTGGGCTTATTCGCGGGCACGCCCGCGCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTa                                                                                            <  1:8274154/150‑1 (MQ=38)
                                                             cggcggATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCCAGTGAGTGGTaca                                                                                          >  2:5015796/1‑150 (MQ=255)
                                                              ggcggATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGAGGTacac                                                                                         >  2:5485545/1‑150 (MQ=255)
                                                               gcggATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACt                                                                                        >  2:4658323/1‑150 (MQ=255)
                                                                   ttCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCACGCCCGCTCCCACCGGTACGGCACCGGCCCTGAGCAGTGTGCAAGCCCTGGGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTcc                                                                                    <  1:6311267/149‑1 (MQ=38)
                                                                   aTCTGGGGTGTGTTGGGTTTAAGCTGTTCCCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACTGGTACTGCACTGGCCCTGAGCTGTGTGCAATCCCTGTGGGCGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTcc                                                                                    <  1:3436496/150‑1 (MQ=38)
                                                                           tGAGTTGGGTTTAAGCTGATCCCGTGTGGGCTTATTCGCGGGCCCGCCCGCCCCCACCGGTACGGCCCCGGCCCGGAGCCGTGTGCAACCCCGGGGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATa                                                                            <  1:587305/150‑1 (MQ=255)
                                                                           tGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGGGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATa                                                                            <  2:3084609/150‑1 (MQ=255)
                                                                               ttGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTa                                                                                            >  2:6012765/1‑130 (MQ=255)
                                                                               ttGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTa                                                                                            <  1:6012765/130‑1 (MQ=255)
                                                                                       aaGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATaa                                                                <  2:4516260/150‑1 (MQ=255)
                                                                                       aaGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCCCGCCCACAGGTACGGCACAGGCCCGGAGCCGTGTGCAATCCCAGTGGGAGCGGGCGTTCCCGCGCATAGGCCATGTGAGTGGTACACAGTCCGCCACATAGCCGGAATATaa                                                                >  1:7852334/1‑150 (MQ=255)
                                                                                        aGCTGATACCGTGTTGGCTTATTCGCGGGCCAGCCCGCTCCCACCGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATaaa                                                               <  2:1825701/150‑1 (MQ=255)
                                                                                        aGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACGGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGCGCGGGCTGGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATaaa                                                               <  1:2755941/150‑1 (MQ=255)
                                                                                         gCTGTTACGGTGTTGGCTTATGGGCGGGCCCGCCCGCCCCCACTGGTAGTGCACCGGCCATGCGCAGGGTGAAACCCCGGTGGGCGGGGGCTTCCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATAAag                                                              <  1:407744/150‑1 (MQ=255)
                                                                                         gCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATAAag                                                              <  2:4575295/150‑1 (MQ=255)
                                                                                           tGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATAAagag                                                            <  2:7669845/150‑1 (MQ=255)
                                                                                           tGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCGCCCACAGGTACTGCACAGGCCCGGGGCCGTGTGCAATCCCTGTGGGTGCGGGCTTGGCCGCGACTAGGGCCTGTGAGTGGGACCCGGTCCGCCACATAGCCGGAATATAAagag                                                            >  1:5599641/1‑150 (MQ=25)
                                                                                            taTACGGTGTTGTCTTTTTCGGGGGCACCCCCGCTCCAACTGGTACTGCACTGGCCCTGCGCTGTGTGCAATCCGTGGGGGCGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATAAAGAgg                                                           <  1:220782/149‑1 (MQ=255)
                                                                                             aTACCGTGTTGGCTTATTCGCGGGCAAGCCCGCGCCCCCAGGTACTGCACAGGCCCGGAGCAGTGTGCAAACCCGGTGGGAGCGGGCTTGCCCGCGAATATGCCAAGTGAGTGGTGAACTGTCCGCCACATAGCCGGGAGATAGAGAggg                                                          >  1:7902981/1‑150 (MQ=37)
                                                                                               cccGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACCGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATAAAGAGGGGa                                                        <  2:6610147/149‑1 (MQ=255)
                                                                                                ccGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATAAAGAGGGGAg                                                       <  2:7877556/150‑1 (MQ=255)
                                                                                                ccGTGTGGGCTTTTTCGGGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGTGATGTGTGCAATCCCTGTGGGAGGGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATAAAGAGGGGAg                                                       <  1:1096350/150‑1 (MQ=38)
                                                                                                   tgtTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATAAAGAGGGGAGGcc                                                    >  2:8853012/1‑150 (MQ=255)
                                                                                                    gtTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATAAAGAGGGGAGGCCa                                                   <  2:1560658/150‑1 (MQ=255)
                                                                                                      tGGCTCATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATAAAGAGGGGAGGCCAgg                                                 <  2:7528147/150‑1 (MQ=255)
                                                                                                          ctATTCGCGGGCAAGCCCGCTCCCACAGGGATTGCACAGGCCCTGAGCAgt                                                                                                                                                <  2:7408862/50‑1 (MQ=14)
                                                                                                          ctATTCGCGGGCAAGCCCGCTCCCACAGGGATTGCACAGGCCCTGAGCAgt                                                                                                                                                >  1:7408862/2‑51 (MQ=14)
                                                                                                            aTTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATAAAGAGGGGAGGCCAGGCTGTAc                                           <  2:6577478/150‑1 (MQ=255)
                                                                                                             ttCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATATAGAGGGGAGGCCAGGCTGTAcc                                          <  2:7711963/150‑1 (MQ=255)
                                                                                                               cgcgGGCAAGCCCGCTCCCACAGGTACAGCACACGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGGGTGCCCGCCAATAGGCCAAGTGAGTGGTACACAGTCCGCCACATAGGCGGTCTATAAAGATGGGAGGCCAGGCTGTACCGc                                        >  1:5237999/1‑150 (MQ=37)
                                                                                                                gcgGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATAAAGAGGGGAGGCCAGGCTGTACCGCa                                       <  2:3724882/150‑1 (MQ=255)
                                                                                                                  gggCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATAAAGAGGGGAGGCCAGGCTGTACCGCAGc                                     >  2:8914756/1‑150 (MQ=255)
                                                                                                                  gggCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATAAAGAGGGGAGGCCAGGCTGTACCGCAGc                                     <  2:7267321/150‑1 (MQ=255)
                                                                                                                    gCAAGCCCGCTCCCACAGGTACTGCACAGGCCCAGAGCAGTGTGCAATCCCAGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACAAAGCCTGAATATAAAGAGGGGAGGCCAGGCTGTACCGCAGCCt                                   >  1:6951858/1‑150 (MQ=255)
                                                                                                                     cAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGTAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATAAAGAGGGGAGGCCAGGCTGTACCGCAGCCTg                                  >  1:8841991/1‑150 (MQ=255)
                                                                                                                          ccGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATAAAGAGGGGAGGCCAGGCTGTACCGCAGCCTGTAGcc                             >  2:8363422/1‑150 (MQ=255)
                                                                                                                              tCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATAAAGAGGGGAGGCCAGGCTGTACCGCAGCCTGTAGCCTccc                         >  2:2430955/1‑150 (MQ=255)
                                                                                                                              tCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATAAAGAGGGGAGGCCAGGCTGTACCGCAGCCTGTAGCCTccc                         >  2:6972003/1‑150 (MQ=255)
                                                                                                                              tCCCACAGGTACTGCACAGGCCCTGAGCAGGGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATAAAGAGGGGAGGCCAGGCTGTACCGCAGCCTGTAGCCTccc                         <  1:5993346/150‑1 (MQ=255)
                                                                                                                                ccACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATAAAGAGGGGAGGCCAGGCTGTACCGCAGCCTGTAGCCTCCCCg                       >  1:6691696/1‑150 (MQ=255)
                                                                                                                                ccACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATAAAGAGGGGAGGCCAGGCTGTACCGCAGCCTGTAGCCTCCCCg                       >  1:9045473/1‑150 (MQ=255)
                                                                                                                                  acaGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATAAAGAGGGGAGGCCAGGCTGTACCGCAGCCTGTAGCCTCCCCGAg                     >  1:4662601/1‑150 (MQ=255)
                                                                                                                                  acaGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATAAAGAGGGGAGGCCAGGCTGTACCGCAGCCTGTAGCCTCCCCGAg                     >  2:4981319/1‑150 (MQ=255)
                                                                                                                                   caGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATAAAGAGGGGAGGCCAGGCTGTACCGCAGCCTGTAGCCTCCCCGAGc                    >  2:1036197/1‑150 (MQ=255)
                                                                                                                                     ggTCCTGCACAGGCCCTGCGCTGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATAAAGAGGGGAGGCCAGGCTGTACCGCAGCCTGTAGCCTCCCCGAGCGt                  <  1:6542329/150‑1 (MQ=255)
                                                                                                                                     ggTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATAAAGAGGGGAGGCCAGGCTGTACCGCAGCCTGTAGCCTCCCCGAGCGt                  >  1:8378924/1‑150 (MQ=255)
                                                                                                                                         cTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATAAAGAGGGGAGGCCAGGCTGTACCGCAGCCTGTAGCCTCCCCGAGCGTTTCt              <  1:4550609/150‑1 (MQ=255)
                                                                                                                                         cTGCACAGGCCCTGAGCAGTGTGCAAGCCCTTTGGGCGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATAAAGAGGGGAGGCCAGGCTGTACCGCAGCCTGTAGCCTCCCCGAGCGTTTCt              <  1:4984513/150‑1 (MQ=255)
                                                                                                                                          tGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGGGTGGTACACTGTCCGCCACATAGCCTGAATATAAAGAGGGGAGGCCAGGCTGTACCGCAGCCTGTAGCCTCCCCGAGCGTTTCTg             >  1:8663587/1‑150 (MQ=255)
                                                                                                                                          tGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATAAAGAGGGGAGGCCAGGCTGTACCGCAGCCTGTAGCCTCCCCGAGCGTTTCTg             >  2:8607278/1‑150 (MQ=255)
                                                                                                                                            cacaGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGGGAGTGGTACACTGTCCGCCACATAGCCTGAATATAAAGAGGGGAGGCCAGGCTGTACCGCAGCCTGTAGCCTCCCCGAGCGTGTCTGGa           >  2:3886318/1‑150 (MQ=255)
                                                                                                                                               aGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATAAAGAGGGGAGGCCAGGCTGTACCGCAGCCTGTAGCCTCCCCGAGCGTTTCTGGACTa        >  1:5983663/1‑150 (MQ=255)
                                                                                                                                                 ggCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATAAAGAGGGGAGGCCAGGCTGTACCGCAGCCTGTAGCCTCCCCGAGCGTTTCTGGACTaa       <  2:5599641/150‑1 (MQ=255)
                                                                                                                                                     cTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATAAAGAGGGGAGGCCAGGCTGTACCGCAGCCTGTAGCCTCCCCGAGCGTTTCTGGACTAAGCAt   <  1:3008730/150‑1 (MQ=255)
                                                                                                                                                     cTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATAAAGAGGGGAGGCCAGGCTGTACCGCAGCCTGTAGCCTCCCCGAGCGTTTCTGGACTAAGCAt   >  1:1715137/1‑150 (MQ=255)
                                                                                                                                                     cTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATAAAGAGGGGAGGCCAGGCTGTACCGCAGCCTGTAGCCTCCCCGAGCGTGTCTGGACTAAGCa  >  1:6527447/1‑149 (MQ=255)
                                                                                                                                                      tGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATAAAGAGGGGAGGCCAGGCTGTACCGCAGCCTGTAGCCTCCCCGAGCGTTTCTGGACTAAGCATg  >  1:2377232/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
CCTGCAGAAGGAGGCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGTGCAATCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCATGTGAGTGGTACACTGTCCGCCACATAGCCTGAATATAAAGAGGGGAGGCCAGGCTGTACCGCAGCCTGTAGCCTCCCCGAGCGTTTCTGGACTAAGCATG  >  NC_002947/413890‑414188

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: