Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,040,623	0	G	A	62.9%	38.7 / 57.7	63	noncoding (107/120 nt)	rrfA	5S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base G (12/11);  new base A (23/16);  total (35/28)
Fisher's exact test for biased strand distribution p-value = 7.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGC  >  NC_000913/4040385‑4040853                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                       tGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGAAACTGCCAGGc                                                                                                                                                                                                                             <  2:468966/250‑1 (MQ=14)  gAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGGAGGAAACTGCCAGGCAt                                                                                                                                                                                                                           >  2:464087/1‑251 (MQ=14)  gAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAAGCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGAAACTGCCAGGCa                                                                                                                                                                                                                            <  2:393617/250‑1 (MQ=14)   aGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTCGTGTGGGGTCTCCTCATGCGAGAGAAGGAAACTGCCAGGCATc                                                                                                                                                                                                                          >  2:182322/1‑251 (MQ=14)          cGGTACTAATGAACCGTGAGGCTTAACCCTACAACGCTGAAGGAGTTTTGGCGGATTGAGAGAAGATTTTCCGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAACCAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGGGTGGGGTCTCCTCATGCGAGAGTAGGAAACTGACAGGCATCAAATaa                                                                                                                                                                                                                    >  2:574875/1‑250 (MQ=11)            ggtcTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTGTCAGCCTGATACAGATTAAAGCAGAACGCAGAAGCGGTCTGATAAAACAGAAGTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGAGGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGAAACTGCCAGGCATCAAATaaacc                                                                                                                                                                                                                 <  2:354585/248‑3 (MQ=255)             tACTAATAAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGAAACTGCCAGGCATCAAATAAAACg                                                                                                                                                                                                                >  1:515798/1‑251 (MQ=18)              aCTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGAAACTGCCAGGCATCAAATAAAACGa                                                                                                                                                                                                               <  2:145254/251‑1 (MQ=21)               cTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACATGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATTCGAGAGTAGGAAACTGCCAGGCATCAAATAAAACGa                                                                                                                                                                                                               >  1:480717/1‑250 (MQ=21)                 aaTGCACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGAAACTGCCAGGCATCAAATAACACGAAAg                                                                                                                                                                                                            <  2:94184/251‑1 (MQ=18)                 aaTGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGAAACTGCCAGGCATCAAATAAAACGAAAg                                                                                                                                                                                                            <  2:187687/251‑1 (MQ=25)                  aTGAACCGTGAGGCTTATCCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGTGGTCTCCTCATGCGAGAGTAGGAAACTGCCAGGCATCAAATAAAACGAAAgg                                                                                                                                                                                                           <  1:269843/251‑1 (MQ=255)                    gAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGGTGTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACTCCGTAGCGCCGATGGTAGTGTGGGGTATCCTCATGCGAGAGTAGGAAACTGCCAGGCAGCAAATAAAACGAAAGGCt                                                                                                                                                                                                         <  2:26466/251‑1 (MQ=255)                     aaCCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTGTCAGCCTGATACAGAGTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATGTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAATTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCTTGCGAGAGTAGGAAACTGCCAGGCATCAAATAAAACGAAAGGCTc                                                                                                                                                                                                        <  1:350503/251‑1 (MQ=255)                        cGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGAAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGt                                                                                                                                                                                                     <  1:295553/251‑1 (MQ=25)                                   cccTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGTGTCTCCTCATGCGAGAGTAGGAAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTgg                                                                                                                                                                                          <  2:357406/250‑1 (MQ=255)                                     cTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGAAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGc                                                                                                                                                                                        >  1:385401/1‑251 (MQ=38)                                         cAACGCCGTAGGTGTTGTGGCGGCTTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATTAAACAGAATTTGCCTGGCGGCAGTAGCGCGGCGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGGAGTGGTGGGTCTCCTCATGCGAGAGGAGGAAACCTCCCGGAATGAAATAAAACGCAAGGCTCAGTCGGACGACTGGGCCttt                                                                                                                                                                                    >  2:142421/1‑251 (MQ=11)                                         cAACGCCGAAGATGTTTTGGCGGAT‑GAGAGAAGCTTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATCAAACAGAATTTGCCTGGCGGCAGTAGCGCGGAGGTCCCACCTGACCCCATGCCGCACTCAGAAGTGAAACGCCGTTGCGCCGCTGGAAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGACGGAAGCCTGGGCCTTTc                                                                                                                                                                                   >  2:333278/1‑251 (MQ=11)                                          aaCGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTTGGGTCTCCTCAGGCGAGAGTAGGAAACGGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTc                                                                                                                                                                                   <  1:586313/251‑1 (MQ=255)                                                 aaGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGAAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTAt                                                                                                                                                                            <  2:357622/251‑1 (MQ=38)                                                      gtTTTGGAGGAT‑GAGAGAAGATTTTCAGCCTGATACTGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTATTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGTCAGCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCtgttg                                                                                                                                                                      <  2:275838/251‑1 (MQ=255)                                                        tttGGCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTAAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGTGCCTTTCGGTTTATCtgttgt                                                                                                                                                                     <  1:127355/250‑1 (MQ=9)                                                           ggCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGttgtttgt                                                                                                                                                                 >  1:17633/1‑251 (MQ=14)                                                                aTTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAAAGCAGACGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCAGAACTCAGACGGGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGAAAATGCCAGGCATCAAATAAAACGAAAGGCTCAGACGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGt                                                                                                                                                             >  1:537563/1‑251 (MQ=11)                                                                     gagaAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGAAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACg                                                                                                                                                        >  1:187957/1‑251 (MQ=18)                                                                      agaagaTTTTCAGCATGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGAAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGc                                                                                                                                                       <  1:284119/251‑1 (MQ=18)                                                                       gaagaTTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCGCAGTCGGAAGACTGGGCCTTTCGTTTTAACTGTTGTTTGTCGGTGAACGct                                                                                                                                                      <  2:298488/251‑1 (MQ=18)                                                                       gaagaTTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAA‑CAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGAAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGctc                                                                                                                                                     >  1:169541/1‑251 (MQ=18)                                                                         agcttttCAGACTGATACAGATTAAATCAGAACGCAGACGAGGTCTGATAAAACAGAATTTGCCTGGAGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAACCGCCGTAGCGCCGATGGTAGTGGGGGGTCTCCTCATGCGAGCGTAGGAAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAAAGctct                                                                                                                                                    >  2:372929/4‑251 (MQ=11)                                                                         agaTTTTCAGCCTGCTACAGATTAAATCAGAACGCAGAAGCGGTCTGATCAAACAGAATTTGCCTGGCGGCAGTAGAGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGCAGTGACACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGAAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGACGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGctct                                                                                                                                                    >  2:521771/1‑251 (MQ=11)                                                                           aTTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTcc                                                                                                                                                  <  2:441105/251‑1 (MQ=33)                                                                               tCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGACGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGAGAGAGTAGGAAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGAGAACGCTCTCCTg                                                                                                                                                >  1:520117/1‑249 (MQ=16)                                                                                cAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTTATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGTCCCCACCTGACCCCATGCCGAACTCAGAAGTTAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGTAACTGCCTTTCTTGAAAGAAATCGTAAGGCTCTTTCGTAGTACTGGGCCTTTGGTTTGATCTGCTTTTGGTCGGTGAACGttcttttgagt                                                                                                                                             <  2:537563/251‑8 (MQ=2)                                                                                cAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGt                                                                                                                                             >  1:421315/1‑251 (MQ=33)                                                                                   ccTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTAGCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAgg                                                                                                                                          <  1:151801/251‑1 (MQ=18)                                                                                      gATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGAc                                                                                                                                        <  1:140941/250‑1 (MQ=33)                                                                                      gATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGAAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGa                                                                                                                                         <  1:294761/249‑1 (MQ=18)                                                                                       aTACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGAAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACa                                                                                                                                       >  2:19688/1‑250 (MQ=18)                                                                                                   tCAGAACGCAGAAGCGGTCTGATAAACCAGAATTTGCCTGGCGGAAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTGTCCTGAGTAGGACAAATCCGCCggg                                                                                                                            >  2:263495/1‑249 (MQ=18)                                                                                                       aaCGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGa                                                                                                                      <  1:142029/251‑1 (MQ=33)                                                                                                              aaGCGGTCTGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTCGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGCACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAAc                                                                                                               >  2:503030/1‑251 (MQ=18)                                                                                                              aaGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCCCAGTCGGAAGACTTGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCCCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAAc                                                                                                               >  1:293565/1‑251 (MQ=18)                                                                                                                gCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGCAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGaaagt                                                                                                             >  1:268218/1‑248 (MQ=18)                                                                                                                 cGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGAAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGtt                                                                                                            >  2:537840/1‑251 (MQ=18)                                                                                                                  ggTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGAAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTg                                                                                                           >  1:99715/1‑251 (MQ=18)                                                                                                                             aCAGAATTTGCCTGGCTGCAGTAGCGCGGCGGTCCCACCTGACCCCCTGCCGAACTCAGAAGTGCAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCCTGCGAGAGTAGGAAACTGCCAGGCATCAAATAAACCGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACgg                                                                                                >  2:258223/1‑251 (MQ=11)                                                                                                                              cAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGACGGGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGATCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGCAGCAACGGc                                                                                               >  2:524628/1‑251 (MQ=16)                                                                                                                                gAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGACGATCTCAGAAGTGAAACGCCGTAGCGCCGATGGTATGGTGGTGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGccc                                                                                             <  1:116976/251‑1 (MQ=255)                                                                                                                                 aaTTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCg                                                                                            >  2:159560/1‑251 (MQ=33)                                                                                                                                       ccTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGCACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGACCGTTGCGAAGCAACGGCCCGGAgggg                                                                                      >  2:391186/1‑250 (MQ=18)                                                                                                                                                           ggtCCCTCCTGACCCCATGCCTAACTCAGAAGTGAAACGCCGTAGCGCCTATGGTAGTGTGGGGTCTCCTCATGCGATAGTAGGGAACTGCCAGGCATCAAATAAAGCGAAAGGCTCAGTCGGAAGACTGTGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCGGAGGAGGACAAATCCGGCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCTCGGAGGGTGGCTGGCTGGAcgcccgcc                                                                   <  2:484832/250‑1 (MQ=11)                                                                                                                                                                       ccccATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGAAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGCAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAgg                                                      >  1:375849/1‑251 (MQ=18)                                                                                                                                                                                gAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGAAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGAGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAAtt                                             >  1:255636/1‑251 (MQ=18)                                                                                                                                                                                 aaCTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGAAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGCAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATCCACTGCCAGGCATCAAAtt                                             >  1:92667/1‑250 (MQ=16)                                                                                                                                                                                                 gccgaagcgcctgtgttAGTGTGGGGTTTCATCATGCGAGTGTGGGGAATTGCAAGGCATAAAAAAAAAGGAAGGGCTCAGTCGGAAGACTGGCCCTTTCTTTTTATCTGTTGTTTGTCGGGGAACGCTCTCTGGAGAAGGACAAATCCGCCGGGAGGGGATTTGAATGTTGCGAAGCAACGCCCGGGGGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCt                            <  1:521347/235‑1 (MQ=2)                                                                                                                                                                                                        cgcCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGAAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGAt                     >  1:101322/1‑251 (MQ=18)                                                                                                                                                                                                         gcCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCCTAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGCCGGATg                    >  2:113072/1‑251 (MQ=18)                                                                                                                                                                                                            gatggttgcgcggagtcacctcATGCGAGCGTAGGGAACTGCCAGGCATCAAATAAGACGAAACGCTCAGTCGGAAGACTGGGCCTTTCGATTTATCTGTTGTTTGTCGGTGAACGATCTCCTGAGTAGGCCAAATCAGCCGGGAGCGGCTTTGAACGTTGCGACGCAACGGCCCGGAGGGCGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCATATTAAGCAGAAAGCCAACCGTACGGATGGcc                 >  2:402993/19‑251 (MQ=11)                                                                                                                                                                                                            gATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGAAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCg                                                                                                         >  2:510793/1‑163 (MQ=25)                                                                                                                                                                                                            gATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGAAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCg                                                                                                         <  1:510793/163‑1 (MQ=25)                                                                                                                                                                                                                   gggtgGGGTCCCCTCATGCGAGAGTAGGGAACTGCCAGGAATAAAAAAAAAAGAAGGGCTAATTCGGAAGCCTGGGCCTTTCGTTTTATGTTTTTTTTGTCGGTGAACGCTCTCCTGAGTAGAACAAATCCGCGGGGGGGGGATTTGAACGTTGGGAAGAAACGGCCCGGGGGGTGGGGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGc          <  1:195107/249‑1 (MQ=11)                                                                                                                                                                                                                              cctcATGCGAGAGTAGGAAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCCTAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGgcgc  >  1:575151/1‑248 (MQ=37)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                       TGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGC  >  NC_000913/4040385‑4040853

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