Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,002,194:1 +G 15.0% intergenic (‑158/+42) PP_5252 ← / ← PP_5253 hypothetical protein/arylesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,002,1941.G15.0% 738.2 / 154.2 227intergenic (‑158/+42)PP_5252/PP_5253hypothetical protein/arylesterase
Reads supporting (aligned to +/- strand):  ref base . (96/97);  new base G (20/14);  total (116/111)
Fisher's exact test for biased strand distribution p-value = 3.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

Only 100 of 241 total aligned reads displayed.

CCGACTAGGGGCCGGCAGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCACGAATCAG  >  NC_002947/6002045‑6002343
                                                                                                                                                      |                                                                                                                                                     
 cGACTAGGGGCCGGCAGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACaa                                                                                                                                                      >  1:5691449/1‑150 (MQ=255)
 cGACTAGGGGCCGGCAGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACaa                                                                                                                                                      <  2:3714199/150‑1 (MQ=255)
 cGACTAGGGGCCGGCAGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACaa                                                                                                                                                      <  1:1873512/150‑1 (MQ=255)
 cGACTAGGGGCCGGCAGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACaa                                                                                                                                                      <  2:6937993/150‑1 (MQ=255)
  gACTAGGGGCCGGCAGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAc                                                                                                                                                     <  1:5741216/150‑1 (MQ=255)
  gACTAGGGGCCGGCAGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAc                                                                                                                                                     <  1:3857896/150‑1 (MQ=255)
   aCTAGGGGCCGGCAGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACagac                                                                                                                                                     <  2:8579575/150‑4 (MQ=255)
    cTAGGGGCCGGCAGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACagacg                                                                                                                                                    <  2:6242935/150‑5 (MQ=255)
      aGGGGCCGGCAGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACagacgtt                                                                                                                                                  <  2:6174609/150‑7 (MQ=255)
      aGGGGCCGGCAGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCt                                                                                                                                                 <  1:4121719/150‑1 (MQ=255)
         ggCCGGCAGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTa                                                                                                                                              <  1:8447780/150‑1 (MQ=255)
         ggCCGGCAGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTa                                                                                                                                              <  1:1726694/150‑1 (MQ=255)
            cGGCAGGGGGCCTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGa                                                                                                                                           <  2:2421173/150‑1 (MQ=255)
            cGGCAGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGa                                                                                                                                           >  2:1361181/1‑150 (MQ=255)
             ggCAGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTATTCGCGGGTAAACCCGCTCCTACAGACGTTGTAGGAg                                                                                                                                          <  1:3358715/150‑1 (MQ=255)
              gCAGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGACGTTGTAGGAGc                                                                                                                                         <  2:320257/150‑1 (MQ=255)
              gCAGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGACGTTGTAGGAGc                                                                                                                                         <  1:3451037/150‑1 (MQ=255)
              gCAGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGACGTTGTAGGAGc                                                                                                                                         <  2:5118068/150‑1 (MQ=255)
              gCAGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGc                                                                                                                                         <  2:3609149/150‑1 (MQ=255)
                aGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGACGTTGTAGGAGCgg                                                                                                                                       <  1:6244046/150‑1 (MQ=255)
                 gggggACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGTGGGTAAACCCGCTCCTACAGACGTTGTAGGAGCgg                                                                                                                                       <  1:6594200/149‑1 (MQ=255)
                 gggggACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGACGTTGTAGGAGCggg                                                                                                                                      <  2:2718821/150‑1 (MQ=255)
                 gggggACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCggg                                                                                                                                      <  2:3793970/150‑1 (MQ=255)
                 gggggACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCggg                                                                                                                                      <  1:8724852/150‑1 (MQ=255)
                  ggggACTGCGAAGCCTGCCGGTTGATGCACGGCCTGGGCTTGAATGGCTGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGt                                                                                                                                     <  1:7776114/150‑1 (MQ=255)
                   gggACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGtt                                                                                                                                    <  2:5287945/150‑1 (MQ=255)
                        tGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAc                                                                                                                                                     >  2:2115350/1‑128 (MQ=255)
                        tGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAc                                                                                                                                                     <  1:2115350/128‑1 (MQ=255)
                        tGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTAccc                                                                                                                               >  2:8061817/1‑150 (MQ=255)
                        tGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAAACCGCTCCAACAACGTCTGTAGGAGCGGGTTTAacc                                                                                                                               >  1:1188191/1‑150 (MQ=255)
                        tGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAAACCCCTCCAACAACCTCTGTAGGAGCGGGTTTAacc                                                                                                                               >  1:423862/1‑150 (MQ=255)
                        tGCGAAGCATGCCGGTTGACGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGACGTTGTAGGAGCggg                                                                                                                                      >  2:7510970/1‑143 (MQ=255)
                          cGAAGCATGCCGGTTGAGGCAAGGCCTGGGCTTGAATGGCTGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCcgc                                                                                                                             <  1:7978708/150‑1 (MQ=255)
                           gAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCcgcg                                                                                                                            >  2:4250417/1‑150 (MQ=255)
                           gAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCcgcg                                                                                                                            >  2:7265708/1‑150 (MQ=255)
                              gCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAAt                                                                                                                         <  2:1723249/150‑1 (MQ=255)
                                  gCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAAAGttt                                                                                                                     >  2:194719/1‑150 (MQ=255)
                                   ccGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTc                                                                                                                    <  2:1449836/150‑1 (MQ=255)
                                    cGGTTTATGCATGGCCTGGGCTTTATTGGTTGTGCGGTGGGTTTAAGAAGGTTCGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGTTAAACCCGCTCATACAACGTCTGTAGGTGGGGTTTTACCCGCGAATGTTTcc                                                                                                                   <  1:2325374/150‑1 (MQ=255)
                                       ttGATGCAAGGCCTGGGCTTGATTGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCGGGTAAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCcac                                                                                                                <  1:6018573/150‑1 (MQ=255)
                                           tGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCACCTAAAACGACTGTAGGAGCGGGTTTACCCGCGCATGTTTCCCACaaat                                                                                                            >  2:7372285/1‑148 (MQ=255)
                                            gCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCAACCACGTCTGTAGGAGCGGGTTTACCCGCGAATGTGTCCCACAAATa                                                                                                           >  1:5847013/1‑150 (MQ=255)
                                              acggCCTGGGCTTGAATGGCTGGGCGGTGGGTTGAAGAAGATTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGTTAAACCCGCTCCTACAACGTCTGTAGGAGGGGTTTTACCCGCGAATGTTTCCCACAATTAcc                                                                                                         <  1:1559207/148‑1 (MQ=255)
                                                 gCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCCCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCCCt                                                                                                      >  2:1522196/1‑150 (MQ=255)
                                                     gggCTTTATTGTTTGTGCGGTGGGTTGAAGTAGTTTAGCCGTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGTTAAACCCGCTCCTACAACGTCTGTTGGAGGGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAg                                                                                                  <  1:7437014/150‑1 (MQ=255)
                                                           gAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGc                                                                                            >  2:5406145/1‑150 (MQ=255)
                                                              tGGATGTGGGGTGGTTTTAAGAAGGTTCGCCCTGGTCACGCTTTGTGCTGCCTGGGCCGGCCTCTTCGCGGTTAAACCCGCTCCTACAACGTCTGTCGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAg                                                                                         <  1:3745719/150‑1 (MQ=255)
                                                               ggATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGc                                                                                        >  2:7933770/1‑150 (MQ=255)
                                                                    tgCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGaa                                                                                                   >  1:7147052/1‑134 (MQ=255)
                                                                     gCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACCGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTcc                                                                                  <  2:8677598/150‑1 (MQ=255)
                                                                      cGGTGGGTGGAAGACGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCt                                                                                 <  1:4900474/150‑1 (MQ=255)
                                                                       ggTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTc                                                                                <  2:6478306/150‑1 (MQ=255)
                                                                        gTGGGTTGTAGAAGGTTAGCCCTGGTCATGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCa                                                                               <  2:7075783/150‑1 (MQ=255)
                                                                          gggTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCAtt                                                                             <  1:5060481/150‑1 (MQ=255)
                                                                           ggTTGAAGAAGGGTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACCATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTc                                                                            >  2:7781115/1‑150 (MQ=255)
                                                                            gTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTAACCGCGAATGTTTACCACAAATACCGCTGAAGGAACGCCAGCACGTCcacaataa                                                                           >  2:2512482/1‑146 (MQ=38)
                                                                             ttGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGACTGTTTCCCACAAATACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAg                                                                          >  2:2269089/1‑150 (MQ=255)
                                                                             ttGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTACCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAg                                                                          >  2:6711071/1‑150 (MQ=255)
                                                                               gaagaaGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAgct                                                                        <  2:6420971/150‑1 (MQ=255)
                                                                                    agGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACCGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCtgctgc                                                                   <  1:882298/150‑1 (MQ=255)
                                                                                     ggTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAAACCGCACCAACAACGTCTGTAGGAGCGGGTTTACCCGCGCATTTTTCCCACAAATACCGCTGAAGGAACGCCAGCAAGTCCTCATTAAGCTGCTgcg                                                                  >  1:2891741/1‑150 (MQ=38)
                                                                                       ttAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCACCAAGTCCTCATTCAGCTGCTGCGCa                                                                >  2:7769924/1‑150 (MQ=255)
                                                                                        ttgCCCTGGTCAAGATTTTTGGTGCCTGGGCGTTCGTCTTGGTGTTTAAACCCGATCATACAACTTATTTTGTAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCTAGTCCTCATTCAGCTGCTGCGCAt                                                               <  1:4563837/148‑1 (MQ=255)
                                                                                        ccgCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCGGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCAt                                                               <  1:330223/148‑1 (MQ=255)
                                                                                          gCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGa                                                             <  2:4009398/150‑1 (MQ=255)
                                                                                           cccTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACCTCTGGAGGAGCGGGGTTACCCCCCAATGTTTCCCACAATAACCGCTGAACGAACGCCCGCAAGACCTCATTCAGCTGCGGCGCATGAg                                                            >  2:5979272/1‑150 (MQ=38)
                                                                                                      cTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCaaaa                                                 <  2:7047823/150‑1 (MQ=255)
                                                                                                       tttttGCTGCCTGGGCCGGCCTCTTCGCGGGTAAAACCCCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAAATAACGCTGTAGGAACGCCAGCAAGTCCGCATTCAGCTGCGGCGCATGAGTCACAGCAAAAc                                                >  1:691346/1‑150 (MQ=37)
                                                                                                                     gccggccTCTTCGCGGGTAAAACCCCTCCTAAAACGTCGGTAGGAGCGGGTTTACCCCCGAATGTTTCCCACAATTACCCCTGAAGGAACGCCAGCAAGTCCTCCTTCAGCTGCTGCGCATGAGGCACAGCAAAACCGGGCGGCGCCCCg                                  >  2:8260449/1‑150 (MQ=38)
                                                                                                                      ccggccTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCgg                                 <  1:6066645/150‑1 (MQ=255)
                                                                                                                      ccggccTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCgg                                 <  2:6537756/150‑1 (MQ=255)
                                                                                                                          ccTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATa                             <  1:1447694/150‑1 (MQ=255)
                                                                                                                           ctctTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATac                            <  1:5221929/150‑1 (MQ=255)
                                                                                                                              ttCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACAcc                         >  2:3297426/1‑150 (MQ=255)
                                                                                                                               tCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCt                        >  2:6663301/1‑150 (MQ=255)
                                                                                                                               tCGCGGGTAAAACCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCt                        >  2:4096077/1‑150 (MQ=255)
                                                                                                                                  cgGGAAAACCCGCTCAAACATGGTCTTTTGGTTCTGGTTTAACCGCTATTGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCTCTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCa                     <  1:2241445/150‑1 (MQ=255)
                                                                                                                                         aCCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGgc              <  1:122289/150‑1 (MQ=255)
                                                                                                                                         aCCAGCTCCTACAACGTCTGTAGTAGCGGGTTTACACGCGAATGTTTCCCACAATTACAGCTGAAGGAACGCAAGCAAGTCCTCATTCATCTGCTGAGCATGAGTCACTGCAAAACCTAGCGGCGCCCCGGCATACACCTTCAGATCGgc              >  2:8441454/1‑150 (MQ=255)
                                                                                                                                          cccGCTCCTACAGACGTTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGgcg             >  2:4842654/1‑150 (MQ=255)
                                                                                                                                           ccGCTCCTACAGACGTTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCGAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGgcgc            >  2:3119230/1‑150 (MQ=255)
                                                                                                                                           ccGCTCCTACAGACGTTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATacac                          >  1:1042364/1‑136 (MQ=255)
                                                                                                                                           ccGCTCCTACAGACGTTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATacac                          <  2:1042364/136‑1 (MQ=255)
                                                                                                                                           ccGCTCCTACAGACGTTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGgcgc            >  1:2565080/1‑150 (MQ=255)
                                                                                                                                           ccGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGgcgc            >  1:3089401/1‑150 (MQ=255)
                                                                                                                                           ccGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGgcgc            >  2:718053/1‑150 (MQ=255)
                                                                                                                                           ccGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGgcgc            >  1:852345/1‑150 (MQ=255)
                                                                                                                                           ccGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTATCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCCTTCAGCTGCTGCGCATGAGTCCCTGCAAAACCCTGCGGCGCCCCGGCATCCACCTACAGCTCGgcgc            >  2:8181972/1‑150 (MQ=255)
                                                                                                                                             gCTCCTACAGACGTTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACCCCTGCAGCTCGGCGCCa          >  2:1642688/1‑150 (MQ=39)
                                                                                                                                             gCTCCTACAGACGTTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCa          >  1:4777658/1‑150 (MQ=255)
                                                                                                                                             gCTCCTACAGACGTTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTACAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCa          >  1:1464739/1‑150 (MQ=255)
                                                                                                                                             gCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCa          >  2:876434/1‑150 (MQ=255)
                                                                                                                                             gCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCa          >  2:7604646/1‑150 (MQ=255)
                                                                                                                                              cTCCTACAGACGTTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCAc         >  1:2174603/1‑150 (MQ=255)
                                                                                                                                                ccTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCACGa       >  2:4537975/1‑150 (MQ=255)
                                                                                                                                                 cTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAgctg                                                                       <  1:8235671/85‑1 (MQ=255)
                                                                                                                                                  tACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCACGAAt     >  1:5939489/1‑150 (MQ=255)
                                                                                                                                                   aCAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCACGAATc    <  2:2891741/150‑1 (MQ=255)
                                                                                                                                                    cAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCAc                                                        >  1:4442456/1‑97 (MQ=255)
                                                                                                                                                    cAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCACGAATCa   <  1:7703996/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
CCGACTAGGGGCCGGCAGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCACGAATCAG  >  NC_002947/6002045‑6002343

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: