Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,943,340	C→T	noncoding (1533/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,943,340	0	C	T	100.0%	127.5 / NA	57	noncoding (1533/1542 nt)	rrsC	16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (29/28);  total (29/28)

TTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCAC  >  NC_000913/3943257‑3943473                                                                                    |                                                                                                                                      ttCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATTACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAgg                                                                               >  1:894884/1‑140 (MQ=11) ttCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATTACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAgg                                                                               >  2:367020/1‑140 (MQ=11)   cGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAGGGTAACCGTAGGGGAACCTGCGGTTGGATTACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGaa                                                                             >  2:321471/1‑140 (MQ=17)   cGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATTACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGaa                                                                             >  2:1146093/1‑140 (MQ=17)    gggAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATTACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAg                                                                            >  2:917562/1‑140 (MQ=17)         ggCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATTACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGaaa                                                                       >  2:951163/1‑140 (MQ=17)           cgcTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATTACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAg                                                                     >  2:926086/1‑140 (MQ=17)               tACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATTACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAg                                                                 <  1:926086/140‑1 (MQ=18)                  cACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATTACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCg                                                              <  2:1133763/140‑1 (MQ=18)                   aCTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATTACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGt                                                             >  1:166504/1‑140 (MQ=17)                      ttGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATTACCTCCTTGCATTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAACAAGGCGTCtt                                                          <  2:1178822/140‑1 (MQ=255)                       tgtgATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATTACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTg                                                         >  1:177145/1‑140 (MQ=17)                           aTTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATTACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGaa                                                     <  2:970407/140‑1 (MQ=17)                           aTTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATTACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGaa                                                     <  1:855419/140‑1 (MQ=18)                           aTTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATTACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGaa                                                     <  1:1089722/140‑1 (MQ=18)                           aTTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATTACCTCCTTACCTTAAAGAAGCGTTCTTT‑CAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAg                                                    <  1:972910/140‑1 (MQ=16)                            ttCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATTACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGc                                                        <  1:1091759/136‑1 (MQ=17)                            ttCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATTACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGc                                                        >  2:1091759/1‑136 (MQ=17)                            ttCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATTACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAg                                                    <  1:917562/140‑1 (MQ=18)                            ttCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATTACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAg                                                    <  1:203352/140‑1 (MQ=18)                             tCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATTACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTTCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGc                                                   <  1:277622/140‑1 (MQ=17)                              cATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATTACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCa                                                  <  2:911932/140‑1 (MQ=18)                                tGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATTACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGa                                                <  2:894884/140‑1 (MQ=18)                                tGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATTACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGa                                                <  1:566826/140‑1 (MQ=18)                                 gACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATTACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGa                                                      <  2:896854/133‑1 (MQ=18)                                 gACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATTACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGa                                                      >  1:896854/1‑133 (MQ=17)                                 gACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATTACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGAc                                               >  2:417303/1‑140 (MQ=17)                                      gggTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATTACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATa                                          >  1:1002836/1‑140 (MQ=17)                                      gggTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATTACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATa                                          <  1:417303/140‑1 (MQ=17)                                         tGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATTACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGt                                       <  2:854008/140‑1 (MQ=17)                                         tGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATTACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGt                                       <  2:854009/140‑1 (MQ=17)                                           aaGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATTACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTcc                                     <  1:1146093/140‑1 (MQ=17)                                            aGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATTACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTccc                                    >  2:397698/1‑140 (MQ=17)                                             gTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATTACCCCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTcccc                                   >  2:385368/1‑140 (MQ=17)                                                gTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATTACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGt                                       >  1:891550/1‑133 (MQ=17)                                                gTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATTACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGt                                       <  2:891550/133‑1 (MQ=17)                                                 tAACAAGGTAACCGTAAGGGAACCTGCGGTTGGATTACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCg                               >  1:910989/1‑140 (MQ=17)                                                      aGGTAACCGTAGGGGAACCTGCGGTTGGATTACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTag                          <  2:263051/140‑1 (MQ=17)                                                      aGGTAACCGTAGGGGAACCTGCGGTTGGATTACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTag                          <  1:1050297/140‑1 (MQ=17)                                                      aGGTAACCGTAGGGGAACCTGCGGTTGGATTACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTag                          <  2:166504/140‑1 (MQ=17)                                                      aGGTAACCGTAGGGGAACCTGCGGTTGGATTACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTag                          <  2:1127754/140‑1 (MQ=17)                                                      aGGTAACCGTAGGGGAACCTGCGGTTGGATTACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTag                          <  1:732779/140‑1 (MQ=17)                                                      aGGTAACCGTAGGGGAACCTGCGGTTGGATTACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTag                          <  1:528234/140‑1 (MQ=17)                                                        gTAACCGTAGGGGAACCTGCGGTTGGATTACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTagag                        >  2:1107726/1‑140 (MQ=17)                                                        gTAACCGTAGGGGAACCTGCGGTTGGATTACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTagag                        >  1:1093308/1‑140 (MQ=17)                                                            ccGTAGGGGAACCTGCGGTTGGATTACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTag                          <  2:69387/134‑1 (MQ=17)                                                            ccGTAGGGGAACCTGCGGTTGGATTACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTag                          >  1:69387/1‑134 (MQ=17)                                                               tAGGGGAACCTGCGGTTGGATTACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTag                          >  2:961754/1‑131 (MQ=17)                                                               tAGGGGAACCTGCGGTTGGATTACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTag                          <  1:961754/131‑1 (MQ=17)                                                               tAGGGGAACCTGCGGTTGGATTACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAgg                 >  1:243261/1‑140 (MQ=17)                                                                 ggggAACCTGCGGTTGGATTACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTaga                         <  2:173803/130‑1 (MQ=17)                                                                 ggggAACCTGCGGTTGGATTACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTaga                         >  1:173803/1‑130 (MQ=17)                                                                 ggggAACCTGCGGTTGGATTACCCCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGac               <  1:385368/140‑1 (MQ=17)                                                                  gggAACCTGCGGTTGGATTACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGaca              >  1:658509/1‑140 (MQ=17)                                                                  gggAACCTGCGGTTGGATTACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGaca              >  1:884384/1‑140 (MQ=17)                                                                   ggAACCTGCGGTTGGATTACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGacac             >  2:1089536/1‑140 (MQ=17)                                                                       ccTGCGGTTGGATTACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGcc         >  2:865711/1‑140 (MQ=17)                                                                           cGGTTGGATTACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCttt     >  1:644430/1‑140 (MQ=17)                                                                            ggTTGGATTACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTc    >  1:889400/1‑140 (MQ=17)                                                                              ttGGATTACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCAc  <  1:1089536/140‑1 (MQ=17)                                                                                    |                                                                                                                                      TTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCAC  >  NC_000913/3943257‑3943473

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