Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 385,479 C→G 7.2% intergenic (+16/+19) PP_0320 → / ← ltaE methyl‑accepting chemotaxis transducer/low specificity l‑threonine aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947385,4790CG 7.2% 577.9 / 15.0 185intergenic (+16/+19)PP_0320/ltaEmethyl‑accepting chemotaxis transducer/low specificity l‑threonine aldolase
Reads supporting (aligned to +/- strand):  ref base C (95/75);  new base G (4/11);  total (99/86)
Fisher's exact test for biased strand distribution p-value = 3.38e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 197 total aligned reads displayed.

CGGACGGTGCAGGGGATTTCCGATTTGGCGCGTGAGACGGCGGCGGAGGTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATGAAGCGCGCGCCACCGCTGCCGATAAAGGTATAGAAGCGCCAGCCCTTGCCCCGCA  >  NC_002947/385332‑385628
                                                                                                                                                   |                                                                                                                                                     
 ggACGGTGCAGGGGATTTCCGATTTGGCGCGTGAGACGGCGGCGGAGGTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCAcggc                                                                                                                                                    <  2:2363786/150‑1 (MQ=255)
 ggACGGTGCAGGGGATTTCCGATTTGGCGCGTGAGACGGCGGCGGAGGTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCAcggc                                                                                                                                                    <  2:4444906/150‑1 (MQ=255)
 ggACGGTGCAGGGGATTTCCGATTTGGCGCGTGAGACGGCGGCGGAGGTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCAcggc                                                                                                                                                    <  1:1845987/150‑1 (MQ=255)
 ggACGGTGCAGGGGATTTCCGATTTGGCGCGTGAGACGGCGGCGGAGGTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCAcggc                                                                                                                                                    <  1:4678618/150‑1 (MQ=255)
  gACGGTGCAGGGGATTTCCGATTTGGCGCGTGAGACGGCGGCGGAGGTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCAcggcg                                                                                                                                                   >  2:1860961/1‑150 (MQ=255)
  gACGGTGCAGGGGATTTCCGATTTGGCGCGTGAGACGGCGGCGGAGGTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCAcggcg                                                                                                                                                   <  1:4591978/150‑1 (MQ=255)
    cGGTGCAGGGGATTTCCGATTTGGCGCGTGAGACGGCGGCGGAGGTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACggcggc                                                                                                                                                 <  1:2746173/150‑1 (MQ=255)
    cGGTGCAGGGGATTTCCGATTTGGCGCGTGAGACGGCGGCGGAGGTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACggcggc                                                                                                                                                 <  2:4275300/150‑1 (MQ=255)
       tGCAGGGGATTTCCGATTTGGCGCGTGAGACGGCGGCGGAGGTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGcccc                                                                                                                                              <  2:5029211/150‑1 (MQ=255)
         cAGGGGATTTCCGATTTGGCGCGTGAGACGGCGGCGGAGGTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCgg                                                                                                                                            <  1:4796131/150‑1 (MQ=255)
          aGGGGATTTCCGATTTGGCGCGTGAGACGGCGGCGGAGGTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGAAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGc                                                                                                                                           <  2:1985930/150‑1 (MQ=255)
              gATTTCCGATTTGGCGCGTGAGACGGCGGCGGAGGTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCtctc                                                                                                                                       >  2:3636908/1‑150 (MQ=255)
                  tCCGATTTGGCGCGTGAGACGGCGGCGGAGGTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTc                                                                                                                                   >  2:51895/1‑150 (MQ=255)
                  tCCGATTTGGCGCGTGAGACGGCGGCGGAGGTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTc                                                                                                                                   >  2:2809403/1‑150 (MQ=255)
                  tCCGATTTGGCGCGTGAGACGGCGGCGGAGGTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTc                                                                                                                                   >  2:2809582/1‑150 (MQ=255)
                   ccGATTGGGCGCGTGAGACGGCGGCGGAGGTGCACGGGTGCCGCGAGGAGGGCCAGGCGTTGCGTGGGTGGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCa                                                                                                                                  <  2:1376879/150‑1 (MQ=255)
                       tttGGCGCGTGAGACGGCGGCGGAGGTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGcca                                                                                                                              <  2:5387055/150‑1 (MQ=255)
                             gcgTGAGACGGCGGCGGAGGTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATg                                                                                                                        >  1:687518/1‑150 (MQ=255)
                             gcgTGAGACGGCGGCGGAGGTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATg                                                                                                                        >  1:3622264/1‑150 (MQ=255)
                              cgTGAGACGGCGGCGGAGGTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGa                                                                                                                       <  1:927806/150‑1 (MQ=255)
                                  agaCGGCGGCGGAGGTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGt                                                                                                                   >  1:3818918/1‑150 (MQ=255)
                                    aCGGCGGCGGAGGTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGc                                                                                                                 <  2:3385973/150‑1 (MQ=255)
                                      ggcggcggAGGTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCgg                                                                                                               >  1:497289/1‑150 (MQ=255)
                                      ggcggcggAGGTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCgg                                                                                                               >  1:2289426/1‑150 (MQ=255)
                                      ggcggcggAGGTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCgg                                                                                                               >  1:1807939/1‑150 (MQ=255)
                                            ggAGGTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATcc                                                                                                         >  1:248521/1‑150 (MQ=255)
                                            ggAGGTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATcc                                                                                                         <  2:1771458/150‑1 (MQ=255)
                                             ctggTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATccg                                                                                                        >  2:824692/3‑150 (MQ=255)
                                               ggTGCAGGGGTGGCGCGCGGTGTGCCTGGGGTTGGGTGGGTTGGCTTTTCACCTGGCGCGGCTGTTGGGTGGTTTCCGGCTATCGCTCGCTGGGGCCGCCCGGCGCCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATccgcc                                                                                                      <  1:109266/150‑1 (MQ=255)
                                               ggTGCAAGGGTGCCGCGTGGAGTGCCTGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCTGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATccgcc                                                                                                      <  2:2614911/150‑1 (MQ=255)
                                                 tGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGTTGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCcgccgc                                                                                                    <  1:2778294/150‑1 (MQ=255)
                                                 tGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCcgccgc                                                                                                    <  2:5584007/150‑1 (MQ=255)
                                                  gCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCgccgcc                                                                                                   >  1:1345678/1‑150 (MQ=255)
                                                  gCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCgccgcc                                                                                                   >  2:3166905/1‑150 (MQ=255)
                                                   cAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCa                                                                                                  >  1:4708991/1‑150 (MQ=255)
                                                    aaGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAg                                                                                                 >  2:3962618/1‑150 (MQ=255)
                                                      gggTGCCGCGAGGTGTGCCTGGGGTTGCGTGGGTTGGCTGATGACCTGGCGGGGCAGTTGGGTGGGTTCAGGCTCTATATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCt                                                                                               <  2:2828763/150‑1 (MQ=255)
                                                       ggTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTc                                                                                              >  2:4861572/1‑150 (MQ=255)
                                                            cgcgAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCa                                                                                         >  1:2310154/1‑150 (MQ=255)
                                                              cgAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCAcc                                                                                       <  2:2772424/150‑1 (MQ=255)
                                                              cgAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTATAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCAcc                                                                                       <  2:2772496/150‑1 (MQ=255)
                                                                 ggagTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCcgc                                                                                    >  1:3913865/1‑150 (MQ=255)
                                                                   tgTGCCTGGGGTTGGGTGGTTTGGCTTTGGACCGGGCGCGGCTGTTGGTGGGTTTCCGGCTCTCGATCGCTGGGGCCGCCGGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCgcgc                                                                                  <  1:3301470/149‑1 (MQ=255)
                                                                   agTGCCAGGCGTTGCGTGGGTTGGCTGTTGACCTGGCGCGGCTGTTGGGTGGGTACAGGCTCTAGATCGCAGGGGCCGCAGGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCgcgc                                                                                  <  1:1287257/150‑1 (MQ=255)
                                                                   agTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCgcgc                                                                                  <  2:3742443/150‑1 (MQ=255)
                                                                         aGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCTGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCttcttc                                                                            <  1:5152159/150‑1 (MQ=255)
                                                                         aGGCGTGGGGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCttcttc                                                                            <  2:544814/150‑1 (MQ=255)
                                                                          ggCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTCGAGCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCg                                                                           <  1:5079060/150‑1 (MQ=255)
                                                                          ggCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCg                                                                           <  1:5626460/150‑1 (MQ=255)
                                                                          ggCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCg                                                                           <  2:4094877/150‑1 (MQ=255)
                                                                          ggCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCg                                                                           <  2:3551733/150‑1 (MQ=255)
                                                                            ggTTGGGTGGGTGGTTTGATGACGTGGCGGGGCTGTTGGGTGGGTTCTGGCTATTGACCGCTGGGGCCGCCCGGCGCCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGt                                                                         <  1:2402061/149‑1 (MQ=255)
                                                                                  gTGGATTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCa                                                                   >  2:2513246/1‑150 (MQ=255)
                                                                                     ggTTGGCTTTTTCCCTGGCGCGGATGTGGGGTGGGTTCCGGCTCTCGATCGCTGGGGCCGCCGGGGGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGa                                                                <  1:545339/150‑1 (MQ=255)
                                                                                      gTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAg                                                               >  1:3171731/1‑150 (MQ=255)
                                                                                          gCTGATGACCTGGCGCGGCAGTTGGGTGGGTTCTGGCTCTAGATCGCTGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGcaca                                                           <  1:4062800/150‑1 (MQ=255)
                                                                                            tgatgaCCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATg                                                         >  2:867081/1‑150 (MQ=255)
                                                                                            tgatgaCCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATg                                                         >  2:2293501/1‑150 (MQ=255)
                                                                                            tgatgaCCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATg                                                         >  2:4000051/1‑150 (MQ=255)
                                                                                             gatgaCCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATGa                                                        >  1:2568046/1‑150 (MQ=255)
                                                                                             gatgaCCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATGa                                                        >  1:3542739/1‑150 (MQ=255)
                                                                                              atgaCCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATGaa                                                       <  2:638245/150‑1 (MQ=255)
                                                                                              atgaCCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATGaa                                                       <  2:4458339/150‑1 (MQ=255)
                                                                                              atgaCCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATGaa                                                       <  2:2754732/150‑1 (MQ=255)
                                                                                                  ccTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATGAAgcgc                                                   >  2:4082488/1‑150 (MQ=255)
                                                                                                  ccTGGCGCGGCAGAGGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGatat                                                                                                           <  1:366359/94‑1 (MQ=255)
                                                                                                    tgggggggCAGATGGGTGGGTTCTGGCTCTTGTTCGCTGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATGAAgcgcgc                                                 <  1:4421460/144‑1 (MQ=255)
                                                                                                    tGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATGAAgcgcgc                                                 <  2:918144/150‑1 (MQ=255)
                                                                                                          ggCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATGAAGCGCGCGCCAcc                                           >  2:2153227/1‑150 (MQ=255)
                                                                                                          ggCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCACCGGCTCTAGCTCAGCCACCAATGATAGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATGAGGCGCGCGCCAcc                                           >  2:4736907/1‑150 (MQ=255)
                                                                                                           gctgttgggtgggTTCAGGCTCTAGAGCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATGAAGCGCGCGCCACCg                                          <  1:851604/145‑1 (MQ=255)
                                                                                                             tgATGGGTGGTTTCAGGCTATTGTTGGCTGGGGCCGCACGGCGGCCCCGTCTCTCTCTCAGCCACCAATGATCGTGCTGATATCCTCCGCCAGCTCACGCACCCGAGCTTCTTCGGTATCCCACGAGCACATGAAGCGCGCGCCACCgat                                        <  2:1635369/149‑3 (MQ=255)
                                                                                                             aGATGGGTGGTTTCAGGCTCTAGATCGCCGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATGAAGCGCGCGCCACCGCt                                        <  1:4363806/150‑1 (MQ=255)
                                                                                                               ttgggtgtgTTCTGGCTCTCGCGCGGGGGGGCCGCCCGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATGAAGCGCGCGCCACCGCTGc                                      <  1:1860961/149‑1 (MQ=255)
                                                                                                               ttgggtgggTTCAGGCTCTAGTTCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATGAAGCGCGCGCCACCGCTGc                                      <  1:2177501/149‑1 (MQ=255)
                                                                                                                ggggtgggtTCCGGCTCTAGAGCGCCGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATGAAGCGCGCGCCACCGCTGcc                                     <  1:1425962/149‑1 (MQ=255)
                                                                                                                  ggtgggtTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATGAAGCGCGCGCCACCGCTGCCGa                                   >  1:3411349/1‑150 (MQ=255)
                                                                                                                       gtTCCGGCTCTAGATCGCCGGGGCCGCAGGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATGAAGCGCGCGCCACCGCTGCCGATAAAg                              <  1:1084403/150‑1 (MQ=255)
                                                                                                                       gtTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATGAAGCGCGCGCCACCGCTGCCGATaaat                              <  1:5029069/150‑2 (MQ=255)
                                                                                                                       gtTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATGAAGCGCGCGCCACCGCTGCCGATAAAg                              >  2:1052787/1‑150 (MQ=255)
                                                                                                                       gtTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATGAAGCGCGCGCCACCGCTGCCGATAAAg                              <  2:4545347/150‑1 (MQ=255)
                                                                                                                        ttCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATGAAGCGCGCGCCACCGCTGCCGATAAAgg                             >  1:3090593/1‑150 (MQ=255)
                                                                                                                           aGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATGAAGCGCGCGCCACCGCTGCCGATAAAGGtct                          >  2:5408939/1‑148 (MQ=255)
                                                                                                                           aGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATGAAGCGCGCGCCACCGCTGCCGATAAAGGtat                          >  2:553321/1‑150 (MQ=255)
                                                                                                                           aGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATGAAGCGCGCGCCACCGCTGCCGATAAAGGtat                          >  1:2633325/1‑150 (MQ=255)
                                                                                                                             gCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATGAAGCGCGCGCCACCGCTGCCGATAAAGGTATAg                        >  2:3797417/1‑150 (MQ=255)
                                                                                                                                 tAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATGAAGCGCGCGCCACCGCTGCCGATAAAGGTATAGAAgc                    >  2:5051714/1‑150 (MQ=255)
                                                                                                                                 tAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATGAAGCGCGCGCCACCGCTGCCGATAAAGGTATAGAAgc                    >  2:2399681/1‑150 (MQ=255)
                                                                                                                                  aGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATGAAGCGCGCGCCACCGCTGCCGATAAAGGTATAGAAgcg                   >  1:1728164/1‑150 (MQ=255)
                                                                                                                                     tCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATGAAGCGCGCGCCACCGCTGCCGATAAAGGTATAGAAGCGCCa                <  2:300054/150‑1 (MQ=255)
                                                                                                                                     tCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATGAAGCGCGCGCCACCGCTGCCGATAAAGGTATAGAAGCGCCa                <  1:4156582/150‑1 (MQ=255)
                                                                                                                                     tCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATGAAGCGCGCGCCACCGCTGCCGATAAAGGTATAGAAGCGCCa                >  1:4922567/1‑150 (MQ=255)
                                                                                                                                        ccggggCCGCCGTGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATGAAGCGCGCGCCACCGCTGCCGATAAAGGTATAGAAGCGCCAGcc             >  1:2032081/3‑150 (MQ=255)
                                                                                                                                           gggCCGCCGTGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATGAAGCGCGCGCCACCGCTGCCGATAAAGGTATAGAAGCGCCAGCCCtt          >  1:2405783/1‑150 (MQ=255)
                                                                                                                                             gCCGCAGGGCGGCCCCGGCTCTCGCTCAGCCACCATTGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATGAAGCGCGCGCCACCGCTGCCGATAAAGGTATAGAAGCGCCAGCCCTTGc        <  1:5411351/150‑1 (MQ=255)
                                                                                                                                             gCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATGAAGCGCGCGCCACCGCTGCCGATAAAGGTATAGAAGCGCCAGCCCTTGc        >  2:4695803/1‑150 (MQ=255)
                                                                                                                                              ccGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCAc                                                                  >  1:356530/1‑91 (MQ=255)
                                                                                                                                               cGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATGAAGCGCGCGCCACCGCTGCCGATAAAGGTATAGAAGCGCCAGCCCTTGccc      <  2:1498711/150‑1 (MQ=255)
                                                                                                                                                gCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATGAAGCGCGCGCCACCGCTGCCGATAAAGGTATAGAAGCGCCAGCCCTTGcccc     >  1:2849718/1‑150 (MQ=255)
                                                                                                                                                gCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATGAAGCGCGCGCCACCGCTGCCGATAAAGGTATAGAAGCGCCAGCCCTTGcccc     >  2:4931152/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CGGACGGTGCAGGGGATTTCCGATTTGGCGCGTGAGACGGCGGCGGAGGTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATGAAGCGCGCGCCACCGCTGCCGATAAAGGTATAGAAGCGCCAGCCCTTGCCCCGCA  >  NC_002947/385332‑385628

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: