Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,056,115	0	A	C	38.3%	73.1 / 11.6	48	G75G (GGT→GGG)	glnL	protein histidine kinase NtrB
Reads supporting (aligned to +/- strand):  ref base A (24/5);  new base C (0/18);  total (24/24)
Fisher's exact test for biased strand distribution p-value = 7.37e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ACGCTTTGCTGAGCAGCTGCGCCGCGCCACGTAAACCGCCAAGCGGATTTTTAATCTCATGTGCCAGGCCGCGCACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCC‑GGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACC‑AGCGTCACTTC‑GTTATCGGTAAAA‑CCTTGCC‑CCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCAGGTTGTCATCGATTAACAAAATACTGTTAATCAGCGAGTTGAGGATCTGCCCAGCATCGGGCTGCGTGCCTGTTGCCATAAAGCAGTCTCCT  >  NC_000913/4055865‑4056352                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                aCGCTTTGCTAAGCAGCTGCGCCTCACCACGTAAACCGCCAAGCGGATTTTTAATCTCATGTGCCAGGCCGCGCACTAAATCACGGGCAGCAACCTGCTGTGCGTTCTTTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATTCCTTCC‑GGCATACGCTTGGCCGTCACAGAAAGGATATGCGATCGCCCGTCGATTACC‑AGCGTCCCTTC‑GTTATCGGTaaaa                                                                                                                                                                                                                                                 <  2:65658/251‑1 (MQ=255)   gCTTTGCTGAGCAGCTGCGCCGCGCCACGTAAACCGCCAAGCGGATTTTTAATCTCATGTGCCAGGCCGCGCACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCC‑GGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACC‑AGCGTCACTTC‑GTTATCGGTAAAA‑cc                                                                                                                                                                                                                                              >  2:31246/1‑251 (MQ=255)                         gcccACGTAAACCGCCAAGCGGTTTTTTAATCTCTGGTGCAAGGCCGCGCACTAAACCGCGGGCAGCATCCTGCTGGGGGTGCTGTAGCTGTTCCTGATTTAAGCGGCGCGGGTTATCCATCGGGGCCATCTCCAGCAGGATCATGCCGTCC‑GGCATGCGCTGGGCCGTCACAGAAAGGATAGGCGAGCGCCCGTCGATGCCC‑GGCTTCATTTC‑GTTAGCG‑TAAAATCCTTGCC‑CCGCCTCCAGACTTTCt                                                                                                                                                                                                                       <  2:447594/249‑1 (MQ=255)                           ccACGTAACCCGCCAAGCGGATTTTTAATCTCATGTGCCGGGCCGCGCACTAAATCACGGGCAGCAACCTGCTGGGCTTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCCCCAGCAGTATCATGCCGCCC‑GGCATACGCTGGCCCGCCACAGAAAGGATAGGCGAGCGCCCGTCGATGACC‑AGCGTCACTTC‑GTTATGGGTAAAC‑CCTTGCC‑CCGCCTCCAGACTTTCTTg                                                                                                                                                                                                                     <  1:416103/251‑1 (MQ=255)                                 aaccccccAAGCGGATTTTTAATCTCATGTGCCAGGCCGCGCACTAAATCACGGGCAGCACCCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGCCGCTGTTTTTCCATCGGAGCCATCTCCAGGAGGATCATGCCGCCC‑GGCATACGCTGGGCCGTCACAGAAAGTATATGCGAGCGCCCGTTGATGACC‑AGCCTCATTTT‑TTTATCGGTAAAA‑CCTTGCC‑CCGCCTCCAGACTTTCTTGCATCAg                                                                                                                                                                                                               <  1:353796/245‑1 (MQ=255)                                  acccGCCAAGCGGATTTTTAATCTCATGTGCCGGGCCGCGCAATAAATCCCGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGGGGCGCTGGTTATCCATCGGAGCCATCCCCAGCAGGATCATGCCGTCC‑GGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGCCC‑AGCGTCACTTC‑GTTACCGGTAAAC‑CCTTGCC‑CCGCCTCCAGACTTTCTTGCATCAGc                                                                                                                                                                                                              <  1:386638/249‑1 (MQ=255)                                      cccAAGCGTATTTTTAATCTAATGTGCCAGGCCGCGCACTAAATCCCGGGCAGCACCCTGCGGGGCGTGCTGTAGCTGTTCCTGATTTAAGCGGCGCTGTTTATCCAGCGGAGCCTTCTCCAGCAGGATCATGCCGCCC‑GGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGCCGATGACC‑AGCGTCACTTC‑GTTATCGGTAAAC‑CCTTGCC‑CCGCCTCCAGACTTTCTTGCATCAGCTCGa                                                                                                                                                                                                          <  1:77903/250‑1 (MQ=255)                                         aaGCGGATTTTTAATCTCATGTGCCAGGCCGCGCACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCC‑GGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACC‑AGCGTCACTTC‑GTTATCGGTAAAA‑CCTTGCC‑CCGCCTCCAGACTTTCTTGCATCAGCTCGatat                                                                                                                                                                                                       >  2:300543/1‑251 (MQ=255)                                                tttttAATCTCATGTGCCAGGCCGCGCACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCC‑GGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACC‑AGCGTCACTTC‑GTTATCGGTAAAA‑CCTTGCC‑CCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATg                                                                                                                                                                                                 >  1:164966/1‑250 (MQ=255)                                                 ttttAATCTCATGTGCCAGGCCGCGCACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGCAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCC‑GGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACC‑AGCGTCACTTC‑GTTATCGGTAAAA‑CCTTGCC‑CCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATg                                                                                                                                                                                                 >  1:114925/1‑249 (MQ=255)                                                      aTCTCATGTGCCAGGCCGCGCACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCC‑GGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACC‑AGCGTCACTTC‑GTTATCGGTAAAA‑CCTTGCC‑CCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGt                                                                                                                                                                                           >  1:404298/1‑250 (MQ=255)                                                          cgggtcccggtcaccgcATTAAACCACGGGAACCAACCTGCTGGGCGGGTTGGAGCTGTTCCTTACTAAGGCGGCGGTGGTTATCCATGGGACCCATCTCCGGCAGTATCATGC‑GTCCCGCCATCCGCTGGGCCGTCACAGAAAGGAAATGCGAGCGCCCGTCGATGCCC‑AGCGTCACTGC‑TTTATCGGTAAAT‑CCTTGCC‑CCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGc                                                                                                                                                                                        <  2:431700/235‑1 (MQ=255)                                                           aTGTGCCAGGCCGCGCACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCC‑GGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACC‑AGCGTCACTTC‑GTTATCGGTAAAA‑CCTTGCC‑CCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTc                                                                                                                                                                                      >  2:189278/1‑250 (MQ=255)                                                             gtgCCAGGCCGCGCACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGATTTAAGCGGCGCTGGTTATCCATCGGGGCCTTCTCCAGCAGGATCATGCCGTCC‑GGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACC‑AGCGTCACTTC‑GTTATCGGTAAAC‑CCTTGCC‑CCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCaa                                                                                                                                                                                    <  1:89458/250‑1 (MQ=255)                                                              ggCCAGGCCGCGCACTAAACCACGGGCAGCAACCTGCGGGGGGTGCTGTAGCTGTTCCTGATTTAAGGGGCGCTGGTTACCCATGGGAGCCATCTCCAGCAGGATCATGCCGTCC‑GGCATACGCTGGGCCGTCACAAAAAGGATATGCGAGCGCCCGTCAATGACC‑AGCGTCACTTC‑TTTATCGGTAAAC‑CCTTGCC‑CCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACa                                                                                                                                                                                  <  1:267673/250‑1 (MQ=255)                                                                                 ccACGGGCAGCACCCTGCTGGGCGTGCTGTAGCTGTCCCTGACTTAGGGGGCGCTGTTTATCCATCGGAGCCATCCCCAGCAGGATCATGCCGCCC‑GGCATCCGCTGGGCCGTCAAAGAAAGGATAGGCGAGCGCCCGTCGATGACC‑AGCGTCTCTTC‑TTTATCGGTAAAC‑CCTTGCC‑CCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTAcc                                                                                                                                                               <  1:319498/250‑1 (MQ=255)                                                                                      ggCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCC‑GGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACC‑AGCGTCACTTC‑GTTATCGGTAAAA‑CCTTGCC‑CCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGACGTAGCTCAACAGTTCCGGTAACGGTGTACCAAAc                                                                                                                                                           >  2:424631/1‑250 (MQ=255)                                                                                      ggCAGCAACCTGCTGGGAGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCC‑GGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGCTGACC‑AGCGTCACTTC‑GTTATCGGTAAAA‑CCTTGCC‑CCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAAcc                                                                                                                                                          >  2:339878/1‑250 (MQ=255)                                                                                           caACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCC‑GGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACC‑AGCGTCACTTC‑GTTATCGGTAAAA‑CCTTGCC‑CCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTg                                                                                                                                                     >  2:67792/1‑251 (MQ=255)                                                                                                 gctgGGCGTGCTGTAGCTGTTCCTGACTAAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGAATCATGCCGTCC‑GGCATACGCGGGGCCGTCACAAAAAGGATATGCGAGCGCCCGTCGATGACC‑AGCGTCACTTC‑TTTATCGGTAAAC‑CCTTGCC‑CCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGa                                                                                                                                                <  1:147691/250‑1 (MQ=255)                                                                                                             gaagttGTTCCTGCTTTAGGCGGCGCTGTTTATCCTTCGGACCCTTCTCCAGCAGTATCATGCCGTCG‑GGAATCCTCGGGGCCGTCACAGAAGGGATATGCGAGCGCCCGTCTATGACC‑AGCGTCTTTTC‑GTTATCGGTAAAC‑CTTTGCC‑CCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGc                                                                                                                                   <  1:402091/246‑1 (MQ=255)                                                                                                               aGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCC‑GGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACC‑AGCGTCACTTC‑GTTATCGGTAAAA‑CCTTGCC‑CCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGc                                                                                                                                   >  1:387995/1‑249 (MQ=255)                                                                                                                   gTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCC‑GGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACC‑AGCGTCACTTC‑GTTATCGGTTAAA‑CCTTGCC‑CCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGTGCAgttgt                                                                                                                             >  2:204784/1‑251 (MQ=255)                                                                                                                         gACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCC‑GGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACC‑AGCGTCACTTC‑GTTATCGGTAAAA‑CCTTGCC‑CCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGcgcg                                                                                                                       >  1:201313/1‑251 (MQ=255)                                                                                                                              aaGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCC‑GGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACC‑AGCGTCACTTC‑GTTATCGGTAAAA‑CCTTGCC‑CCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAgg                                                                                                                  >  1:347180/1‑251 (MQ=255)                                                                                                                                                   ggAGCCATCTCCAGCAGGATCATGCCGTCC‑GGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACC‑AGCGTCACTTC‑GTTATCGGTAAAA‑CCTTGCC‑CCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGcc                                                                                               >  2:300566/1‑249 (MQ=255)                                                                                                                                                          cccccAGCATGATCATGCCCTCC‑GGCATCCGTTGGCCCGTCCCAGAAAGGATAGGCCAGCGCCCGGCATTGACC‑ACCGTACCTTC‑TTTATCATTAAAC‑CCTGGCC‑CCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCAGGTTGTCa                                                                                      <  1:173746/248‑1 (MQ=255)                                                                                                                                                               agcagGATCATGCCGTCC‑GGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACC‑AGCGTCACTTC‑GTTATCGGTAAAA‑CCTTGCC‑CCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCAGGTTGTCATCGAt                                                                                 >  2:328677/1‑251 (MQ=255)                                                                                                                                                                gcagGATCATGCCGCCC‑GGCATCCGGTGGGCCGTCAAAGAAAGGATATGCGGGCGCCCGTCGATGACC‑AGCGTCATTTC‑TTTATCGGTAAAC‑CCTGGCC‑CCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCAGGTTGTCATCGAtt                                                                                <  1:378292/251‑1 (MQ=255)                                                                                                                                                                 gggggttgtgcccgccgggCATACGCGGGGCCGACAAAGAAGGGATATGCGGGCCCCCGGCGTTGACC‑GGTCTCATTTC‑GTTACCGGTAAAC‑CCTTGCC‑CCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCAGGTTGTCATCGAtt                                                                                <  2:450204/234‑1 (MQ=255)                                                                                                                                                                  agGATCATGCCGTCC‑GGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACC‑AGCGTCACTTC‑GTTATCGGTAAAA‑CCTTGCC‑CCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCAGGTTGTCATCGATTa                                                                               >  1:153307/1‑250 (MQ=255)                                                                                                                                                                             gccc‑GGCATACGCTGGGCCGTCACAGAAAGGATATGCGCGCGCCCGCCGATGACC‑AGCGTCACTTC‑TTTATCGGAAAAC‑CCTTCCC‑CCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCAGGTTGTCATCGATTAACAAAATACTGt                                                                   <  1:64352/249‑1 (MQ=255)                                                                                                                                                                                c‑GGCATCCGCTGGGCCGCCACAGAAAGGATATGCGAGCGCCCGTCGATGACC‑GGCGCCCTTTC‑GTTATCGGTAAAA‑CCTTGCC‑CCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCAGGTTGTCATCGATTAACAAAATACTGTTaa                                                                <  1:66923/251‑1 (MQ=255)                                                                                                                                                                                c‑GGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACC‑AGCGTCACTTC‑GTTATCGGTAAAA‑CCTTGCC‑CCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCAGGTTGTCATCGATTAACAAAATACTGTTa                                                                 >  2:8077/1‑250 (MQ=255)                                                                                                                                                                                  ggCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACC‑AGCGTCACTTC‑TTTATCGGTAAAC‑CCTTGCC‑CCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAg                                                                                                                                    <  1:213116/182‑1 (MQ=255)                                                                                                                                                                                  ggCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACC‑AGCGTCACTTC‑GTTATCGGTAAAA‑CCTTGCC‑CCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAg                                                                                                                                    >  2:213116/1‑182 (MQ=255)                                                                                                                                                                                           gggtcctacccaaaaaGGATATGGGCGCGCCCGTCGATGACC‑AGCGTCCCTTC‑TTTATCGTTAAAC‑CCTTGCC‑CCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCAGGTTGTCATCGATTAACAAAATACTGTTAATCAGCGAGtt                                                      <  1:231723/238‑1 (MQ=255)                                                                                                                                                                                           ggggCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACC‑AGCGTCTCTTC‑TTTATCGGTAAAC‑CCTTGCC‑CCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCAGGTTGTCATCGATTAACAAAATACTGTTAATCAGCGAGtt                                                      <  1:443666/250‑1 (MQ=255)                                                                                                                                                                                                  ccacaGAAGGGATATGCGACCCCCCGCCG‑TGACCCAGCGTCTCTTC‑TTTACCGGTAAAC‑CCTTGCC‑CCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCAGGTTGTCATCGATTAACAAAATACTGTTAATCAGCGAGTTGAGGATc                                               <  1:71727/250‑1 (MQ=255)                                                                                                                                                                                                      aGAAAGGATATGCGAGCGCCCGTCGATGACC‑AGCGTCACTTA‑GTTATCGGTAAAA‑CCTTGCC‑CAGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTACCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCAGGTTGTCATCGCTTAACAAAATACTGTTAATCAGCGAGTTGAGGATCTGcc                                           >  2:271121/1‑251 (MQ=255)                                                                                                                                                                                                         aaGGTTATGCGGGCGCCCTCCGATCACC‑AGGGTC‑CTTCCTTTACCGGAAAAC‑CCTTCCC‑CCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCAGGTTGTCATCGATTAACAAAATACTGTTAATCAGCGAGTTGAGGATCTGCCCAg                                        <  1:424631/251‑1 (MQ=255)                                                                                                                                                                                                          aGGAGATGCGAGCGCCCGTCGATGAAC‑AGCGTCACTTC‑GTTATCGGTAAAA‑CCTTGCC‑CCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCAGGTGGTCATCGATTAACAAAATACTGTTAATCAGCGAGTTGAGGCTCTGCCCAGc                                       >  1:336541/1‑251 (MQ=255)                                                                                                                                                                                                                       gcCCGGCGATGACC‑AGCGTCACTTC‑GTTATCGGTAAAA‑CCTTGCC‑CCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGCAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCAGGTTGTCATCGATTAACAAAATACTGTTAATCAGCGAGTTGAGGATCTGCCCAGCATCGGGCTGCGTg                          >  1:370110/1‑251 (MQ=255)                                                                                                                                                                                                                                      ggCGTCACTTC‑GTTATCGGTAAAA‑CCTTGCC‑CCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCAGGTTGTCATCGATTAACAAAATACTGTTAATCAGCGAGTTGAGGATCTGCCCAGCATCGGGCTGCGTGCCTGTTGCCATaaa            <  1:303680/250‑1 (MQ=255)                                                                                                                                                                                                                                       gCGTCACTTC‑GTTATCGGTAAAA‑CCTTGCC‑CCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCAGGTTGTCATCGATTAACAAAATACTGTTAATCAGCGAGTTGAGGATCTGCCCAGCATCGGGCTGCGTGCCTGTTGCCATaa             >  2:374180/1‑249 (MQ=255)                                                                                                                                                                                                                                        cGTCACTTC‑GTTATCGGTAAAA‑CCTTGCC‑TCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCAGGTTGTCATCGATTAACAAAATACTGTTAATCAGCGAGTTGAGGATCTGCCCAGCATCGGGCTGCGTGCCTGTTGCCATAAAGc          >  1:50376/1‑251 (MQ=255)                                                                                                                                                                                                                                               tC‑GTTATCGGTAAAA‑CCTTGCC‑CCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCAGGTTGTCATCGATTAACAAAATACTGTTAATCAGCGAGTTGAGGATCTGCCCAGCATCGGGCTGCGTGCCTGTTGCCATAAAGCAGTCTcc   >  2:416655/1‑251 (MQ=255)                                                                                                                                                                                                                                                tcgatctcgttAAAA‑CCT‑GCCTCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTATCTCAACAGTTCCGGTAACGGTGTACCAAACAAGTTGCGGGAGCGTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCAGGTTGTCATCGATTAACAAAATACTGTTAATCAGCGAGGTGAGGATCTGCCCAGCATCGGGCTGCGTGCCTGTTGCCATAAAGCAGtctc    <  1:385330/240‑1 (MQ=255)                                                                                                                                                                                                                                                   ttATCGGTAAAC‑CCTTGCC‑CCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCAGGTTGTCATCGATTAACAAAATACTGTTAATCAGCGAGTTGAGGATCTGCCCAGCATCGGGCTGCGTGCCTGTTGCCATAAAGCAGTCTCCt  <  1:189278/249‑1 (MQ=255)                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                ACGCTTTGCTGAGCAGCTGCGCCGCGCCACGTAAACCGCCAAGCGGATTTTTAATCTCATGTGCCAGGCCGCGCACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCC‑GGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACC‑AGCGTCACTTC‑GTTATCGGTAAAA‑CCTTGCC‑CCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCAGGTTGTCATCGATTAACAAAATACTGTTAATCAGCGAGTTGAGGATCTGCCCAGCATCGGGCTGCGTGCCTGTTGCCATAAAGCAGTCTCCT  >  NC_000913/4055865‑4056352

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑