Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	90,999	T→C	S302S (AGT→AGC)	rsmH →	16S rRNA m(4)C1402 methyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	90,999	0	T	C	98.5%	233.7 / NA	65	S302S (AGT→AGC)	rsmH	16S rRNA m(4)C1402 methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (32/32);  minor base A (1/0);  total (33/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGGCGCTAAAAAGCTCGCTCAACGTGCTGGCCCCGGGTGGGCGGCTTTCGATCATCAGCTTCCACTCGCTGGAAGACCGTATTGTGAAACGTTTTATGCGTGAAAACAGCCGCGGTCCGCAAGTTCCGGCAGGGTTACCGATGACTGAAGAGCAGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGACGGCGGTGACTGTGGTAACCACGGCGCACCATACCCGTTTACTGACCGCTCAGCGCGAACAAC  >  NC_000913/90754‑91239                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                 cAGGCGCTAAAAAGCGCGCTCAACGTGCTGGCCCCGGGTGGGCGGCTTTCGATCATCAGCTTCCACTCGCTGGAAGACCGTATTGTGAAACGTTTTATGCGTGAAAACAGCCGCGGTCCGCAAGTTCCGGCAGGGTTACCGATGACTGAAGAGCAGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTagctcagt                                                                                                                                                                                                                                             <  2:156746/251‑1 (MQ=255)  aGGCGATAAAAAGCTCGCTCAACGTGCTGGCCCCGGGTGGGCGGCTTTCGATCATCAGCTTCCACTCGCTGGAAGACCGTATTGTGAAACGTTTTATGCGTGAAAACAGCCGCGGTCCGCCAGTTCCGGCAGGGTTACCGATGACTGAAGAGCAGCTCAAAAAACTGGGTGGCCGGCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTagctcagtt                                                                                                                                                                                                                                            >  2:19665/1‑251 (MQ=255)   ggCGCTAAAAAGCTCGCTCAACGTGCTGGCCCCGGGTGGGCGGCTTTCGATCATCAGCTTCCACTCGCTGGAAGACCGTATTGTGAAACGTTTGATGCGTGAAAACAGCCGCGGTCCGCAAGTTCCGGCAGGGTTACCGATGACTGAAGAGCAGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCGTAGAACCCTCGTGCCCGTtgctcagttc                                                                                                                                                                                                                                           <  2:538446/251‑1 (MQ=255)          aaaaGCTCGCTCAACGTGCTGGCCCCGGGTGGGCGGCTTTCGATCATCAGCTTCCACTCGCTGGAAGACCGTATTGTGAAACGTTTTATGCGTGAAAACAGCCGCGGTCCGCAAGTTCCGGCAGGGTTACCGATGACTGAAGAGCAGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGCTCAGTTCTGCGTAt                                                                                                                                                                                                                                    >  1:27586/1‑251 (MQ=255)            aaGCTCGCTCAACGTGCTGGCCCCGGGTGGGCGGCTTTCGATCATCAGCTTCCACTCGCTGGAAGACCGTATTGTGAAACGTTTTATGCGTGAAAACAGCCGCGGTCCGCAAGTTCCGGCAGGGTTACCGATGACTGAAGAGCAGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGCTCAGTTCTGCGTATTg                                                                                                                                                                                                                                  <  1:518195/251‑1 (MQ=255)            aaGCTCGCTCAACGTGCTGGCCCCGGGTGGGCGGCTTTCGATCATCAGCTTCCACTCGCTGGAAGACCGTATTGTGAAACGTTTTATGCGTGAAAACAGCCGCGGTCCGCAAGTTCCGGCAGGGTTACCGATGACTGAAGAGCAGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGCTCAGTTCTGCGTATTg                                                                                                                                                                                                                                  <  1:137309/251‑1 (MQ=255)             aGCTCGCTCAACGTGCTGGCCCCGGGTGGGCGGCTTTCGATCATCAGCTTCCACTCGCTGGAAGACCGTATTGTGAAACGTTTTATGCGTGAAAACAGCCGCGGTCCGCAAGTTCCGGCAGGGTTACCGATGACTGAAGAGCAGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGCTCAGTTCTGCGTATTGc                                                                                                                                                                                                                                 <  1:122528/251‑1 (MQ=255)                ccgcccAACGTGCTGGCCCCGGGTGGGCGGCTTTCGATCATCAGCCTCCACTCGCTGGAAGACCGTATTGTGAAACGTTTTATGCGTGAAAACAGCCGCGGTCCGCAAGTTCCGGCAGGGTTACCGATGACTGAAGAGCAGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGCTCAGTTCTGCGTATTGCaga                                                                                                                                                                                                                              <  2:312419/246‑1 (MQ=255)                     cAACGTTCTGGCCCCGGGTGGTCGGCTTTCGATCAACAGCTGCCACTCGCTGGAAGACCGTATTGTGAAACGTTTTATGCGTGAAAACAGCCGCGGTCCGCAAGTTCCGGCAGGGTTACCGATGACTGAAGAGCAGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGCTCAGTTCTGCGTATTGCagagag                                                                                                                                                                                                                           <  2:278018/249‑1 (MQ=255)                       aCGTGCTGGCCCCGGGGGGGCGGCTTTCGATCATCAGCTTCCACTCGCTGGAAGACCGTATTGTGAAACGTTTTATGCGTGAAAACAGCCGCGGTCCGCAAGTTCCGGCAGGGTTACCGATGACTGAAGAGCAGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTAAATGCCGGGCGAAGAAGAGGTGGATGAGAACCCTCGAGCCCGTAGCTCCGTTCTGTGTATTGCAGAGAGGACg                                                                                                                                                                                                                       >  2:437635/1‑251 (MQ=255)                           gCTGGCCCCGGGTGGGCGGTTTTCGATCATCAGCTGCCACTCGCTGGAAGACCGTATTGTGAAACGTTTTATGCGTGAAAACAGCCGCGGTCCGCAAGTTCCGGCAGGGTTACCGATGACTGAAGAGCAGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGCTCAGTTCTGCGTATTGCAGAGAGGACGAAt                                                                                                                                                                                                                    <  2:438264/250‑1 (MQ=255)                                     ggTGGGCGGCTTTCGATCATCAGCTTCCACTCGCTGGAAGACCGTATTGTGAAACGTTTTATGCGTGAAAACAGCCGCGGTCCGCAAGTTCCGGCAGGGTTACCGATGACTGAAGAGCAGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGCTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATcagc                                                                                                                                                                                                         >  1:343626/1‑251 (MQ=255)                                       tGGGCGGCTACCGCTCATCAGCTTCCACTCGCTGGAAGCCCGTATTGTGACACGTTTTATGCGTGAAAACAGCCGCGGTCCGCAAGTTCCGGCAGGGTTACCGATGACTGAAGAGCAGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCCCTAGGCAAGTTAATGCCGGGCGCAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGCTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATcagc                                                                                                                                                                                                         >  2:512053/1‑249 (MQ=255)                                             gcTTTCGATCATCAGCTTCCACTCGCTGGAAGACCGTATTGTGAAACGTTTTATGCGTGAAAACAGCCGCGGTCCGCAAGTTCCGGCAGGGTTACCGATGACTGAAGAGCAGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGCTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGAc                                                                                                                                                                                                 >  1:199641/1‑251 (MQ=255)                                                    atcatcAGCTTCCACTCGCTGGAAGACCGTATTGTGAAACGTTTTATGCGTGAAAACAGCCGCGGTCCGCAAGTTCCGGCAGGGTTACCGATGACTGAAGAGCAGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGCTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGct                                                                                                                                                                                          >  1:419802/1‑251 (MQ=255)                                                      catcaGCTTCCACTCGCTGGAAGACCGTATTGTGAAACGTTTTATGCGTGAAAACAGCCGCGGTCCGCAAGTTCCGGCAGGGTTACCGATGACTGAAGAGCAGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGCTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGctc                                                                                                                                                                                         <  1:118123/250‑1 (MQ=255)                                                      catcaGCTTCCACTAGCTGGAAGACCGTATTGTGAAACGTTTTATGCGTGAAAACAGCCGCGGTCCGCAAGTTCCGGCAGGGTTACCGATGACTGAAGAGCAGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGCTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGctct                                                                                                                                                                                        >  1:204150/1‑251 (MQ=255)                                                          aGCTTCCACTCACTGGAAGACCGAATTGTGAAACGTTTTATGCGTGCAAACAGCCGCGGTCCGCAAGTTACGGCAGGGTTACCGATGACTGAAGAGCAGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACCAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACACTCGTGACCGTAGCTCAGTTCTGCGTATTGCAGAGAGGACGAAGGCATGACCAGCAGAGTGACAGAAGCTCTAAGc                                                                                                                                                                                    >  2:224814/1‑251 (MQ=255)                                                            cTTCCACTCGCTGGAAGACCGTATTGTGAAACGTTTTATGCGTGAAAACAGCCGCGGTCCGCAAGTTCCGGCAGGGTTACCGATGACTGAAGAGCAGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGCTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCa                                                                                                                                                                                   >  1:392064/1‑250 (MQ=255)                                                                    cGCTGGAAGACCTGATTGTGAGACGTTTTATGCGTAAAACCAGCCGCTGTCCGCAAGGTCCGGCAGGGTTACCGATGACTGAAGAGCAGCTCAAAAAACTGGGTGGCTGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGCGCAGTTCTGCGTATTGCGGAGAAGACGAAAGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTaaa                                                                                                                                                                           <  2:357077/250‑1 (MQ=255)                                                                            gACCGTATTGTGAAACGTTTTATGCGTGAAAACAGCCGCGGTCCGCAAGTTCCGGCAGGGTTACCGATGACTGAAGAGCAGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGCTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATg                                                                                                                                                                  <  2:182741/251‑1 (MQ=255)                                                                             aCCGTATTGTGAAACGTTTTATGCGTGAAAACAGCCGCGGTCCGCAAGTTCCGGCAGGGTTACCGATGACTGAAGAGCAGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGCTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGAGAGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATgg                                                                                                                                                                 >  1:345064/1‑251 (MQ=255)                                                                               cGTATTGTGAAACGTTTTATGCGTGAAAACAGCAGCGATCCGCAAGTTTCGCCGGGGTTACCGATGACAGAAGAGCAGCTCAAAAAAGCGGGTAGCCCTCCGCTGCGAGAAATAGGCAAGATAATGCCGGGCGCAGCAGAGGTGGCTGCGAACACCCGTGCCCGTAGCTCAGTTCTGCGTATTGCAGCGAGGACGAATGCCTGATCAGCAGAGTGACCGCAGCTCTAAGCACAGTTAAAGGATCGATggg                                                                                                                                                                >  1:525937/1‑250 (MQ=255)                                                                                    tgtgAAACGTTTTATGCGTGAAAACAGCCGCGGTCCGCAAGTTCCGGCAGGGTTACCGATGACTGAAGAGCAGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGCTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGcc                                                                                                                                                           >  2:95375/1‑250 (MQ=255)                                                                                      tgAAACGGTTTATGCGTGAAAACAGCCGCGGTCCGAAAGTTCCGGAAGGGTTACCGATGACTGAAGAGCAGCTCAAAAAACTGGGTGGCCGTCAGCTGCGATCACTAGGCAAGTTAATGCCGGACGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGCTCAGTTCGGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGCTATAAGCAAGGTTAAAGGATAGATGGGAAGccgcg                                                                                                                                                        >  1:449285/1‑248 (MQ=255)                                                                                        aaccgaattATGCGTGAAAACAGCCTCGGTGCGCAAGTTCGAGCGGGGTTACAGATGACTGAAGAGCAGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGACCGTGGATCAGTTCTGCGAATTGCAGCGGGGACGAACGCATGATCAGCAGAGTGACAGAAACTCTAAGCAAGGTTAAAGGGTCGATGGGACGCCACg                                                                                                                                                        >  2:363329/8‑249 (MQ=255)                                                                                             ttttATGCGTGAAAACCGCCGCGGTCCGCAAGATCAGGCAGGGTTACCGATGACTGAAGAGCCGCTCAAAAAACTGGGTCGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGCTCAGTTCTGCGGATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGATCTAAGCAACGTTAACGGATCGATGGGAAGCCACGAGCGcc                                                                                                                                                  >  2:201158/1‑250 (MQ=255)                                                                                               ttATGCGTGAAAACAGCCGCGGTGCGCAAGTTCCGGCAGGGTTACCGATGACTGAAGAGCAGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGACGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGCTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATg                                                                                                                                               >  1:424968/1‑251 (MQ=255)                                                                                                   gCGTGAAAACAGCCGCGGTACGCAAGTGCCGGCAGGGTTACCGATGACTGAAGAGCAGCTCAAAAAACTGGGGGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGCTCAGTTCTGCATAGTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACATAAGCTCTAAGCAAAGTGAAAGGATCGATGGGAAGCCACGAGCGCCTTGCAtt                                                                                                                                           <  2:194917/251‑1 (MQ=255)                                                                                                    cGTGAAAACAGCCGCGGTCCGCAAGTTCCGGCAGGGTTACCGATGACTGAAGAGCAGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGCTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTg                                                                                                                                          <  1:358861/251‑1 (MQ=255)                                                                                                          aaCAGCCGCGGTCCGCAAGTTCCGGCAGGGTTACCGATGACTGAAGAGCAGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGCTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGgt                                                                                                                                    <  2:473201/251‑1 (MQ=255)                                                                                                             aGCCGCGGTCCGCAAGTTCCGGCAGGGTTACCGATGACTGAAGAGCAGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGCTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGtt                                                                                                                                 <  1:127306/251‑1 (MQ=255)                                                                                                               ccGCGGTCCGCAAGTTCCGGCAGGGTTACCGATGACTGAAGAGCAGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGCTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTAt                                                                                                                               >  1:138454/1‑251 (MQ=255)                                                                                                                 gcgGTCCGCAAGTTCCGGCAGGGTTACCGATGACTGAAGAGCAGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGCTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCg                                                                                                                             <  1:203378/251‑1 (MQ=255)                                                                                                                              tCCGGCAGGGTTACCGATGACTGAAGAGCAGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGCTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCttt                                                                                                                 <  1:120741/250‑1 (MQ=255)                                                                                                                                    aGGGTTACCGATGACTGAAGAGCAGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGCTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGAtt                                                                                                          >  1:12075/1‑251 (MQ=255)                                                                                                                                     gggTTACCGATGACTGAAGAGCAGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGCTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGAttt                                                                                                         <  1:291530/251‑1 (MQ=255)                                                                                                                                        ttACCGATGACTGAAGAGCAGCTCAAAAAACTGGGTGGCCGTCAGCCGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGCTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCCGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGACACGAGCGACATGCATTGCCTGGCGTTATCGGGAACGAGCTTTTGCGATTTggg                                                                                                      >  2:359921/1‑251 (MQ=255)                                                                                                                                              aTGACTGGAGAGCAGCGCAAAAAACCGGGTGGCCGTCAGCTGCGAGCACTAGGAAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGCGCAGGTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAATCGACGTGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTg                                                                                                <  2:458212/251‑1 (MQ=255)                                                                                                                                                   tGAAGAGCAGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCAGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGCTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGATCTAAGCAAAGTTAAAGGATCGCTGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCAct                                                                                           >  2:484297/1‑251 (MQ=255)                                                                                                                                                          cagcTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGCTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTg                                                                                    <  1:19665/251‑1 (MQ=255)                                                                                                                                                          cagcTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGCTCAGTTCTGCGTATTGCAGAGAGGACGAATGCAGGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTGGGGAAGCTGCCACTCTGCCTg                                                                                    <  2:27586/251‑1 (MQ=255)                                                                                                                                                               cAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGTCTGAGAACCCTCGTGCCCGTAGCTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGTGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTc                                                                                 <  1:517507/249‑1 (MQ=255)                                                                                                                                                               cAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGCTCAGTTCTGCGTATTGCAGAGATGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTc                                                                                 <  1:363329/249‑1 (MQ=255)                                                                                                                                                               cAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGCTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCAt                                                                               >  2:31161/1‑251 (MQ=255)                                                                                                                                                                        gggTGGCCGGCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGCTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCattat                                                                      <  2:392064/251‑1 (MQ=255)                                                                                                                                                                         ggTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGCTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCattatt                                                                     <  1:111038/251‑1 (MQ=255)                                                                                                                                                                            ggCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGCTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTg                                                                  >  2:59354/1‑251 (MQ=255)                                                                                                                                                                                  cAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGACGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGCTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTCAGCACAGTTAAAGGATCGATGGGAAGCCACGAGCGCCCTGCATGGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGCAGCTGcc                                                                                              >  2:182401/1‑217 (MQ=255)                                                                                                                                                                                  cAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGCTCAGTTCTGCGTATTGCAGAGATGACGAATGCATGAGCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGcc                                                                                              <  1:182401/217‑1 (MQ=255)                                                                                                                                                                                  cAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGCTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGAcggcg                                                            <  2:105767/251‑1 (MQ=255)                                                                                                                                                                                    gCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGCTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGACGGCGGt                                                          <  2:199641/251‑1 (MQ=255)                                                                                                                                                                                         gAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGCACCCTCGTGCCCGTAGCTCCGTTCTGCGTATTGCAGGGAGGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAgc                                                                                                                                                     >  2:528400/1‑155 (MQ=255)                                                                                                                                                                                         gAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGCTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAgc                                                                                                                                                     <  1:528400/155‑1 (MQ=255)                                                                                                                                                                                              cTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGCTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTg                                                                            <  1:391700/223‑1 (MQ=255)                                                                                                                                                                                              cTAGGCAAGTTAATGCCGGGCGAAGAAGAGGGGGCTGAGAACCCTCGTGCCCGTAGCTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGAATGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGGGTTATCGGTGACGGTCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTg                                                                            >  2:391700/1‑223 (MQ=255)                                                                                                                                                                                                aGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGCTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGACGGCGGTGACCGTGGTAAc                                              >  1:162882/1‑251 (MQ=255)                                                                                                                                                                                                    aaGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGCTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGACGGCGGGGACTGTGGGAACCACg                                          >  2:450001/1‑251 (MQ=255)                                                                                                                                                                                                       ttAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGCTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGACGGCGGTGACTGTGGTAACCACGgcg                                       >  1:371520/1‑251 (MQ=255)                                                                                                                                                                                                           tGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGCTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGACGGCGGTGACTGTGGTAACCACGGCGCAcc                                   >  1:474923/1‑251 (MQ=255)                                                                                                                                                                                                                ggCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGCTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGACGGCGGTGACTGTGGTAACCACGGCGCACCATAcc                              >  1:461466/1‑251 (MQ=255)                                                                                                                                                                                                                 gCGAAGAAGAGGGGGCGTAGAACCCTCGTGCCCGTAGCTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGAGCGATGGGAAGCCACGAGCTCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGGTCATTTGCATTATTTTGACGGCGGTGACTGTGGTAACCACGGCGCATCATAccc                             <  2:70402/251‑1 (MQ=255)                                                                                                                                                                                                                    aagaagAGGTGGCTGAGAACCCTCGTGCCCGTAGCTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGACGGCGGTGACTGTGGTAACCACGGCGCACCCTACCCGtt                          >  1:360764/1‑251 (MQ=255)                                                                                                                                                                                                                    aagaagAGGGGGCTGAGAACCCTCGTGCCCGTAGCTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGACGGCGGTGACTGTGGTAACCACGGCGCACCATACCCGtt                          <  2:204150/251‑1 (MQ=255)                                                                                                                                                                                                                        agagGTGGCTGAGAACCCTCGTGCCCGTAGCTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTg                                                                                                               >  1:130335/1‑162 (MQ=255)                                                                                                                                                                                                                        agagGTGGCTGAGAACCCTCGTGCCCGTAGCTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTg                                                                                                               <  2:130335/162‑1 (MQ=255)                                                                                                                                                                                                                                      aCCCTCGTGCCCGTAGCTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGACGGCGGGGACTGTGGTAACCACGGCGCACCATACCCGTTTACTGACCGCTCAGcgcg        >  1:199510/1‑251 (MQ=255)                                                                                                                                                                                                                                            gTGCCCGTAGCTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGACGGCGGTGACTGTGGTAACCACGGCGCACCATACCCGTTTACTGACCGCTCAGCGCGaacaac  <  1:396547/251‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                 CAGGCGCTAAAAAGCTCGCTCAACGTGCTGGCCCCGGGTGGGCGGCTTTCGATCATCAGCTTCCACTCGCTGGAAGACCGTATTGTGAAACGTTTTATGCGTGAAAACAGCCGCGGTCCGCAAGTTCCGGCAGGGTTACCGATGACTGAAGAGCAGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGACGGCGGTGACTGTGGTAACCACGGCGCACCATACCCGTTTACTGACCGCTCAGCGCGAACAAC  >  NC_000913/90754‑91239

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