Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,942,092	C→A	noncoding (285/1542 nt)	rrsC →	16S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,942,092	0	C	A	100.0%	126.3 / NA	57	noncoding (285/1542 nt)	rrsC	16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (27/30);  total (27/30)

TAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGACACGGTCCAGACTCCTACGGGAGGCAGCAGTGGGGA  >  NC_000913/3941957‑3942170                                                                                                                                        |                                                                               tAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCACTAg                                                                            >  1:1249836/1‑140 (MQ=17)    ctactGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCACTAGCTg                                                                         >  1:446536/1‑140 (MQ=17)    ctactGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCACTAGCTg                                                                         <  1:597789/140‑1 (MQ=17)      actGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCTGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCACTAGCTGGt                                                                       >  2:533658/1‑140 (MQ=17)      actGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCACTAGCTGGt                                                                       >  2:1271984/1‑140 (MQ=17)        tGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCACTAGCTGGTCt                                                                     >  1:182845/1‑140 (MQ=17)        tGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCACTAGCTGGTCt                                                                     <  1:1230364/140‑1 (MQ=17)        tGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCACTAGCTGGTCt                                                                     >  1:178227/1‑140 (MQ=17)        tGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCACTAGCTGGTCt                                                                     >  1:135144/1‑140 (MQ=17)           aaaCGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCACTAGCTGGTCTgag                                                                  <  2:506871/140‑1 (MQ=17)             aCGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCACTAGCTGGTCTGagag                                                                <  2:483612/140‑1 (MQ=17)               ggTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCACTAGCTGGTCTGAGAGGa                                                              >  2:932373/1‑140 (MQ=17)                 tAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCACTAGCTGGTCTGAGAGGATg                                                            >  2:1021688/1‑140 (MQ=17)                 tAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCACTAGCTGGTCTGAGAGGATg                                                            <  1:950686/140‑1 (MQ=17)                 tAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCACTAGCTGGTCTGAGAGGATg                                                            >  2:539849/1‑140 (MQ=17)                 tAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCACTAGCTGGTCTGAGAGGATg                                                            >  2:950686/1‑140 (MQ=17)                   gCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCACTAGCTGGTCTGAGAGGATGAc                                                          <  1:289992/140‑1 (MQ=17)                       aTACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCACTAGCTGGTCTGAGAGGATGACCAGc                                                      <  1:397942/140‑1 (MQ=17)                         aCCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCACTAGCTGGTCTGAGAGGATGACCAGCca                                                    >  2:145882/1‑140 (MQ=17)                           cGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCACTAGCTGGTCTGAGAGGATGACCAGCcaca                                                  >  1:1247212/1‑140 (MQ=17)                             cATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCACTAGCTGGTCTGAGAGGATGACCAGCCACACt                                                >  1:973213/1‑140 (MQ=17)                             cATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCACTAGCTGGTCTGAGAGGATGACCAGCCACACt                                                >  2:318010/1‑140 (MQ=17)                             cATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCACTAGCTGGTCTGAGAGGATGACCAGCCACACt                                                >  1:1357541/1‑140 (MQ=17)                             cATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATCAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCACTAGCTGGTCTGAGAGGATGACCAGCCACACt                                                >  1:1224056/1‑140 (MQ=17)                              aTAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCACTAGCTGGTCTGAGAGGATGACCAGCCACACTg                                               >  1:488786/1‑140 (MQ=17)                              aTAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCACTAGCTGGTCTGAGAGGATGACCAGCCACACTg                                               >  2:681928/1‑140 (MQ=17)                               tAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGT‑GGGTAACGGCTCACCTAGGCGACGATCACTAGCTGGTCTGAGAGGATGACCAGCCACACTGGa                                             >  1:939609/1‑140 (MQ=12)                               tAACGTCGCAAGACCAAAGAGGGGGACCTTCGCGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCACTAGCTGGTCTGAGAGGATGACCAGCCACACTgg                                              <  2:1357541/140‑1 (MQ=17)                                  cGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCACTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAAc                                           >  1:877794/1‑140 (MQ=17)                                     cGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCACTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTga                                        >  1:719151/1‑140 (MQ=17)                                     cGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCACTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTga                                        >  1:144189/1‑140 (MQ=17)                                       cAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCTTCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCACTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTgaga                                      >  1:172706/1‑140 (MQ=17)                                        aaGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCACTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGac                                     <  1:1092633/140‑1 (MQ=17)                                        aaGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCACTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGac                                     <  1:165196/140‑1 (MQ=17)                                         aGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCACTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGaca                                    <  1:509495/140‑1 (MQ=17)                                         aGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCACTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGaca                                    <  2:719151/140‑1 (MQ=17)                                          gACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCACTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGacac                                   <  2:910948/140‑1 (MQ=17)                                              aaaGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCACTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGACACGGTc                               <  2:811255/140‑1 (MQ=17)                                               aaGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCACTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGACACGGTcc                              >  1:785768/1‑140 (MQ=18)                                                agagGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCACTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGACACGGTCCa                             <  1:318010/140‑1 (MQ=17)                                                     gggACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCACTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGACACGGTCCAGACTc                        >  1:144175/1‑140 (MQ=18)                                                        aCCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCACTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGACACGGTCCAGACTCCTa                     <  2:969748/140‑1 (MQ=17)                                                          cTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCACTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGACACGGTCCAGACTCCTACg                   <  2:1134712/140‑1 (MQ=17)                                                          cTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCACTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGACACGGTCCAGACTCCTACg                   <  1:539849/140‑1 (MQ=17)                                                             cGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCACTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGACACGGTCCAGACTCCTACGGGa                <  2:785768/140‑1 (MQ=17)                                                             cGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCACTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGACACGGTCCAGACTCCTACGGGa                <  1:1271984/140‑1 (MQ=17)                                                               ggCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCACTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGACACGGTCCAGACTCCTACGGGAgg              <  1:932373/140‑1 (MQ=17)                                                               ggCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCACTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGACACGGTCCAGACTCCTACGGGAgg              <  1:1325880/140‑1 (MQ=17)                                                                 ccTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGT‑GGGTAACGGCTCACCTAGGCGACGATCACTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGACACGGTCCAGACTCCTACGGGAGgcag           <  2:939609/140‑1 (MQ=12)                                                                    ctTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCACTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGACACGGTCCAGACTCCTACGGGAGgcagca         >  1:663327/1‑140 (MQ=18)                                                                    ctTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCACTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGACACGGTCCAGACTCCTACGGGAGgcagca         <  1:1145092/140‑1 (MQ=17)                                                                    ctTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCACTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGACACGGACCAGACTCCTACGGGAGgcagca         <  2:1015913/140‑1 (MQ=17)                                                                         cttCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCACTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGACACGGTCCAGACTCCTACGGGAGGCAGCAGTggg    <  2:172706/138‑1 (MQ=17)                                                                         cATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCACTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGACACGGTCCAGACTCCTACGGGAGGCAGCAGTggg    <  2:744766/140‑1 (MQ=17)                                                                         cATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCACTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGACACGGTCCAGACTCCTACGGGAGGCAGCAGTggg    <  2:825352/140‑1 (MQ=17)                                                                         cATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCACTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGACACGGTCCAGACTCCTACGGGAGGCAGCAGTggg    <  1:1184289/140‑1 (MQ=17)                                                                           tCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCACTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGACACGGTCCAGACTCCTACGGGAGGCAGCAGTGGGGa  <  1:381104/140‑1 (MQ=17)                                                                                                                                        |                                                                               TAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGACACGGTCCAGACTCCTACGGGAGGCAGCAGTGGGGA  >  NC_000913/3941957‑3942170

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