Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,731,034	G→A	noncoding (124/1542 nt)	rrsG ←	16S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,731,034	0	G	A	100.0%	89.9 / NA	35	noncoding (124/1542 nt)	rrsG	16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (16/19);  total (16/19)

CACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGTTTGATGCTCAAAGAATTAAACTTCGTAATGAATTACGTGTTCACTCTTGAGACTTGGTATTCATTTTTCGTCTTGCGACGTTAAGAATCCGTATCTTCGAGTGCCCACAC  >  NC_000913/2730799‑2731273                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                cACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTATGGCTGGTCATACTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTAACCCACCAACAAGCTACTCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGATACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAAACATTACTCACCCGt                                                                                                                                                                                                                                  >  1:125756/1‑251 (MQ=255) cACTGCTGCCTCCCGTAGGAGTCTGGAACGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGAACAGCTAGGGATCGTAGCCTTGGTGAGAAGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGAAAGAGGACCGAGGGACACCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATACCCATCCATAAGGCAGCTTCCCAAACATTACTCACCCGt                                                                                                                                                                                                                                  >  2:152008/1‑251 (MQ=11)                 aGGAGTCTGGACCGGTTCTCATTTCCAGTGTGGCGGGTCATCCTCTCAGCCCAGCTAGGGAGCGCCGCCTTGGTGGGCCGTTCCCCCACCAACAAGCTACCCCCATCTTGGCACATCCGATGGCTAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCGACCTTTTCCATTAGTTATCCCCCTCCAGCAGGCCGCTTCCCAAACATTACTCACCCGCCCGCCACTCGACAGCg                                                                                                                                                                                                                  <  2:318902/251‑1 (MQ=14)                     gTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTCGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACACGCTAATCCCATCTGGGCACATCCGATGGCACGAGGCCCGAAGGTCCCACTCTTTGGTATTGCGACATTATGCGGAATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGAAGCTTCCCCAACATTACTCACCCGTCCGCCACTCGTCAGCGAagc                                                                                                                                                                                                              >  1:343748/1‑251 (MQ=11)                             cGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAAACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAgct                                                                                                                                                                                                      >  1:299492/1‑251 (MQ=37)                                   tcAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAAACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTcc                                                                                                                                                                                                >  1:44489/1‑251 (MQ=38)                                                       cATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAAACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCAt                                                                                                                                                                            >  2:10484/1‑251 (MQ=38)                                                         tCCTCTCAGACCAGCAAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAAACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCAtgt                                                                                                                                                                          >  1:319417/1‑251 (MQ=255)                                                              tcAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAAACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGtt                                                                                                                                                                       >  1:270554/1‑249 (MQ=39)                                                               cAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAAACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAgg                                                                                                                                                                    <  2:270554/251‑1 (MQ=38)                                                                                ggagcctgggggAGCCGTTCCCCCACCAACAAGCAAACCCCATCTGGACACATCCGATGGCAAGGGGCCCGAAGGTCCCCCACTTTGGTCTTGCGACATTATGCGGTGTTAGCTACCGTTTCCAGAAGTTATCCCCCTCCATCAGGCAGCCTCCCAAACATCGCCCCCCCGGCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGCGTGAGGCCTGCCGCCAGCGTTCa                                                                                                                                                   <  2:33354/240‑1 (MQ=255)                                                                                             gccttaaCCCCACCACCGAGTTAACCCCATCCGGGAACATCCGATGGCAGAAGGCCCAAAGGTCCCCCTCTTTGGTCTTTCGACTTTATGCGGCCTTAGCTACCGTTTCCCTTAGTTACCCCCCTCCAACAGGCAGCGTCCCAAACATTCCTCACCCGTCCGCCACTCGTCAGCGAAGCAGCACGCTTCTTCCTGTTACCGTTCGACTTGCATGTTTTAGGCCCGCCGCCAGCGTTCAATCTGAGCCATGa                                                                                                                                      <  2:255089/245‑1 (MQ=2)                                                                                                          caacaaGCTATTCCCATTCGGGCGCAGCCGAGGGCAGGAGGCCCGAAGGTCCCCCGCTTTGGTCTTGCGCCATTATGCGGTATTAGCTCCCGGTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAAACATTCCTCCCCCGTCCGCCACGCGTCAGCTAAGCATCAAGCTGCTTCCTGTTACCTTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAAGCTGAGTCATGATCAAACCCTTCaa                                                                                                                         <  2:245768/251‑1 (MQ=14)                                                                                                                 cTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGAGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAAACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTaaa                                                                                                                   >  2:156426/1‑250 (MQ=255)                                                                                                                 cTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCAGAAGCTCCCCCTCTTTGTTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCACCCTCCATCAGGCAGCTTCACAAACATTCCTAACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTGAGGCATGCCGCCAGCGTTCACTCTGAGCCATGCTCAAACTCTTCACTTTaaa                                                                                                                   >  1:168480/1‑250 (MQ=14)                                                                                                                   aaTCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAAACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTaaacgt                                                                                                                >  2:258572/1‑248 (MQ=255)                                                                                                                       ccATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAAACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGTTTGa                                                                                                            <  2:51000/251‑1 (MQ=38)                                                                                                                             gggCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTCCGACATTATGCGGTATTAGCTCCCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAAACATTACTCCCCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGTTTGATGCTCa                                                                                                      <  1:318066/251‑1 (MQ=255)                                                                                                                                cacaTCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGAATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAACCACTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGTTTGATGCTCAAAg                                                                                                   <  1:248620/251‑1 (MQ=255)                                                                                                                                                               ccTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAAACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGTTTGATGCTCAAAGAATTAAACTTCGTAATGAATTACGTGTTCAc                                                                    >  1:21462/1‑251 (MQ=37)                                                                                                                                                                  ttttGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTACCCCCCTCCATCAGGCAGCTTCCCAAACATTACCCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGTTTGATGCTCAAAGAATTAAACTTCGTAATGAATTACGTGTTCAct                                                                   <  2:319417/248‑1 (MQ=255)                                                                                                                                                                    ttGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAAACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTc                                                                                                                                                                                     <  1:154935/133‑1 (MQ=255)                                                                                                                                                                    ttGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAAACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTc                                                                                                                                                                                     >  2:154935/1‑133 (MQ=255)                                                                                                                                                                        gcTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCACCTCCATCAGGCAGCTTCCCAAACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGTTTGATGCTCAAAGAATTAAACTTCGTAATGAATTACGTGTTCACTCTTGAGAc                                                           <  1:54578/250‑1 (MQ=37)                                                                                                                                                                                cATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAAACATTACTCACCCGGCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAAGCTGAGCCAGGATCAAACTCTTCAAGTTAAAAGTTTGATGCTCAAAGAATTAAACGTCGTAATGAATTACGTGTTCACTCTTGAGACTTGGTAtt                                                   <  1:325822/251‑1 (MQ=255)                                                                                                                                                                                    aTGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAACCATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGTTTGATGCTCAAAGAATTAAACTTCGTAATGAATTACGTGTTCACTCTTGAGACTTGGTATTCAtt                                               <  2:21462/251‑1 (MQ=21)                                                                                                                                                                                     tGCGGTATTAGCTACCGTTTCCAGTAGTTATCACCCTCCATCAGGAAGCTTCACAAAAATTACTCCCCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGCTCGACTTGCATGTGTTAGGCATGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGTTTGATGCTCAAAGAATTAAACTTCGTAATGAATTACGTGTTCACTCTTGAGACTTGGTATTCAttt                                              >  2:23094/1‑251 (MQ=11)                                                                                                                                                                                                   cctgttCCATTATTTATCCCCCTCCATCAGGCCCCTTCCCAACCATTCCCCCCCCGTCCGCCACTCGTCAGCGCAGCAGCAATCTGCTTCCTGTTACCGTCCGACTTGCATGTGTTAGGCCTGCCGCCAGTGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGTTTGATGCTCAAAGAATTAAACTTCGTAATGAATTACGTGTTCACTCTTGAGACTTGGTATTCATCTTTCGTCTTGCGACg                                <  2:61173/247‑1 (MQ=11)                                                                                                                                                                                                   ccGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAAACATTACGCTCTCGTCCGCCACGCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGTGTGATGCTCAAAGAATTAACCTTCGTAATGAATTACGTGTTCACTCTTGAGACTTGGTATTCATTTTTCGTCTTGCGAcc                                <  2:275659/251‑2 (MQ=11)                                                                                                                                                                                                                  accccccTCCATCAGGCAGCTTCCCAAACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGTTTGATGCTCAAAGAATTAAACTTCGTAATGAATTACGTGTTCACTCGTGAGACTTTGTATTCATTTTTCGTCTTGCGACGTTAAGAATCCGTATc                 <  1:303891/249‑1 (MQ=255)                                                                                                                                                                                                                     ccccTCCATCAGGCGGCTTCCCAAACATTCCCCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCTTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAAGCTGAGCCATGATCAAACTCTTCAATTTAAAAGTTTGATGCTCAAAGAATTAAGCTTCGTAATGAATTACGTGTTCACTCATGGGACGTGGTAGTCAGTTTGCGTCTGTCGACGTTAGGaagctttatcttc              <  2:314380/251‑12 (MQ=11)                                                                                                                                                                                                                     ccccTCCATCAGGCAGCTTCCCAAACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGTTTGATGCTCAAAGAATTAAACTTCGTAATGAATTACGTGTTCACTCTTGAGACTTGGTATTCATTTTTCGTCTTGCGACGTTAAGAATCCGTATCTTc              >  2:35230/1‑251 (MQ=21)                                                                                                                                                                                                                            accAGGCGGCTTCCCAACCATTACTCACCCGGCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGGTAGTCCTGCCGCCAGCGTTCAATCGGAGCCATGATCAAACTCTTCAATTTAAAAGTTTGATGCTCAAAGAATTAAACTTCGTAATGAATTACGGGTTCACTCTTGAGACTTTGTAGTCATTTTTCGTCTTGCGTCGTTAAGAATCCGTATCTTCGAGTGcc       <  2:201240/249‑1 (MQ=11)                                                                                                                                                                                                                             ccAGGCAGCTTCCCAAACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGTTTGATGCTCAAAGAATTAAACTTCGTAATGAATTACGTGTTCACTCTTGAGACTTGGTATTCATTTTTCGTCTTGCGACGTTAAGAATCCGTATCTTCGAGTGccc      <  2:44489/250‑1 (MQ=21)                                                                                                                                                                                                                                 gCAGCTTCCCAAACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTATTACCGTTCGAATTGCATGTGTTAGGCCTGACGCCAGCGTTCAATCTGAGACATGATCAAACTCTTCAATTTAACAGTTTGATGCTCAAAGAATTAACCTTCGTAATGAATTACGTGTTCACTCTTGAGACTTGGTATTCATTTTTCGTCTTGCGACGTTAAGAATCCGTATCTTCGAGTGCCcaccc  >  2:243122/1‑249 (MQ=11)                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                CACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTTAAAAGTTTGATGCTCAAAGAATTAAACTTCGTAATGAATTACGTGTTCACTCTTGAGACTTGGTATTCATTTTTCGTCTTGCGACGTTAAGAATCCGTATCTTCGAGTGCCCACAC  >  NC_000913/2730799‑2731273

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