Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,342,049	0	A	G	36.7%	90.6 / 77.4	79	intergenic (+68/‑263)	ribAB‑II/PP_3814	3,4‑dihydroxy‑2‑butanone‑4‑phosphate synthase/polyamine ABC transportersubstrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (25/25);  new base G (13/16);  total (38/41)
Fisher's exact test for biased strand distribution p-value = 8.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CCGCTGAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGC‑GGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCG‑CCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGTCTGCAGGCCGGC  >  NC_002947/4341901‑4342198                                                                                                                                                      |                                                                                                                                                       ccgcTGAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGC‑GGGTTTACCCGCGAATGCGAGGGTGacac                                                                                                                                                       <  2:2643395/150‑1 (MQ=255) ccgcTGAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGC‑GGGTTTACCCGCGAATGCGAGGGTGacac                                                                                                                                                       <  2:2643438/150‑1 (MQ=255)   gcTGAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGC‑GGGTTTACCCGCGAATGCGAGGGTGAcacca                                                                                                                                                     <  2:56576/150‑1 (MQ=255)    cTGAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGC‑GGGTTTACCCGCGAATGCGAGGGTGAcaccac                                                                                                                                                    <  2:1940871/150‑1 (MQ=255)    cTGAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGC‑GGGTTTACCCGCGAATGCGAGGGTGAcaccac                                                                                                                                                    <  1:2850497/150‑1 (MQ=255)       aaaTACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGC‑GGGTTTACCCGCGAATGCGAGGGTGACACCACCGa                                                                                                                                                 <  1:2172947/150‑1 (MQ=255)        aaTACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGC‑GGGTTTACCCGCGAATGCGAGGGTGACACCACCGAc                                                                                                                                                <  2:647668/150‑1 (MQ=255)        aaTACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGC‑GGGTTTACCCGCGAATGCGAGGGTGACACCACCGAc                                                                                                                                                <  2:1032520/150‑1 (MQ=255)        aaTACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGC‑GGGTTTACCCGCGAATGCGAGGGTGACACCACCGAc                                                                                                                                                <  2:1721113/150‑1 (MQ=255)        aaTACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGC‑GGGTTTACCCGCGAATGCGAGGGTGACACCACCGAc                                                                                                                                                <  1:3604845/150‑1 (MQ=255)        aaTACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGC‑GGGTTTACCCGCGAATGCGAGGGTGACACCACCGAc                                                                                                                                                <  2:3766095/150‑1 (MQ=255)                gggCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGC‑GGGTTTACCCGCGAATGCGTCGGTGGTGTCACCCTCGCATTcgc                                                                                                                                        <  1:782749/150‑1 (MQ=25)                gggCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGC‑GGGTTTACCCGCGAATGCGTCGGTGGTGTCACCCTCGCATTcgc                                                                                                                                        <  1:782716/150‑1 (MQ=25)                gggCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGC‑GGGTTTACCCGCGAATGCGTCGGTGGTGTCACCCTCGCATTcgc                                                                                                                                        <  2:4534829/150‑1 (MQ=25)                gggCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGC‑GGGTTTACCCGCGAATGCGTCGGTGGTGTCACCCTCGCATTcgc                                                                                                                                        <  2:4534856/150‑1 (MQ=25)                gggCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGC‑GGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTcgc                                                                                                                                        <  1:1569759/150‑1 (MQ=255)                gggCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGC‑GGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTcgc                                                                                                                                        <  1:1407709/150‑1 (MQ=255)                 ggCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGC‑GGGTTTACCCGCGAATGCGTCGGTGGTGTCACCCTCGCATTcgcg                                                                                                                                       <  2:4506945/150‑1 (MQ=25)                 ggCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGC‑GGGTTTACCCGCGAATGCGTCGGTGGTGTCACCCTCGCATTcgcg                                                                                                                                       <  1:1457526/150‑1 (MQ=25)                 ggCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGC‑GGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTcgcg                                                                                                                                       <  1:4721664/150‑1 (MQ=255)                  gCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGC‑GGGTTTACCCGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCgg                                                                                                                                      <  2:4010621/150‑1 (MQ=25)                  gCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGC‑GGGTTTACCCGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCgg                                                                                                                                      <  2:4324355/150‑1 (MQ=25)                  gCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGC‑GGGTTTACCCGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCgg                                                                                                                                      <  1:1336823/150‑1 (MQ=25)                  gCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGC‑GGGTTTACCCGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCgg                                                                                                                                      <  1:3651926/150‑1 (MQ=25)                  gCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGC‑GGGTTTACCCGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCgg                                                                                                                                      <  1:4081558/150‑1 (MQ=25)                  gCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGC‑GGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCgg                                                                                                                                      <  2:2174158/150‑1 (MQ=255)                  gCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGC‑GGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCgg                                                                                                                                      <  2:1662522/150‑1 (MQ=255)                  gCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGC‑GGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCgg                                                                                                                                      <  2:1546091/150‑1 (MQ=255)                  gCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGC‑GGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCgg                                                                                                                                      <  2:3253714/150‑1 (MQ=255)                  gCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGC‑GGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCgg                                                                                                                                      <  2:1289912/150‑1 (MQ=255)                  gCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCAGGGCCCTGTAGGAGC‑GGGTTTACCCGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCgg                                                                                                                                      <  2:2778071/150‑1 (MQ=25)                  gCTGGCGGGCTAGGAGCTGGAGGTGGTGCAGCGCATCCCGGTCGACCCGGGGATGGAGGGCTGAGAACACACACGCAATTGATTTAGCGGGGCCCTGTAGGAGC‑GGGTTTACCCGCGACTGCGTCGGTGGTGTCACCCTCGCATTCGCgg                                                                                                                                      <  2:5288483/150‑1 (MQ=21)                   cTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGC‑GGGTTTACCCGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCggg                                                                                                                                     <  2:2792648/150‑1 (MQ=25)                   cTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGC‑GGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCggg                                                                                                                                     <  1:1543815/150‑1 (MQ=255)                   cTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGC‑GGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCggg                                                                                                                                     <  1:4336183/150‑1 (MQ=255)                   cTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGC‑GGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCggg                                                                                                                                     <  1:2081999/150‑1 (MQ=255)                        gggCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGC‑GGGTTTAAACCCCAAATCGAGGGTGTCAACACCCACGCAATAGCGGGTAAAc                                                                                                                                >  2:2932752/1‑150 (MQ=255)                          gCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGC‑GGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAAccc                                                                                                                              <  1:542088/150‑1 (MQ=255)                                        gtggtgCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGGTTTACCC‑CGCAAGCGAGGGGGGCACCACCCCCGCCTTGGCGGGGAAAACCCCCCCGACACGGcc                                                                                                                >  1:5305598/1‑150 (MQ=255)                                                           tCGACCCGGGGTGTGTTGGCTTATAAAACAAAACCATTTTTTTTTGGGGGGCCGTGTCGGTGG‑GGTTTTACCCGCTAATGCTTGGGTTACACCACCTACGATTTGGGGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAggg                                                                                             <  2:2688981/150‑1 (MQ=255)                                                                    ggatggatGGCTTATAACACACCAGCAATTTATTTAGCGGGGCCGTGTAGGAGC‑GGGTTTACCGGCGATTGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGAt                                                                                    <  2:3109140/150‑1 (MQ=37)                                                                                        acaAGCAATTGATTTAGCGGGGCCCTGTAGGAGC‑GGGTTTACCCGCGAATGCGAGGGTGaca                                                                                                                                                        >  1:2138778/1‑62 (MQ=255)                                                                                        acaAGCAATTGATTTAGCGGGGCCCTGTAGGAGC‑GGGTTTACCCGCGAATGCGAGGGTGaca                                                                                                                                                        <  2:2138778/62‑1 (MQ=255)                                                                                          cagaatttGTTTTGGGGGGGCCGTGTCGGCGC‑GGTTTTACCCGCGAAGGGGGGGGGCCCACCCCCGCCGCTGTCGCGGGAAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCG‑CCCCATACCTCtt                                                             <  2:3225226/144‑1 (MQ=255)                                                                                          aaGCATTTGTTTACGCGGGGCCCGGTCGGAGG‑GGTTTTACCCGCGACTGGGTGGGTGAAACCACCCACGCTTTCGCGGCTAACCCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCG‑CCCCATACCTCtt                                                             <  2:2743145/150‑1 (MQ=255)                                                                                                          ggggCCCTGTAGGAGCGGGGTTTAACC‑CCAATGCGAGGGTGACACCACCGCCCCATTCGCGGGTGAAACCCCGCCTACAGGGCCCGTGGTGGTTTTTTCACGGCGGTGGGATGGCCCGGCGCCCCCAAACCACTGGCCAAAA‑CTTGGAAt                                             >  2:3457639/1‑150 (MQ=255)                                                                                                                      gAGC‑GGGTTTACCCCCGAATTCCACGGTGACACCACCGACGCATTCGCGGGTAAACCCCCTCCTAAACGGCCCGTGGTGCTTTTGACAGGGCGGAGTGGTGGCCCGTCG‑CCCCATAACTCTTGCCAAAATCTTGGAATACCAAAATAAGa                                 >  2:295192/1‑150 (MQ=37)                                                                                                                      gAGC‑GGGTTTAACCGCGAATGCGAGGGTGACACCACCGACCCATTCGCGGGTAAACCCGCTCCTACAGGGCCCGTGGTGCTTTTGACAGGGCGGACTGAAGGCCCGTCG‑CCCCATAACGCGTGCCCAAATCATGGAAAACCCTAATATGa                                 >  2:5163749/1‑150 (MQ=25)                                                                                                                        gC‑GGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCG‑CCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGata                               >  1:4376162/1‑150 (MQ=255)                                                                                                                           gggTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCG‑CCCCATACCTCTTGCCCAAATCTTGGAAGACCATAATATGATATc                             >  2:4859130/1‑150 (MQ=255)                                                                                                                             gTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCG‑CCCCATACCTCTTGCCAAAATCTTGGAATACCataat                                     >  2:4933187/1‑140 (MQ=255)                                                                                                                             gTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCG‑CCCCATACCTCTTGCCAAAATCTTGGAATACCataat                                     <  1:4933187/140‑1 (MQ=255)                                                                                                                              tttACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCG‑CCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCg                          >  2:1660626/1‑150 (MQ=255)                                                                                                                                tACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCG‑CCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTa                        >  2:4139735/1‑150 (MQ=255)                                                                                                                                  cccGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCG‑CCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGa                      >  2:3080822/1‑150 (MQ=25)                                                                                                                                   ccGCGAATGCGTCGGTGGTGTCACCCTCGCCTTCGCGGGTAAACCCGCTCCTACAGGGCCCCTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCG‑CCACATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGAc                     >  2:293802/1‑150 (MQ=25)                                                                                                                                   ccGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCG‑CCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGAc                     >  1:3233512/1‑150 (MQ=25)                                                                                                                                   ccGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCG‑CCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGAc                     >  1:5420455/1‑150 (MQ=25)                                                                                                                                   ccGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCG‑CCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGAc                     >  2:1425925/1‑150 (MQ=25)                                                                                                                                   ccGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCG‑CCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGAc                     >  1:4937659/1‑150 (MQ=25)                                                                                                                                   ccGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCG‑CCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGAc                     >  2:25152/1‑150 (MQ=25)                                                                                                                                   ccGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCACTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCG‑CCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGAc                     >  2:3053425/1‑150 (MQ=25)                                                                                                                                   ccGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCG‑CCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATc                             <  1:3268291/142‑1 (MQ=255)                                                                                                                                   ccGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCG‑CCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATc                             >  2:3268291/1‑142 (MQ=255)                                                                                                                                   ccGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCG‑CCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGAc                     >  2:642461/1‑150 (MQ=255)                                                                                                                                   ccGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCG‑CCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGAc                     >  2:4252136/1‑150 (MQ=255)                                                                                                                                   ccGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCG‑CCCCATACCTCTAGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGAc                     >  2:4764571/1‑150 (MQ=255)                                                                                                                                   ccGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCATTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCG‑CCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGAc                     >  2:1623150/1‑150 (MQ=255)                                                                                                                                    cgcgACTGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCG‑CCCCATACCTCTGGCCAAAATCTTGGAATACCAGAATATGATATCCCGTAGAcc                    >  2:986715/1‑150 (MQ=255)                                                                                                                                    cgcgAATGCGTCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCG‑CCCCATACCTCTTGCCAAAATCTTg                                                 <  2:988486/121‑1 (MQ=21)                                                                                                                                    cgcgAATGCGTCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCG‑CCCCATACCTCTTGCCAAAATCTTg                                                 >  1:988486/1‑121 (MQ=21)                                                                                                                                    cgcgAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCG‑CCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGAcc                    >  1:945952/1‑150 (MQ=255)                                                                                                                                    cgcgAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCG‑CCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGAcc                    >  2:1657187/1‑150 (MQ=255)                                                                                                                                     gcgAATGCGTCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCg                                                                               >  1:479320/1‑91 (MQ=14)                                                                                                                                     gcgAATGCGTCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCg                                                                               <  2:479320/91‑1 (MQ=14)                                                                                                                                     gcgAATGCGTCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCG‑CCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCt                   >  1:2198506/1‑150 (MQ=25)                                                                                                                                     gcgAATGCGTCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCG‑CCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCt                   >  2:4392466/1‑150 (MQ=25)                                                                                                                                     gcgAATGCGTCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCG‑CCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCt                   >  1:5347462/1‑150 (MQ=25)                                                                                                                                     gcgAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCG‑CCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCt                   >  1:3194736/1‑150 (MQ=255)                                                                                                                                       gAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCG‑CCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTg                 >  2:5345427/1‑150 (MQ=255)                                                                                                                                        aaTGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCG‑CCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGa                <  1:5345427/150‑1 (MQ=255)                                                                                                                                          tGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCG‑CCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTCGACCTTGAGt              >  1:221229/1‑150 (MQ=255)                                                                                                                                          tGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCG‑CCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGt              >  1:2224058/1‑150 (MQ=255)                                                                                                                                          tGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCG‑CCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGt              >  2:2488479/1‑150 (MQ=255)                                                                                                                                          gtcggtggTGTCACCCTCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCG‑CCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGt              >  2:3116588/7‑150 (MQ=39)                                                                                                                                                ggTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCG‑CCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGTCTGCAg        >  1:1619151/1‑150 (MQ=255)                                                                                                                                                    acacCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCG‑CCCCATACCTCTTGCCAAAATCTTGGAATACCataa                                      <  1:5290526/116‑1 (MQ=255)                                                                                                                                                    acacCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCG‑CCCCATACCTCTTGCCAAAATCTTGGAATACCataa                                      >  2:5290526/1‑116 (MQ=255)                                                                                                                                                      accaccGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCG‑CCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGTCTGCAGGCCGGc  >  2:5183866/1‑150 (MQ=255)                                                                                                                                                      accaccGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCG‑CCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGTCTGCAGGCCGGc  >  1:1194462/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                       CCGCTGAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGC‑GGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCG‑CCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGTCTGCAGGCCGGC  >  NC_002947/4341901‑4342198

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: ‑