Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,040,623	0	G	A	50.0%	‑6.0 / 58.7	40	noncoding (107/120 nt)	rrfA	5S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base G (9/11);  new base A (9/11);  total (18/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.49e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

ATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGG‑CTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGA  >  NC_000913/4040379‑4040856                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                           aTGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGAAACTGc                                                                                                                                                                                                                                      <  2:21512/251‑1 (MQ=14)             aaCCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGGTGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCTGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGGCCCTCCTTACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGAAACTGCCAGGCATCaa                                                                                                                                                                                                                            <  1:258638/249‑1 (MQ=11)                  gTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGAAACTGCCAGGCATCAAATAAAAc                                                                                                                                                                                                                     >  1:347262/1‑251 (MQ=18)                           aaCCGTGAGGCTTAACCTTACAACGCCGAAGGTTTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCTGTCTGATAAAACAGAATTGGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCTCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGAAACTGCCAGGCATCAAATAAAACGAAAGG‑CTc                                                                                                                                                                                                           <  1:167218/251‑1 (MQ=255)                                     cTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGCGAGTAGGAAACTGCCAGGCATCAAATAAAACGAAAGG‑CTCAg                                                                                                                                                                                                         >  1:280566/1‑243 (MQ=25)                                     cTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGGCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGAAACTGCCAGGCATCAAATAAAACGAAAGT‑CTCAg                                                                                                                                                                                                         <  2:280566/243‑1 (MQ=255)                                                         gttgtgttGGCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGGCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTACTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGG‑CTCAGTCGGAAGACTGGGCCTTTCGTTTTATCtg                                                                                                                                                                            <  2:175451/245‑1 (MQ=9)                                                                 ggCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGG‑CTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGttgtttgt                                                                                                                                                                    >  1:91304/1‑251 (MQ=14)                                                                    cggatGAGAGAAGATTTTCAGCCTGATACAGATTTAATCAGAACGCATAAGCGGTCTGATAAAACAGAATTTGCCTGTCGGCAGTAGCGCGGTTGTCCCTCCTGACCCCATGCCGAACTGAGAAGTGAAGCGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATGAAATAATACGAAAGG‑CTGATTCGGAAGACTGGGCCGTTCGTTTTATCTGTTGTTTGTCg                                                                                                                                                                  <  2:402947/247‑1 (MQ=11)                                                                      gatGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGG‑CTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGt                                                                                                                                                                <  2:185052/249‑1 (MQ=17)                                                                           gagaAGATTTTCCGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGG‑CTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACg                                                                                                                                                           >  2:20598/1‑251 (MQ=18)                                                                           gagaAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGG‑CTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACg                                                                                                                                                           >  1:185678/1‑251 (MQ=33)                                                                            agaagaTTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGAAACTGCCAGGCATCAAATAAAACGAAAGG‑CTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGc                                                                                                                                                          >  1:401402/1‑251 (MQ=18)                                                                             gaagaTTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGaaa                                                                                                                                                                                                                 >  2:65215/1‑196 (MQ=32)                                                                             gaagaTTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGaaa                                                                                                                                                                                                                 <  1:65215/196‑1 (MQ=32)                                                                                   gttCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCTTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGG‑CTCAGTCGGAAGACTGGGCCTTGCGTGTTATCTGGTGTTTGTCGGTGAACGCTCTCCTg                                                                                                                                                   <  2:400488/250‑1 (MQ=16)                                                                                    tgcAGCCTGATACAGATTAAAGCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGAAACTGCCAGGCATCAAATAAAACGAAAGG‑CTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTg                                                                                                                                                   <  2:236343/248‑1 (MQ=18)                                                                                      cAGCCTGATACAGACTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGATCTGCCAGGCATCAAATAAATCGAAAGG‑CTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGATCGCTCTCCTGAg                                                                                                                                                 <  2:151767/250‑1 (MQ=16)                                                                                         ccTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGG‑CTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAgg                                                                                                                                             <  1:279893/251‑1 (MQ=33)                                                                                           tGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGAAACTGCCAGGCATCAAATAAAACGAAAGG‑CTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGAc                                                                                                                                           >  1:305739/1‑251 (MQ=18)                                                                                                  gATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGG‑CTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATccg                                                                                                                                    >  1:21792/1‑251 (MQ=33)                                                                                                     tAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGAAACTGCCAGGCATCAAATAAAACGAAAGG‑CTCAGTCGGAAGACTGGGCCTTTCGTTTTAGCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATccgcc                                                                                                                                  <  1:114797/250‑1 (MQ=18)                                                                                                          cAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGG‑CTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAg                                                                                                                             <  1:95966/250‑1 (MQ=33)                                                                                                           aGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGAAACTGCCAGGCATCAAATAAAACGAAAGG‑CTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGCGTAGGACAAATCCGCCGGGAGCg                                                                                                                           >  2:43951/1‑251 (MQ=18)                                                                                                               cGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGTGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGG‑CTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGGCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGAtt                                                                                                                       <  1:24705/251‑1 (MQ=18)                                                                                                                    aaGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGAAACTGCCAGGCATCAAATAAAACGAAAGG‑CTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAAc                                                                                                                  <  1:132774/251‑1 (MQ=18)                                                                                                                     aGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGAAACTGCCAGGCATCAAATAAAACGAAAGG‑CTCAGTCGGAAGACTGGGc                                                                                                                                                                                           <  1:307097/177‑1 (MQ=21)                                                                                                                     aGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGAGGTCCCACCTGACCCCAGGACGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGAAACTGCCAGGCATCAAATAAAACGAAAGG‑CTCAGTCGGAAGACTGGGc                                                                                                                                                                                           >  2:307097/1‑177 (MQ=255)                                                                                                                      gCGGTCTGAGAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGTCCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGG‑CTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGt                                                                                                                <  2:220951/251‑1 (MQ=18)                                                                                                                         gTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGG‑CTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGc                                                                                                             >  1:285653/1‑251 (MQ=33)                                                                                                                               tAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGG‑CTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAga                                                                                                        >  1:99233/1‑249 (MQ=18)                                                                                                                                   aCAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGAAACTGCCAGGCATCAAATAAAACGAAAGG‑CTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACgg                                                                                                   <  1:390191/251‑1 (MQ=18)                                                                                                                                     aGAATTTGCCTGGCGGCAGGAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGAAACTGCCAGGCATCAAATAAAACGAAAGG‑CTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGcc                                                                                                 <  1:350538/251‑1 (MQ=18)                                                                                                                                                               gtggtCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGAAACTGCCAGGCATCAAATAAAACGAAAGG‑CTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGAcgcccgc                                                                       <  2:347262/251‑1 (MQ=18)                                                                                                                                                                                cATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGG‑CTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCAt                                                      >  2:84486/1‑251 (MQ=33)                                                                                                                                                                                           cAGAAGAGAAAAGCCGTAGGGCCGCTGGTAGAGTGGCGACTCCTCATGCGAGAGTAGGGAACTACCAGGCATCAAATAACACCACAGG‑CTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGCGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGCGGCGGGCGGGACGACCGCCATAAACTGCCAGGCATCAACTTAAGc                                            >  2:41461/1‑250 (MQ=11)                                                                                                                                                                                              aagggaacccccGTACCGCCGTTGGTGGTGTGGGGCCCCCTCAGGCGAGGGGAGGGACCTGCCAGTCAAAAAAAAAAAC‑AAAGGTCTGTGTCGGAAGTCTGGGCTTTTTGTTTAATTTGTGGTTTGTGGGTGAACGCTCTCGTGAGTATGAAAAATCCGCGGGGAGGGGATTTAAAGGTTGCGAAGCACCGGCCCGGAGGGTGGGGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGaa                                        <  1:267510/241‑1 (MQ=2)                                                                                                                                                                                                      cGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGAAACTGCCAGGCATCAAATAAAACGAAAGG‑CTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAAAGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCCTCACCTACCAAGCATCAAATTAAGCAGAAGGCCCTcc                                >  2:306916/1‑251 (MQ=11)                                                                                                                                                                                                                           gtgGGGTCTCCTCATGCGCGCGTAGGAAACTGCCAGGCATCCAATAAAACGAAAGG‑CTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGCGGCGGGCAGGACGCCCGCCATAAAATGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTgaat           >  2:293646/1‑248 (MQ=11)                                                                                                                                                                                                                                     ctcATGCGAGAGTAGGAAACTGCCAGGCATCAAATAAAACGAAAGG‑CTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGa  >  2:350190/1‑250 (MQ=38)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                           ATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGG‑CTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGA  >  NC_000913/4040379‑4040856

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