Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,944,115	A→G	noncoding (412/2904 nt)	rrlC →	23S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,944,115	0	A	G	100.0%	165.8 / NA	56	noncoding (412/2904 nt)	rrlC	23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (30/26);  total (30/26)

GAATCAGTATGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑GGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGG  >  NC_000913/3943980‑3944198                                                                                                                                         |                                                                                    gAATCAGTATGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑GGGGGGGCCAt                                                                                 >  1:361135/1‑140 (MQ=37)  aaTCAGTATGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑GGGGGGGCCAtc                                                                                >  1:883482/1‑140 (MQ=37)          tgtgtgTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑GGGGGGGCCATCCTCCAAgg                                                                        <  1:263198/140‑1 (MQ=21)          tgtgtgTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑GGGGGGGCCATCCTCCAAgg                                                                        <  2:282617/140‑1 (MQ=21)              tgtTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑GGGGGGGCCATCCTCCAAGGCTaa                                                                    <  1:346742/140‑1 (MQ=21)                ttAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑GGGGGGGCCa                                                                                  <  1:168669/124‑1 (MQ=25)                ttAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑GGGGGGGCCa                                                                                  >  2:168669/1‑124 (MQ=25)                  aGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑GGGGGGGCCATCCTCCAAGGCTAAATAc                                                                >  2:49669/1‑140 (MQ=21)                  aGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑GGGGGGGCCATCCTCCAAGGCTAAATAc                                                                >  1:840415/1‑140 (MQ=21)                  aGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑GGGGGGGCCATCCTCCAAGGCTAAATAc                                                                >  1:482999/1‑140 (MQ=21)                   gTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑GGGGGGGCCATCCTCCAAGGCTAAATACt                                                               <  2:417668/140‑1 (MQ=21)                    tGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑GGGGGGGCCATCCTCCAAGGCTAAATACTc                                                              >  1:44352/1‑140 (MQ=21)                    tGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑GGGGGGGCCATCCTCCAAGGCTAAATACTc                                                              >  1:396811/1‑140 (MQ=21)                        aGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑GGGGGGGCCATCCTCCAAGGCTAAATACTCctga                                                          <  2:361135/140‑1 (MQ=21)                            tCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑GGGGGGGCCATCCTCCAAGGCTAAATACTCctgactga                                                      >  2:891999/1‑140 (MQ=21)                            tCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑GGGGGGGCCATCCTCCAAGGCTAAATACTCctgactga                                                      >  1:804001/1‑140 (MQ=21)                            tCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑GGGGGGGCCATCCTCCAAGGCTAAATACTCctga                                                          >  2:818721/1‑136 (MQ=21)                            tCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑GGGGGGGCCATCCTCCAAGGCTAAATACTCctga                                                          <  1:818721/136‑1 (MQ=21)                             cTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑GGGGGGGCCATCCTCCAAGGCTAAATACTCCtgactgac                                                     <  1:1079037/140‑1 (MQ=21)                             cTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑GGGGGGGCCATCCTCCAAGGCTAAATACTCCtgactgac                                                     <  2:396811/140‑1 (MQ=21)                             cTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑GGGGGGGCCATCCTCCAAGGCTAAATACTCCtgactgac                                                     <  1:363690/140‑1 (MQ=21)                              tGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑GGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGAcc                                                    <  1:927638/140‑1 (MQ=21)                                   aGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑GGGGGGGTCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAg                                               <  2:458404/140‑1 (MQ=19)                                   aGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑GGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAg                                               <  2:482999/140‑1 (MQ=21)                                     gcgcgcGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑GGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTg                                             >  2:492479/1‑140 (MQ=21)                                     gcgcgcGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑GGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTg                                             <  1:891999/140‑1 (MQ=21)                                     gcgcgcGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑GGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTg                                             >  2:159438/1‑140 (MQ=21)                                        cgcgATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑GGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAAc                                          >  1:1040555/1‑140 (MQ=21)                                        cgcgATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑GGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAAc                                          >  1:801244/1‑140 (MQ=21)                                        cgcgATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGCAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑GGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAAc                                          >  1:768708/1‑140 (MQ=21)                                           gATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑GGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAg                                       <  1:49864/140‑1 (MQ=21)                                                aGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑GGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAAcctatcggt                                  <  2:44352/140‑8 (MQ=255)                                                     gACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑GGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAgg                             >  2:44775/1‑140 (MQ=37)                                                       cAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑GGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTg                                             >  2:201048/1‑122 (MQ=37)                                                       cAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑GGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTg                                             <  1:201048/122‑1 (MQ=25)                                                       cAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑GGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGa                           <  1:460304/140‑1 (MQ=21)                                                       cAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑GGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGa                           >  2:812133/1‑140 (MQ=37)                                                       cAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑GGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGa                           <  2:1040555/140‑1 (MQ=21)                                                       cAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑GGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGa                           >  1:684512/1‑140 (MQ=37)                                                       cAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑GGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGa                           >  2:460304/1‑140 (MQ=37)                                                       cAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑GGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGa                           >  2:862617/1‑140 (MQ=37)                                                       cAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑GGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGa                           >  2:533067/1‑140 (MQ=37)                                                        aGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGGCCGATAGTGAACCAGTACCGTGAGGGa                           <  1:579535/140‑1 (MQ=255)                                                          ccccGTACACAAAAATGCACATATTGTGAGCTCGATGAGCAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑GGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAg                        <  2:768708/140‑1 (MQ=21)                                                            ccGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑GGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGc                      >  1:760146/1‑140 (MQ=37)                                                             cGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑GGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCg                     >  1:124393/1‑140 (MQ=37)                                                              gTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑GGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGa                    >  2:945097/1‑140 (MQ=37)                                                               tACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑GGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGaa                   <  1:862617/140‑1 (MQ=21)                                                                acacAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑GGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGaaa                  >  1:704270/1‑140 (MQ=37)                                                                      aaaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑GGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAAcc            <  2:555671/140‑1 (MQ=21)                                                                      aaaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑GGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAAcc            <  2:684512/140‑1 (MQ=21)                                                                      aaaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑GGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAAcc            >  2:138675/1‑140 (MQ=37)                                                                      aaaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑GGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAAcc            <  1:829530/140‑1 (MQ=21)                                                                         tGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑GGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCg         >  2:281433/1‑140 (MQ=37)                                                                           cacaTATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑GGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGc       >  2:1012390/1‑140 (MQ=25)                                                                                atTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑GGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAggg  <  2:760146/140‑1 (MQ=14)                                                                                                                                         |                                                                                    GAATCAGTATGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑GGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGG  >  NC_000913/3943980‑3944198

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