Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,342,055	0	G	C	36.4%	121.9 / 97.5	89	intergenic (+74/‑257)	ribAB‑II/PP_3814	3,4‑dihydroxy‑2‑butanone‑4‑phosphate synthase/polyamine ABC transportersubstrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (32/24);  new base C (15/17);  total (47/42)
Fisher's exact test for biased strand distribution p-value = 3.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01
Rejected as consensus: Frequency below/above cutoff threshold.

AAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGC‑GGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCG‑CCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGTCTGCAGGCCGGCACCTAC  >  NC_002947/4341907‑4342204                                                                                                                                                      |                                                                                                                                                       aaaTACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGC‑GGGTTTACCCGCGAATGCGAGGGTGACACCACCGa                                                                                                                                                       <  1:2172947/150‑1 (MQ=255)  aaTACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGC‑GGGTTTACCCGCGAATGCGAGGGTGACACCACCGAc                                                                                                                                                      <  1:3604845/150‑1 (MQ=255)  aaTACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGC‑GGGTTTACCCGCGAATGCGAGGGTGACACCACCGAc                                                                                                                                                      <  2:1032520/150‑1 (MQ=255)  aaTACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGC‑GGGTTTACCCGCGAATGCGAGGGTGACACCACCGAc                                                                                                                                                      <  2:1721113/150‑1 (MQ=255)  aaTACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGC‑GGGTTTACCCGCGAATGCGAGGGTGACACCACCGAc                                                                                                                                                      <  2:3766095/150‑1 (MQ=255)  aaTACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGC‑GGGTTTACCCGCGAATGCGAGGGTGACACCACCGAc                                                                                                                                                      <  2:647668/150‑1 (MQ=255)          gggCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGC‑GGGTTTACCCGCGAATGCGTCGGTGGTGTCACCCTCGCATTcgc                                                                                                                                              <  2:4534829/150‑1 (MQ=25)          gggCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGC‑GGGTTTACCCGCGAATGCGTCGGTGGTGTCACCCTCGCATTcgc                                                                                                                                              <  2:4534856/150‑1 (MQ=25)          gggCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGC‑GGGTTTACCCGCGAATGCGTCGGTGGTGTCACCCTCGCATTcgc                                                                                                                                              <  1:782716/150‑1 (MQ=25)          gggCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGC‑GGGTTTACCCGCGAATGCGTCGGTGGTGTCACCCTCGCATTcgc                                                                                                                                              <  1:782749/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:4081558/150‑1 (MQ=25)            gCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGC‑GGGTTTACCCGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCgg                                                                                                                                            <  1:3651926/150‑1 (MQ=25)            gCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGC‑GGGTTTACCCGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCgg                                                                                                                                            <  1:1336823/150‑1 (MQ=25)            gCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGC‑GGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCgg                                                                                                                                            <  2:1546091/150‑1 (MQ=255)            gCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGC‑GGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCgg                                                                                                                                            <  2:3253714/150‑1 (MQ=255)            gCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGC‑GGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCgg                                                                                                                                            <  2:1662522/150‑1 (MQ=255)            gCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGC‑GGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCgg                                                                                                                                            <  2:2174158/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:4336183/150‑1 (MQ=255)             cTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGC‑GGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCggg                                                                                                                                           <  1:2081999/150‑1 (MQ=255)             cTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGC‑GGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCggg                                                                                                                                           <  1:1543815/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)                                                                                    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                                           <  1:4933187/140‑1 (MQ=255)                                                                                                                       gTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCG‑CCCCATACCTCTTGCCAAAATCTTGGAATACCataat                                           >  2:4933187/1‑140 (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                           >  2:25152/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                           >  1:3233512/1‑150 (MQ=25)                                                                                                                             ccGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCACTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCG‑CCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGAc                           >  2:3053425/1‑150 (MQ=25)                                                                                                                             ccGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCG‑CCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATc                                   >  2:3268291/1‑142 (MQ=255)                                                                                                                             ccGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCG‑CCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATc                                   <  1:3268291/142‑1 (MQ=255)                                                                                                                             ccGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCG‑CCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGAc                           >  2:4252136/1‑150 (MQ=255)                                                                                                                             ccGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCG‑CCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGAc                           >  2:642461/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                          >  2:1657187/1‑150 (MQ=255)                                                                                                                              cgcgAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCG‑CCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGAcc                          >  1:945952/1‑150 (MQ=255)                                                                                                                               gcgAATGCGTCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCg                                                                                     <  2:479320/91‑1 (MQ=14)                                                                                                                               gcgAATGCGTCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCg                                                                                     >  1:479320/1‑91 (MQ=14)                                                                                                                               gcgAATGCGTCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCG‑CCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCt                         >  1:2198506/1‑150 (MQ=25)                                                                                                                               gcgAATGCGTCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCG‑CCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCt                         >  1:5347462/1‑150 (MQ=25)                                                                                                                               gcgAATGCGTCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCG‑CCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCt                         >  2:4392466/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                    >  2:2488479/1‑150 (MQ=255)                                                                                                                                    tGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCG‑CCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGt                    >  1:2224058/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)                                                                                                                                                 ccaccGACGCAGTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCG‑CCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGTCTGCAGGCCGGCa       >  2:494634/1‑150 (MQ=255)                                                                                                                                                  caccGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCG‑CCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGTCTGCAGGCCGGCAc      >  1:2140859/1‑150 (MQ=255)                                                                                                                                                   accGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCT‑CCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTATACCTTGAGTCTGCAGGCCGGCAcc     >  2:858879/1‑150 (MQ=255)                                                                                                                                                   accGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCG‑CCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGTCTGCAGGCCGGCAcc     >  1:1290631/1‑150 (MQ=255)                                                                                                                                                     cGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCG‑CCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGTCTGCAGGCCGGCACCTa   >  1:557251/1‑150 (MQ=255)                                                                                                                                                     cGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCG‑CCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGTCTGCAGGCCGGCACCTa   >  2:1800295/1‑150 (MQ=255)                                                                                                                                                      gACGCTTGCGCGGTAAACCCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCG‑CCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGTCTGCAGGCCGGCACCTAc  <  2:1031724/150‑1 (MQ=255)                                                                                                                                                      gACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCG‑CCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGTCTGCAGGCCGGCACCTAc  <  1:3116588/150‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                       AAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGC‑GGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCG‑CCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGTCTGCAGGCCGGCACCTAC  >  NC_002947/4341907‑4342204

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