Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,856,513 A→G 100% D203G (GAT→GGT)  relA → ATP:GTP 3'‑pyrophosphotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,856,5130AG99.6% 913.0 / ‑5.2 256D203G (GAT→GGT) relAATP:GTP 3'‑pyrophosphotransferase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (106/149);  total (107/149)
Fisher's exact test for biased strand distribution p-value = 4.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 256 total aligned reads displayed.

TCGACGACGTGCGCGTGGCGCTGATCAAGCTGGCCGAACGCACGTGCGCAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTC  >  NC_002947/1856366‑1856661
                                                                                                                                                   |                                                                                                                                                    
tCGACGACGTGCGCGTGGCGCTGATCAAGCTGGCCGAACGCACGTGCGCAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTc                                                                                                                                                    <  2:4299013/150‑1 (MQ=255)
tCGACGACGTGCGCGTGGCGCTGATCAAGCTGGCCGAACGCACGTGCGCAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTc                                                                                                                                                    >  2:1552487/1‑150 (MQ=255)
 cgacgaCGTGCGCGTGGCGCTGATCAAGCTGGCCGAACGCACGTGCGCAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTca                                                                                                                                                   <  2:6764584/150‑1 (MQ=255)
  gacgacGTGCGCGTGGCGCTGATCAAGCTGGCCGAACGCACGTGCGCAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTcac                                                                                                                                                  >  2:3515848/1‑150 (MQ=255)
  gacgacGTGCGCGTGGCGCTGATCAAGCTGGCCGAACGCACGTGCGCAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTcac                                                                                                                                                  >  1:4439056/1‑150 (MQ=255)
  gacgacGTGCGCGTGGCGCTGATCAAGCTGGCCGAACGCACGTGCGCAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTcac                                                                                                                                                  >  1:1617373/1‑150 (MQ=255)
        gTGCGCGTGGCGCTGATCAAGCTGGCCGAACGCACGTGCGCAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAg                                                                                                                                            >  1:1828761/1‑150 (MQ=255)
          gcgcGTGGCGCTGATCAAGCTGGCCGAACGCACGTGCGCAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTg                                                                                                                                          <  1:3046646/150‑1 (MQ=255)
             cgTGGCGCTGATCAAGCTGGCCGAACGCACGTGCGCAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGa                                                                                                                                       <  2:2308018/150‑1 (MQ=255)
             cgTGGCGCTGATCAAGCTGGCCGAACGCACGTGCGCAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGa                                                                                                                                       <  1:6708737/150‑1 (MQ=255)
              gTGGCGCTGATCAAGCTGGCCGAACGCACGTGCGCAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAg                                                                                                                                      >  1:1690577/1‑150 (MQ=255)
               tGGCGCTGATCAAGCTGGCCGAACGCACGTGCGCAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTGGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGc                                                                                                                                     <  1:5524031/150‑1 (MQ=255)
               tGGCGCTGATCAAGCTGGCCGAACGCACGTGCGCAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGc                                                                                                                                     <  2:3353565/150‑1 (MQ=255)
                ggCGCTGATCAAGCTGGCCGAACGCACGTGCGCAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCt                                                                                                                                    <  2:601142/150‑1 (MQ=255)
                 gcgcTGATCAAGCTGGCCGAACGCACGTGCGCAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTg                                                                                                                                   <  1:4341451/150‑1 (MQ=255)
                 gcgcTGATCAAGCTGGCCGAACGCACGTGCGCAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTg                                                                                                                                   >  2:1025421/1‑150 (MQ=255)
                   gcTGATCAAGCTGGCCGAACGCACGTGCGCAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGa                                                                                                                                 <  1:3987554/150‑1 (MQ=255)
                   gcTGATCAAGCTGGCCGAACGCACGTGCGCAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGa                                                                                                                                 >  2:599608/1‑150 (MQ=255)
                    cTGATCATGCTGGCCGAACGCACGTGCGCAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGaa                                                                                                                                <  1:1103263/150‑1 (MQ=255)
                    cTGATCAAGCTGGCCGAACGCACGTGCGCAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGaa                                                                                                                                <  2:3650076/150‑1 (MQ=255)
                    cTGATCAAGCTGGCCGAACGCACGTGCGCAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGaa                                                                                                                                <  1:6329340/150‑1 (MQ=255)
                       aTCAAGCTGGCCGAACGCACGTGCGCAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGAc                                                                                                                             <  2:4733107/150‑1 (MQ=255)
                       aTCAAGCTGGCCGAACGCACGTGCGCAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGAc                                                                                                                             <  2:6646846/150‑1 (MQ=255)
                        tCAAGCTGGCCGAACGCACGTGCGCAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGAcc                                                                                                                            >  2:6319782/1‑150 (MQ=255)
                         cAAGCTGGCCGAACGCACGTGCGCAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAATACCt                                                                                                                           <  2:6321453/150‑1 (MQ=255)
                         cAAGCTGGCCGAACGCACGTGCGCAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCt                                                                                                                           <  1:1802697/150‑1 (MQ=255)
                         cAAGCTGGCCGAACGCACGTGCGCAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCt                                                                                                                           <  2:5513543/150‑1 (MQ=255)
                          aaGCTGGCCGAACGCACGTGCGCAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTg                                                                                                                          <  2:2448327/150‑1 (MQ=255)
                           aGCTGGCCGAACGCACGTGCGCAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGTCATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGt                                                                                                                         <  1:4917481/150‑1 (MQ=255)
                                gCCGAACGCACGTGCGCAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTc                                                                                                                    <  1:3603061/150‑1 (MQ=255)
                                 ccGAACGCACGTGCGCAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTcc                                                                                                                   >  1:2138461/1‑150 (MQ=255)
                                 ccGAACGCACGTGCGCAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTcc                                                                                                                   <  1:6473812/150‑1 (MQ=255)
                                 ccGAACGCACGTGCGCAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTcc                                                                                                                   <  1:5824814/150‑1 (MQ=255)
                                  cGAACGCACGTGCGCAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCt                                                                                                                                    >  1:4265868/1‑132 (MQ=255)
                                  cGAACGCACGTGCGCAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGaa                                                                                                                                >  1:3469814/1‑136 (MQ=255)
                                  cGAACGCACGTGCGCAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCg                                                                                                                  >  1:3300197/1‑150 (MQ=255)
                                       gCACGTGCGCAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTAcc                                                                                                             <  2:5333633/150‑1 (MQ=255)
                                         aCGTGCGCAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTg                                                                                                           >  1:4626163/1‑150 (MQ=255)
                                         aCGTGCGCAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTg                                                                                                           >  1:1019460/1‑150 (MQ=255)
                                           gTGCGCAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGa                                                                                                         >  1:6628689/1‑150 (MQ=255)
                                             gcgcAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGc                                                                                                       <  2:6499819/150‑1 (MQ=255)
                                             gcgcAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGc                                                                                                       >  2:5951761/1‑150 (MQ=255)
                                              cgcAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGcc                                                                                                      >  1:5877649/1‑150 (MQ=255)
                                              cgcAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGcc                                                                                                      <  2:4650028/150‑1 (MQ=255)
                                               gcAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGccc                                                                                                     >  1:973857/1‑150 (MQ=255)
                                                cAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCg                                                                                                    >  2:1051146/1‑150 (MQ=255)
                                                 aaTCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGa                                                                                                   >  2:6207232/1‑150 (MQ=255)
                                                  aTCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGAc                                                                                                  <  2:2241671/150‑1 (MQ=255)
                                                  aTCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGAc                                                                                                  <  1:3602447/150‑1 (MQ=255)
                                                    ccGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCa                                                                                                <  2:1386145/150‑1 (MQ=255)
                                                     cGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAg                                                                                               >  2:2908663/1‑150 (MQ=255)
                                                     cGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAg                                                                                               >  2:3573733/1‑150 (MQ=255)
                                                          ggTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACaa                                                                                          <  1:5670048/150‑1 (MQ=255)
                                                          ggTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACaa                                                                                          <  1:5348595/150‑1 (MQ=255)
                                                          ggTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCATGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGCACaa                                                                                          <  2:3679751/150‑1 (MQ=255)
                                                            tCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGATCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGc                                                                                        >  1:5637491/1‑150 (MQ=255)
                                                            tCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTACCGCTACATGGAGCCCGACCAGTACAAGc                                                                                        >  1:2840484/1‑150 (MQ=255)
                                                                 tCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATc                                                                                   <  1:2634589/150‑1 (MQ=255)
                                                                   cgccgATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTa                                                                                             <  2:1176920/138‑1 (MQ=255)
                                                                   cgccgATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTa                                                                                             >  1:1176920/1‑138 (MQ=255)
                                                                           aCGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGt                                                                         >  1:6792275/1‑150 (MQ=255)
                                                                             gagaAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTg                                                                       >  2:928223/1‑150 (MQ=255)
                                                                             gagaAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTg                                                                       >  2:1488993/1‑150 (MQ=255)
                                                                             gagaAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTg                                                                       >  1:760826/1‑150 (MQ=255)
                                                                                    gCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAgc                                                                >  2:677657/1‑150 (MQ=255)
                                                                                     ccTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAgcg                                                               <  1:2336452/150‑1 (MQ=255)
                                                                                     ccTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAgcg                                                               >  2:3217727/1‑150 (MQ=255)
                                                                                         cGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGcggcgg                                                           <  2:208008/150‑1 (MQ=255)
                                                                                          gtgtCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCggcggc                                                          <  1:3116307/150‑1 (MQ=255)
                                                                                           tgtCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCt                                                         <  2:4998018/150‑1 (MQ=255)
                                                                                            gtCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTg                                                        <  1:702424/150‑1 (MQ=255)
                                                                                              cgcgcgTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGa                                                      <  1:1794131/150‑1 (MQ=255)
                                                                                               gcgcgTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGAc                                                     <  1:1386034/150‑1 (MQ=255)
                                                                                                  cgTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACcgc                                                  <  1:5956016/150‑1 (MQ=255)
                                                                                                  cgTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACcgc                                                  <  1:5713008/150‑1 (MQ=255)
                                                                                                  cgTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACcgc                                                  <  1:2386037/150‑1 (MQ=255)
                                                                                                   gTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACcgcg                                                 >  1:6549793/1‑150 (MQ=255)
                                                                                                   gTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACcgcg                                                 >  1:1371922/1‑150 (MQ=255)
                                                                                                     gAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAg                                               >  1:5303818/1‑150 (MQ=255)
                                                                                                      aaGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAgc                                              >  2:5649031/1‑150 (MQ=255)
                                                                                                           ttCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTtc                                         >  1:2404901/1‑150 (MQ=255)
                                                                                                             cGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTtcat                                       >  2:4111874/1‑150 (MQ=255)
                                                                                                              gACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTtcatc                                      >  1:2595421/1‑150 (MQ=255)
                                                                                                              gACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTtcatc                                      >  1:4463579/1‑150 (MQ=255)
                                                                                                                    tACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGAc                                >  1:129817/1‑150 (MQ=255)
                                                                                                                    tACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTTCATAAGCGAc                                >  2:5883207/1‑150 (MQ=255)
                                                                                                                      cgcgCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGt                              <  2:4439056/150‑1 (MQ=255)
                                                                                                                           cgcTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGCTCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGtgatga                         <  2:6118699/150‑1 (MQ=255)
                                                                                                                            gcTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGtgatgac                        >  1:3800015/1‑149 (MQ=255)
                                                                                                                            gcTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGaa                        >  2:5231822/1‑150 (MQ=255)
                                                                                                                            gcTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGaa                        >  2:6405613/1‑150 (MQ=255)
                                                                                                                               cgcgCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCa                     <  2:4578047/150‑1 (MQ=255)
                                                                                                                                gcgcACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAg                    >  1:5481946/1‑150 (MQ=255)
                                                                                                                                   cACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTg                 <  1:3287572/150‑1 (MQ=255)
                                                                                                                                     ccGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCa               <  2:39564/150‑1 (MQ=255)
                                                                                                                                      cGCCTGGGTATCGGTCACATCAAGTTGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAATCAGATCGCCAAGCTGTTGCATGAGCTGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAg              <  2:1653881/150‑1 (MQ=255)
                                                                                                                                         cTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAAc           <  2:1828761/150‑1 (MQ=255)
                                                                                                                                           gggTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGa         <  1:4686759/150‑1 (MQ=255)
                                                                                                                                               aTCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTg     <  2:4409737/150‑1 (MQ=255)
                                                                                                                                                  gATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTc  >  2:5095715/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TCGACGACGTGCGCGTGGCGCTGATCAAGCTGGCCGAACGCACGTGCGCAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTC  >  NC_002947/1856366‑1856661

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