Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	90,432	A→G	E113E (GAA→GAG)	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,432	0	A	G	100.0%	129.7 / NA	38	E113E (GAA→GAG)	rsmH	16S rRNA m(4)C1402 methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (14/24);  total (14/24)

CGCGGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAACGTGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTGGCTGCTGCAACGCCGGTGA  >  NC_000913/90193‑90673                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                  ggcgGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTc                                                                                                                                                                                                                                        <  2:106990/250‑1 (MQ=255)             tcacGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGcgcg                                                                                                                                                                                                                             >  2:10937/1‑250 (MQ=255)              cacGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGCTCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGa                                                                                                                                                                                                                            <  1:45856/250‑1 (MQ=255)               acGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATg                                                                                                                                                                                                                          <  1:169665/251‑1 (MQ=255)               acGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTGTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTCTCCGCGATCTTGGCGTCTCTGCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATTCGCGATg                                                                                                                                                                                                                          <  2:262550/251‑1 (MQ=255)                        ccTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTg                                                                                                                                                                                                                  <  1:77740/250‑1 (MQ=255)                               cAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCg                                                                                                                                                                                                          >  1:94317/1‑251 (MQ=255)                                ggCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGGTCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGGTGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCAGTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGt                                                                                                                                                                                                         <  2:217709/250‑1 (MQ=255)                                     ggCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATTCGCGATGGTCCGCTGGACATGCGTATGGa                                                                                                                                                                                                    <  1:212753/251‑1 (MQ=255)                                                      gCTGGCGATCGATCGCGACACGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCAGCGATTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTACGCTGGACCTGCGTATGGACCCAACCCGTGGGCAGt                                                                                                                                                                                   >  2:260852/1‑251 (MQ=255)                                                                 atcgCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGa                                                                                                                                                                         >  2:19781/1‑250 (MQ=255)                                                                             aGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAAcc                                                                                                                                                            <  2:195774/251‑1 (MQ=255)                                                                                  aTCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCCATCATCCACGGACCGTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCa                                                                                                                                                         <  2:11878/249‑1 (MQ=255)                                                                                            cGAAGACTATTGATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGat                                                                                                                                             <  1:158563/251‑1 (MQ=255)                                                                                                          gatCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTAt                                                                                                                                <  2:154756/250‑1 (MQ=255)                                                                                                               gcgcTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGaaaa                                                                                                                          <  1:188851/251‑1 (MQ=255)                                                                                                               gcgcTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTCCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGACTACAAACCGCAGAAGAAGCCGATATAGCCTGGGTATTGaaa                                                                                                                           >  2:242663/1‑250 (MQ=255)                                                                                                                  cTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCt                                                                                                                       >  2:210173/1‑251 (MQ=255)                                                                                                                             tcCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCCGCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAag                                                                                                              >  2:112854/1‑249 (MQ=255)                                                                                                                              ccACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAagag                                                                                                            <  1:89596/250‑1 (MQ=255)                                                                                                                                  ggACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGttt                                                                                                       <  1:122170/251‑1 (MQ=255)                                                                                                                                    aCCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTg                                                                                                     <  1:204625/251‑1 (MQ=255)                                                                                                                                       tttCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGAATTGAAAACCTATGGTGAAGAGCGTTTTGcc                                                                                                   >  1:163553/1‑250 (MQ=255)                                                                                                                                            ccGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCGCCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAACCCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGc                                                                                             >  2:71359/1‑251 (MQ=255)                                                                                                                                            ccGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACg                                                                                              <  2:34675/250‑1 (MQ=255)                                                                                                                                                         aTACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGcc                                                                                 >  1:30646/1‑250 (MQ=255)                                                                                                                                                           aCGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCAtt                                                                              <  2:10738/251‑1 (MQ=255)                                                                                                                                                                       gcgATCTTATCGGCAAGAGCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTAGGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAATCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCTTTTTGCCAAACGCATTGCCCGCGCCATTGTCGACCGTAAc                                                                  <  2:186496/251‑1 (MQ=255)                                                                                                                                                                              tATCGGCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACcg                                                                <  1:10937/246‑1 (MQ=255)                                                                                                                                                                                   gCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCg                                                      >  2:182940/1‑251 (MQ=255)                                                                                                                                                                                     aaGATAGACGGCATTATCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTCTGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTCTTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGAt                                                    >  2:240976/1‑251 (MQ=255)                                                                                                                                                                                           gACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCg                                              >  2:86249/1‑251 (MQ=255)                                                                                                                                                                                                         cGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACt                                  <  1:210173/249‑1 (MQ=255)                                                                                                                                                                                                          gATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGc                               >  1:175092/1‑251 (MQ=255)                                                                                                                                                                                                                       gtcagacGCAACTTGATGATGCTGGGCGCGGGTTTTCCGTTATGCGCGATGGGCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGAAGAAGAAGCCGATATCGCCTGGGTATTGAAACCCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGACCAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTGgc                   <  2:78662/244‑1 (MQ=255)                                                                                                                                                                                                                                aaCTTGATGATGCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTGGCTGCTGCAACg         <  1:38250/251‑1 (MQ=255)                                                                                                                                                                                                                                    tgatgaTGCTGAGCGCGGCTTTTCCTTTAGGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTGGCTGCTGCAACGcc       <  2:94317/249‑1 (MQ=255)                                                                                                                                                                                                                                       ggatgCTGAGCGCGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCTTATGGACCCAACCCGTGGGCAGTCGGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTGGCTGCTGCAACGCCGGTGa  <  2:144801/250‑1 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                  CGCGGTGGTCACTCACGTCTGATCCTCTCGCAGCTTGGCGAAGAGGGGCGTTTGCTGGCGATCGATCGCGACCCGCAGGCTATCGCCGTTGCGAAGACTATTGATGATCCGCGCTTCTCCATCATCCACGGACCTTTCTCCGCGCTGGGCGAATACGTTGCCGAGCGCGATCTTATCGGCAAGATCGACGGCATTCTCCTCGATCTTGGCGTCTCTTCACCGCAACTTGATGATGCTGAACGTGGCTTTTCCTTTATGCGCGATGGTCCGCTGGACATGCGTATGGACCCAACCCGTGGGCAGTCAGCCGCTGAATGGCTACAAACCGCAGAAGAAGCCGATATCGCCTGGGTATTGAAAACCTATGGTGAAGAGCGTTTTGCCAAACGCATTGCCCGCGCCATTGTCGAGCGTAACCGCGAACAGCCGATGACCCGCACCAAAGAACTGGCGGAAGTCGTGGCTGCTGCAACGCCGGTGA  >  NC_000913/90193‑90673

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