Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,337,723	0	A	C	23.4%	133.3 / 10.9	65	G566G (GGT→GGG)	gyrA	DNA gyrase subunit A
Reads supporting (aligned to +/- strand):  ref base A (41/8);  new base C (0/15);  total (41/24)
Fisher's exact test for biased strand distribution p-value = 3.07e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CGCCATGAAGACTTTCACGCCTTCTTCAAACTCGGTCACTGGCAGGATCGCAGTGATACGTTCGTCCTGCTCCAGCGGCAGCAGGTTGACGATCGGACGACCGCGCGCGCCACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACC‑AGCAGTCGGTCGATAAAGTCTT‑CTTCTTTAATACGTGCG‑GCAGATTTACCTTT‑CCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTACGACGTTTGTCACCGAACTGTTCACGAACCAGCTCCAGCTCTTCACGGATCACTTCCATCAGACGATCGGCGCTACCAAGAATACGCAACAGTTCCGCG  >  NC_000913/2337474‑2337966                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                     cGCCATGAAGACTTTCACGCCTTCTTCAAACTCGGTCACTGGCAGGATCGCAGTGATACGTTCGTCCTGCTCCAGCGGCAGCAGGTTGACGATCGGACGACCGCGCGCGCCACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACC‑AGCAGTCGGTCGATAAAGTCTT‑CTTCTTTAATACGTGCG‑GCAGATTTAc                                                                                                                                                                                                                                                     >  1:263360/1‑251 (MQ=255)     aTGAAGACTTTCACGCCTTCTTCAAACTCGGTCACTGGCAGGATCGCAGTGATACGTTCGTCCTGCTCCAGCGGCAGCAGGTTGACGATCGGACGACCGCGCGCGCCACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACC‑AGCAGTCGGTCGATAAAGTCTT‑CTTCTTTAATACGGGCG‑GCAGATTTACCttt                                                                                                                                                                                                                                                 >  1:282065/1‑251 (MQ=255)     aTGAAGACTGTCACGCCTTCTTCAAACTCGGTCACTGGCAGGATCGCAGTGATACGTTCGTCCTGCTCCAGCGGCAGCAGGTTGACGATCGGACGACCGCGCGCGCCACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACC‑AGCAGTCGGTCGATAAAGTCTT‑CTTCTTTAATACGTGCG‑GCAGATTTACCttt                                                                                                                                                                                                                                                 <  2:353019/251‑1 (MQ=255)       gAAGACTTTCACGCCTTCTTCAAACTCGGTCACTGGCAGGATCGCAGTGATACGTTCGTCCTGCTCCAGCGGCAGCAGGTTGACGATCGGACGACCGCGCGCGCCACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACC‑AGCAGTCGGTCGATAAAGGCTT‑CTTCTTTAATACGTGCG‑GCAGAGTTACCTTT‑cc                                                                                                                                                                                                                                              <  2:370337/251‑1 (MQ=255)         aGACTTTCACGCCTTCTTCAAACTCGGTCACTGGCAGGATCGCAGTGATACGTTCGTCCTGCTCCAGCGGCAGCAGGTTGACGATCGGACGACCGCGCGCGCCACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACC‑AGCAGTCGGTCGATAAAGTCTT‑CTTCTTTAATACGTGCG‑GCAGATTTACCTTT‑CCCg                                                                                                                                                                                                                                            <  2:430598/251‑1 (MQ=255)                cACGCCTTCTTCAAACTCGGTCACTGGCAGGATCGCAGTGATACGTTCGTCCTGCTCCAGCGGCAGCAGGTTGACGATCGGACGACCGCGCGCGCCACGAGTGGCTTCCGGCAACTGATACACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACC‑AGCAGTCGGTCGATAAAGTCTT‑CTTCTTTAATACGTGCG‑GCAGATTTACCTTT‑CCCGCCacgac                                                                                                                                                                                                                                     >  1:14919/1‑251 (MQ=255)                cACGCCTTCTTCAAACTCGGTCACTGGCAGGATCGCAGTGATACGTTCGTCCTGCTCCAGCGGCAGCAGGTTGACGATCGGACGACCGCGCGCGCCACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACC‑AGCAGTCGGGCGATAAAGTCTT‑CTTCTTTAATACGTGCG‑GCAGATTTACCTTT‑CCCGCCacgac                                                                                                                                                                                                                                     >  2:123024/1‑251 (MQ=255)                        cttcAAACTCGGTCACTGGCAGGATCGCAGTGATACGTTCGTCCTGCTCCAGCGGCAGCAGGTTGACGATCGGACGACCGCGCGCGCCACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACC‑AGCAGTCGGTCGATAAAGTCTT‑CTTCTTTAATACGTGCG‑GCAGATTTACCTTT‑CCCGCCACGACGCTGCGCt                                                                                                                                                                                                                             >  1:430343/1‑251 (MQ=255)                         ttcAAACTCGGTCACTGGCAGGATCGCAGTGATACGTTCGTCCTGCTCCAGCGGCAGCAGGTTGACGATCGGACGACCGCGCGCGCCACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACC‑AGCAGTCGGTCGATAAAGTCTT‑CTTCTTTAATACGTGCG‑GCAGATTTACCTTT‑CCCGCCACGACGCTGCGCtt                                                                                                                                                                                                                            >  1:192568/1‑251 (MQ=255)                           cAAACTCGGTCACTGGCAGGATCGCAGTGATACGTTCGTCCTGCTCCAGCGGCAGCAGGTTGACGATCGGACGACCGCGCGCGCCACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACC‑AGCAGTCGGACGATAAAGTCTT‑CTTCTTTAATACGTGCG‑GCAGATTTACCTTT‑CCCGCCACGACGCTGCGCTTCg                                                                                                                                                                                                                          >  1:16855/1‑251 (MQ=255)                                  ggTCACTGGCAGGATCGCAGTGATACGTTCGTCCTGCTCCAGCGGCAGCAGGTTGACGATCGGACGACCGCGCGCGCCACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGCGAAGCACAGAATATGGTCGTGAGTGTTCGCCACC‑AGCAGTCGGTCGATAAAGTCTT‑CTTCTTTAATACGTGCG‑GCAGATTTACCTTT‑CCCGCCACGACGCTGCGCTTCGTATTc                                                                                                                                                                                                                     >  2:438527/1‑249 (MQ=255)                                         ggCAGGATCGCAGTGATACGTTCGTCCTCCACCAGCGGCAGCAGATTGACGATCGGACGACCGCGCGCGCCACGAGTGCCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGACTATGGTCGTGAGTGTTCGCCACC‑AGCAGTCGGTCGATAAAGTCTT‑CTTCTTTAATACGTGCG‑GCAGATTTACCTTT‑CCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCgg                                                                                                                                                                                                            >  1:341490/1‑251 (MQ=255)                                            aGGATCGCAGTGATACGTTCGTCCTGCTCCAGCGGCAGCAGGTTGACGATCGGACGACCGCGCGCGCCACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACC‑AGCAGTCGGTCGATAAAGTCTT‑CTTCTTTAATCCGTGCG‑GCAGATTTACCTTT‑CCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTg                                                                                                                                                                                                         >  2:431554/1‑251 (MQ=255)                                             ggATCGCAGTGATACGTTCGTCCTGCTCCAGCGGCAGCAGGTTGACGATCGGACGACCGCGCGCGCCACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACC‑AGCAGTCGGTCGATAAAGTCTT‑CTTCTTTAATACGTGCG‑GCAGATTTACCTTT‑CCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTgt                                                                                                                                                                                                        >  2:216588/1‑250 (MQ=255)                                                     gTGATACGTTCGTCCTGCTCCAGCGGCCGCAGGTTGACGCTAGGACGACCGCGCGCGCCACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATCTGGTCGTGAGTGTTCGCCACC‑AGCAGTCGGTCGATAAAGTCTT‑CTTCTTTAATACGTGCG‑GCAGATTTACCTTT‑CCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAAc                                                                                                                                                                                                >  2:388435/1‑251 (MQ=255)                                                     gTGATACGTTCGTCCTGCTCCAGCGGCAGCAGGTTGACGATCGGACGACCGCGCGCGCCACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGTGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACC‑AGCAGTCGGTCGATAAAGTCTT‑CTTCTTTAATACGTGCG‑GCAGATTTACCTTT‑CCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAAc                                                                                                                                                                                                >  2:235369/1‑251 (MQ=255)                                                     gTGATACGTTCGTCCTGCTCCAGCGGCAGCAGGTTGACGATCGGACGACCGCGCGCGCCACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACC‑AGCAGTCGGTCGATAAAGTCTT‑CTTCTTTAATACGTGCG‑GCAGATTTACCTTT‑CCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAAc                                                                                                                                                                                                >  1:299818/1‑251 (MQ=255)                                                           cgttctccctccTCCGGCGGCAGCAGGTTGACGATCGGCCGCCCGCGCGCGCCACGAGTGGCTTCCGGAAACTGATAAACTTTCATCGAAGAGACGCGCCCCCGGCTGGAGAAGCACAAAATATGGTCGTGAGTGTTCGCCCCC‑AGCAGCCGGTCGATAAAGTCTT‑CTTCTTTAATCCGGGCG‑GCAGTTTTACCTTT‑CCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGcc                                                                                                                                                                                          <  1:297400/240‑1 (MQ=255)                                                                gTCCTGCTCCAGCGGCAGCAGGTTGACGATCGGACGACCGCGCGCGCCACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACC‑AGCAGTCGGTCGATAAAGTCTT‑CTTCTTTAATACGTGCG‑GCAGATTTACCTTT‑CCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGt                                                                                                                                                                                     >  1:285625/1‑251 (MQ=255)                                                                            gggCAGCGGGTTGACGATCGGACGACCGCGCGCGCCACGAGTGGCTTCCGGCAACTGATAATCTTTCAGCGAATAGACGCGACCGCGGCGGGAAAAGCACAAAATATGGTCGTGAGTGTTCCCCACC‑AGCAGTCGGTCGATAAAGTCTT‑CTTTTTTAATACGTGCG‑GAAGATTTCCCTTT‑CCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCa                                                                                                                                                                         <  1:389976/250‑1 (MQ=255)                                                                                    ggTTGACGATAGGACGACCGCGCGCGCCACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACC‑AGCAGTCGGTCGATAAAGTCTT‑CTTCTTTAATACGTGCG‑GCAGATTTACCTTT‑CCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTgagagccgttgaagtgcagc                                                                                                                                                                 >  1:75970/1‑236 (MQ=255)                                                                                     gTTGACGATCGGACGACCGCGCGCGCCACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACC‑AGCAGTCGGTCGATAAAGTCTT‑CTTCTTTAATACGTGCG‑GCAGATTTACCTTT‑CCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACAt                                                                                                                                                                >  1:419086/1‑251 (MQ=255)                                                                                         aCGATCGGACGACCGCGCGCGCCACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACC‑AGCAGTCGGTCGATAAAGTCTT‑CTTCTTTAATACGTGCG‑GCAGATTTACCTTT‑CCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTc                                                                                                                                                            >  1:285954/1‑251 (MQ=255)                                                                                           gATCGGACGCCCGCGCGCGCCGCGATTGGCTTCCGGCATCTGATAAACTTTCATCGAATAGCCGCGCCCACGGCTGGAGAAGCACAGAATAGGGTCGTGAGTGTTCGCCCCC‑AGCTGTCGGTCGATAAAGTCTT‑CTTTTTTAATACGTGCG‑GCAGATTTACCTTT‑CCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCt                                                                                                                                                          <  1:424455/251‑1 (MQ=255)                                                                                               ggACGACCGCGCGCGCCACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACC‑AGCAGTCGGTCGATAAAGTCTT‑CTTCTTTAATACGTGCG‑GCAGATTTACCTTT‑CCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGt                                                                                                                                                      >  1:322982/1‑251 (MQ=255)                                                                                                        gccgcgcCACGGGTGTCTTCCGGCAACTGATACACTTACATCGAAGAGCCGCGCCCGCGGCGGGAGAAGCACAGAATAGGGTGGGGTGTGT‑CCCCACCCGGCAGTCGGTAAATAAAGTCTT‑TTTTTTAAAGAGGGGGG‑GAAGATTTCCCTTT‑CCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATc                                                                                                                                             <  2:280137/249‑1 (MQ=255)                                                                                                               cACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACC‑AGCAGTCGGTCGATAAAGTCTT‑CTTCTTTAATACGTGCG‑GCAGATTTACCTTT‑CCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAgg                                                                                                                                      >  1:235979/1‑251 (MQ=255)                                                                                                                    gTGGCTCCCGGCAACTGATAAACTTTCAACAAATAGACGCGACCGCGGCTGGAGAAGCAAAGATTATGGTGGGGAGTTTCCCCCCCC‑AGCGGTCGGTCTATAAAGTCTT‑CTTTTTTAATACGTGCG‑GCATATTTCCCTTC‑CCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCTCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGAt                                                                                                                                 <  1:367579/251‑1 (MQ=255)                                                                                                                          tcggacAACTAATAAACTTTCATCGAATGGACGCGCCCCCGGCTGAAAAAGCCAAGAATATGTTGGTGAGTGTCCGCCCCA‑AGCAGTCGGTCAATAAGTTCTT‑CTTCTTTAATAGGGGCG‑GCAAATTTCCCTTT‑CCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCt                                                                                                                             <  1:193271/244‑1 (MQ=255)                                                                                                                          tCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACC‑AGCAGTCGGTCGATAAAGTCTT‑CTTCTTTAATACGTGCG‑GCAGATTTACCTTT‑ACCGCCACGACGCTGCGCTTCGTATTCAGATAGCGGCTGATACTTAACGTAGCCCTGGTGCGAGAGCGTCACGAACACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTgc                                                                                                                           >  1:165815/1‑251 (MQ=255)                                                                                                                           cgggCAACGGATAACCTTTCATCAAATGGACGCGACCACGGTGGGAGAAGCAAAGAATAGGGCCGTGAGTGTTCCCCCCC‑AGCAGTCGGTGGATAAAGTCTT‑CTTTTTTAATACGTGCG‑GAAGATTTCCTTTT‑CCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTgcg                                                                                                                          <  1:192544/249‑1 (MQ=255)                                                                                                                               aaagtGATAAACTTCCAACGAAGAGCCGCGACCACGGCGGGAGAAGCACAGAATATGGTCTGGGGTGTTCGCCCCC‑AGAATTCGGTCAATAAGGTCTT‑CTTCTTTAATCCGTGCG‑GCAGATTTACCTTC‑CCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGt                                                                                                                      <  1:375468/247‑1 (MQ=255)                                                                                                                                  cTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACC‑AGCAGTCGGTCGATAAAGTCTT‑CTTCTTTAATACGTGCG‑GCAGATTTACCTTT‑CCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTgg                                                                                                                   >  1:343015/1‑251 (MQ=255)                                                                                                                                    gATAAACTTTCATCGAATAGACGCGCCCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACC‑AGCAGTCGGTCGATAAAGTCTT‑CTTCTTTAATACGTGCG‑GCAGATTTACCTTT‑CCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCg                                                                                                                 >  1:82232/1‑251 (MQ=255)                                                                                                                                        aaCTTTCTCCGAATAGACGCGACCAGGGCTGGAAAAGCACAGAATAGGTTCGTGAGTGTCCGCCACC‑AGCAGTCGGTCGATAAAGTCTT‑TTTTTTTAATACGTGCG‑GCAGATTTCCCTTT‑CCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGa                                                                                                             <  1:106756/251‑1 (MQ=255)                                                                                                                                              cATCGAATAGAAGCGACCACCGCTGGAGAAACACAGAATATGGTAGTGAGAGCTCGCCACA‑AGCAGTCGGTAGATAAAGACTT‑CTTCTTTAATACGTGCG‑GCAGATTTACATTT‑ACCGCCAAGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTACCACCGGTGAGAGAGCGTCACAACCCCATCTTCCGAGGAGCTCAGATCTTCCAGGTTGAGGTCAGCGCTGTTTGCGGGGAtgtc                                                                                                         >  2:361962/1‑246 (MQ=255)                                                                                                                                                    agagccgcgcccgcggcggaagaAGCACAAAATATGGCGGGGAGTGTCCCCCCCC‑ATGAGTCGGTCGAAAAAGTCTT‑TTTTTTTAATACGGGCG‑GAAGATTTCCCTTC‑CCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTacgac                                                                                                 <  1:132313/231‑1 (MQ=255)                                                                                                                                                    aTAGACGCGACCCGGGCTGGAGAAGCACAGAATTTGGTCGGGATTGTCCCCCCCC‑AGCAGTCGGTCGATAAAGTCTT‑CTTTTTTAATCCGTGCG‑GCAGATTTCCCTTT‑CCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTacgac                                                                                                 <  1:191150/251‑1 (MQ=255)                                                                                                                                                    aTAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACC‑AGCAGTCGGTCGATAAAGTCTT‑CTTCTTTAATACGTGCG‑GCAGATTTACCTTT‑CCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTacgac                                                                                                 >  1:30996/1‑251 (MQ=255)                                                                                                                                                            gACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACC‑AGCAGTCGGTCGATAAAGTCTT‑CTTCTTTAATACGTGCG‑GCAGATTTACCTTT‑CCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGCTCTTCCAGGTTGATGTCTGCGCTGTTGGCGGGGGTTTCAGGACGACGTTTGTcc                                                                                         >  2:207079/1‑250 (MQ=255)                                                                                                                                                             aCCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACC‑AGCAGTCGGTCGATAAAGTCTT‑CTTCTTTAATACGTGCG‑GCAGATTTACCTTT‑CCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTAAGTACGACGTTTGTCAc                                                                                        >  2:73227/1‑251 (MQ=255)                                                                                                                                                                aCGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACC‑AGCAGTCGGTCGATAAAGTCTT‑CTTCTTTAATACGTGCG‑GCAGATTTACCTTT‑CCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTACGACGTTTGTCACCg                                                                                      >  1:91338/1‑250 (MQ=255)                                                                                                                                                                         gaAGCACAGAATATGGTCGTGAGTGTTCGCCACC‑AGCAGTCGGTCGATAAAGTCTT‑CTTCTTTAATACGTGCG‑GCAGATTTACCTTT‑CCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTACGACGTTTGTCACCGAACTGTTCAc                                                                            >  1:62184/1‑251 (MQ=255)                                                                                                                                                                          aaGCAAAGAATATGGTCGTGAGTGTACGCCACC‑AGCAGTCGGTCGCTAAAGTCTT‑CTTCTTTAATACGTGCG‑GCTGATTTACCTTT‑CCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTGGATGTCTGCGCTGTTGGCGGTGCTTTCAGTACGACGTTTGTCACCGAACTGATCACg                                                                           >  2:41901/1‑251 (MQ=255)                                                                                                                                                                            gCACAGAATATGGTCGTGAGTGTCCGCCCCC‑AGCAGTCGGTCGATAAATTCTT‑CTTCTTTAATACGTGCG‑GCAGATTTCCCTTT‑CCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTACGACGTTTGTCACCGAACTGTTCACGaa                                                                         <  1:12632/251‑1 (MQ=255)                                                                                                                                                                              acaGAATATGGTCGTGAGTGTTCGCCACC‑AGCAGTCGGTCGATAAAGTCTT‑CTTCTTTAATACGTGCG‑GCAGATTTACCTTT‑CCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTACGACGTTTGTCACCGAACTGTTCACGAAcc                                                                       >  1:79387/1‑251 (MQ=255)                                                                                                                                                                                 gAATATGGTCGTGAGTGTTCGCCACC‑AGCAGTCGGTCGATAAAGTCTT‑CTTCTTTAATACGTGCG‑GCAGATTTACCTTT‑CCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGGTTTCAGTACGACGTTTGTCACCGAACTGTTCACGAAccagc                                                                    >  1:15938/1‑251 (MQ=255)                                                                                                                                                                                 aaacaggtcgttgtttgtTCGCCACC‑AGCCGCAGGTCGATAATTGCTT‑CTTCTTTAATTCGCGCG‑GCAGATTTATCTTT‑CCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCTCGACTACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTACGAAGTTTGTCACCGAACTGTTCACGAAccagc                                                                    <  1:380420/236‑1 (MQ=255)                                                                                                                                                                                     atGGTCGTGAGTGTTCGCCACC‑AGCAGTCGGTCGATAAAGTCTT‑CTTCTTTAATACGTGCG‑GCAGATTTACCTTT‑CCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTACGACGTTTGTCACCGAACTGTTCACGAAccagctccc                                                                >  2:6617/1‑250 (MQ=255)                                                                                                                                                                                     atGGTCGTGAGTGTTCGCCACC‑AGCAGTCGGTCGATAAAGTCTT‑CTTCTTTAATACGTGCG‑GCAGATTTACCTTT‑CCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGCGGTTTCAGGACGTAGCTTGTCACCGAACTGGTCATGGAccagatcc                                                                 >  2:187111/1‑250 (MQ=255)                                                                                                                                                                                        gTCGTGAGTGTTCGCCACC‑AGCAGTCGGTCGATAAAGTCTT‑CTTCTTTAATACGTGCG‑GCAGATTTACCTTT‑CCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTACGACGTTTGTCACCGTACTGTTCACGAAccagctccagct                                                             >  1:50682/1‑251 (MQ=255)                                                                                                                                                                                                 gtTCGCCCCC‑AGCAGTCGGTCGATAAAGTCTT‑CTTCTTAAATGCGTGGG‑GAAGTTTTACCTTT‑CCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTACGACGTTTGTCACCGAACTGTTCACGAACCAGCTCCAGCTCTTCACgg                                                     <  1:153029/250‑1 (MQ=255)                                                                                                                                                                                                  ttCGCCACC‑AGCAGTCGGTCGATAAAGTCTT‑CTTCTTTAATACGTGCG‑GCAGATTTACCTTT‑CCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTACGACGTTTGTCACCGAACCGTTCACGAACCAGCTCCAGCTCCTCACGGAt                                                   >  1:183745/1‑251 (MQ=255)                                                                                                                                                                                                      ccccactccatgcggtcgATAAAG‑CTTGCTTGTTTAATCCTTGCG‑GAAAATTT‑CCTTTACCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTACGACGTCTGTCACCGAACTGTTCACGAACCAGCTCCAGCTCTTCACGGATCAc                                                <  2:55700/239‑1 (MQ=255)                                                                                                                                                                                                            ggcattcggtcaATAAAGTTTT‑CTTCTTTAATAC‑TGGGTGAAGATTTCCCTTC‑CCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTACGACGTTTGTCACCGAACTGTTCACGAACCAGCTCCAGCTCTTCACGGATCACTTCCAt                                          <  1:187846/246‑1 (MQ=255)                                                                                                                                                                                                                 tcggtcgATAAAGTCTT‑CTTCTTTAATACGTGCG‑GCAGATTTACCTTT‑CCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTTGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTCCGACGTTTGTCACCGACCTGTTCCCGAACCAGCTCCAGCTCTTCCCGGATCACTTCCATCAGAc                                     >  2:306521/1‑251 (MQ=255)                                                                                                                                                                                                                  gggtcgATAAAGTCTT‑CTTTTTTAATCGGTGCG‑GAATATTTCCCTCT‑CCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTACGACGTTTGTCACCGAACTGTTCACGAACCAGCTCCAGCTCTTCACGGATCACTTCCATCAGACg                                    <  1:431554/250‑1 (MQ=255)                                                                                                                                                                                                                         tAAAGTCTT‑CTTCTTTAATACGTGCG‑GCAGATTTACCTTT‑CCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGTTCAGATCTTCCAGGTTGATGTCTGAGCTGTTGGCGGTGAGTTCAGTACGACGTTTGTCACCGAACTGTTCACGAACCAGCTCCAGCTCTTAACGGATCACTTCCATCAGACGATCGgcg                             >  1:403278/1‑251 (MQ=255)                                                                                                                                                                                                                                      cttTAATACGTGCG‑GCAGATTTACCTTT‑CCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTACGACGTTTGTCACCGAACTGTTCACGAACCAGCTCCAGCTCTTCACGGATCACTTCCATCAGCCGATCGGCGCTACCAAGAAt                  >  2:88467/1‑250 (MQ=255)                                                                                                                                                                                                                                        ttAATCCGGGCG‑GAAGATTTCCCTTC‑CCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTACGACGTTTGTCACCGAACTGTTCACGAACCAGCTCCAGCTCTTCACGGATCACTTCCATCAGACGATCGGCGCTACCAAGAAt                  <  2:285954/248‑1 (MQ=255)                                                                                                                                                                                                                                             gcGTGCG‑GCAGATTTCCCTTC‑CCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTACGACGTTTGTCACCGAACTGTTCACGAACCAGCTCCAGCTCTTCACGGATCACTTCCATCAGACGATCGGCGCTACCAAGAATACGCAACa          <  1:282773/250‑1 (MQ=255)                                                                                                                                                                                                                                               gTGCG‑GCAGATTTACCTTT‑CCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTACGACGTTTGTCACCGAACTGTTCACGAACCAGCTCCAGCTCTTCACGGATCACTTCCATCAGACGATCGGCGCTACCAAGAATACGCAACAGt        >  1:277522/1‑251 (MQ=255)                                                                                                                                                                                                                                                 ggg‑gCAGTTTTCCCTTC‑CCCGCCACGACGCTGCGCTTCGTATGCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTACGACGTTTGTCACCGAACTGTTCACGAACCAGCTCCAGCTCTTCACGGATCACTTCCAGCAGACGATCGGCGCTACCAAGAAGACGTAACAGTTc      <  2:411129/249‑1 (MQ=255)                                                                                                                                                                                                                                                  gg‑gCATATTTCCCTTT‑CCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTACGACGTTTGTCACCGAACTGTTCACGAACCAGCTCCAGCTCTTCACGGATCACTTCCATCAGACGATCGGCGCTACCAAGAATACGCAACAGTTcc     <  1:123024/250‑1 (MQ=255)                                                                                                                                                                                                                                                      cAGATTTACCTTT‑CCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTACGACGTTTGTCACAGAACTGTTCACGAACCAGCTCCAGCTCTTCACGGATCACTTCCCTCCGACGATCGGCGCTACCAAGAATACGCAACAGTTCcgcg  >  2:115192/1‑251 (MQ=255)                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                     CGCCATGAAGACTTTCACGCCTTCTTCAAACTCGGTCACTGGCAGGATCGCAGTGATACGTTCGTCCTGCTCCAGCGGCAGCAGGTTGACGATCGGACGACCGCGCGCGCCACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACC‑AGCAGTCGGTCGATAAAGTCTT‑CTTCTTTAATACGTGCG‑GCAGATTTACCTTT‑CCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTACGACGTTTGTCACCGAACTGTTCACGAACCAGCTCCAGCTCTTCACGGATCACTTCCATCAGACGATCGGCGCTACCAAGAATACGCAACAGTTCCGCG  >  NC_000913/2337474‑2337966

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