Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,037,929	(G)6→7	noncoding (411/2905 nt)	rrlA →	23S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,037,923	1	.	G	100.0%	90.7 / NA	42	noncoding (405/2905 nt)	rrlA	23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (21/21);  total (21/21)

CCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑GGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAA  >  NC_000913/4037786‑4037972                                                                                                                                           |                                                  ccAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑gg                                                 <  2:987592/140‑1 (MQ=32) ccAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑gg                                                 >  1:49573/1‑140 (MQ=32)   agagCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑gggg                                               >  2:228181/1‑140 (MQ=32)    gagCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑ggggg                                              <  2:49573/140‑1 (MQ=32)     agCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑gggggg                                             >  2:735131/1‑140 (MQ=33)     agCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑gggggg                                             >  2:1102484/1‑140 (MQ=33)     agCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑gggggg                                             <  2:1589394/140‑1 (MQ=32)          gAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑gggggggacca                                        >  2:1326124/1‑135 (MQ=17)          gAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑gggggggacca                                        >  1:434329/1‑135 (MQ=18)          gAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑gggggggacca                                        >  2:1012354/1‑135 (MQ=17)          gAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑gggggggacca                                        >  2:1236826/1‑135 (MQ=17)          gAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑gggggggacca                                        <  1:1390368/140‑6 (MQ=17)          gAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑gggggggacca                                        <  1:1102484/140‑6 (MQ=17)               aGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCAtcctc                                    >  1:1605989/1‑140 (MQ=17)               aGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCAtcctc                                    >  2:520041/1‑140 (MQ=18)               aGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCAtcctc                                    <  1:49838/140‑1 (MQ=17)                gtgtgtgtgtTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCAtcctcc                                   >  1:1037413/1‑140 (MQ=18)                 tgtgtgtgtTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCa                                  <  2:1279205/140‑1 (MQ=17)                 tgtgtgtgtTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCa                                  <  1:748571/140‑1 (MQ=17)                    gagtgtTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAgg                               <  1:483699/138‑1 (MQ=12)                     tgtgtTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGc                              >  2:338023/1‑140 (MQ=17)                     tgtgtTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGc                              >  2:4618/1‑140 (MQ=17)                      gtgtTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCt                             >  1:814122/1‑140 (MQ=18)                       tgtTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTa                            >  2:656333/1‑140 (MQ=18)                       tgtTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTa                            >  2:324259/1‑140 (MQ=17)                       tgtTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTa                            >  1:268296/1‑140 (MQ=17)                        gtTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTaa                           <  1:1014456/140‑1 (MQ=17)                        gtTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTaa                           <  1:228181/140‑1 (MQ=17)                        gtTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTaa                           <  2:1243299/140‑1 (MQ=17)                           aGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATa                        >  1:222899/1‑140 (MQ=17)                           aGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGCCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATa                        >  1:1187458/1‑140 (MQ=12)                            gTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCCCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATAc                       <  1:766148/140‑1 (MQ=12)                             tGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACt                      >  2:1550972/1‑140 (MQ=17)                                 aGCGTCTGGAAAGGCGTGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGGCACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCctg                  >  1:1596362/1‑140 (MQ=12)                                 aGCGTCTGGAAAGGCGTGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTaa                           <  1:584713/131‑1 (MQ=17)                                 aGCGTCTGGAAAGGCGTGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTaa                           >  2:584713/1‑131 (MQ=17)                                 aGCGTCTGGAAAGGCGTGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCctg                  >  1:1577618/1‑140 (MQ=18)                                  gCGTCTGGAAAGGCGTGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCctga                 <  1:1012354/140‑1 (MQ=17)                                   cGTCTGGAAAGGCGTGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCctgac                >  1:1398828/1‑140 (MQ=17)                                      cTGGAAAGGCGTGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCctgactga             >  1:484223/1‑140 (MQ=18)                                      cTGGAAAGGCGTGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCctgactga             >  2:946769/1‑140 (MQ=17)                                       tGGAAAGGCGTGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCCtgactgac            >  1:1269605/1‑140 (MQ=18)                                       tGGAAAGGCGTGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCCtgactgac            <  2:655832/140‑1 (MQ=17)                                       tGGAAAGGCGTGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCCtgactgac            <  1:4618/140‑1 (MQ=17)                                       tGGAAAGGCGTGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCCtgactgac            <  2:533073/140‑1 (MQ=17)                                       tGGAAAGGCGTGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCCtgactgac            <  2:1318402/140‑1 (MQ=17)                                       tGGAAAGGCGTGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCCtgactgac            >  1:952273/1‑140 (MQ=17)                                       tGGAAAGGCGTGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCCtgactgac            <  2:1486590/140‑1 (MQ=17)                                       gggAAAGGCGTGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCCGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCCtgactgac            <  1:1082039/139‑1 (MQ=12)                                        ggAAAGGCGTGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGAcc           <  1:1414435/140‑1 (MQ=17)                                        ggAAAGGCGTGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGAcc           <  2:1501594/140‑1 (MQ=17)                                        ggAAAGGCGTGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGAcc           <  1:237662/140‑1 (MQ=17)                                           aaGGCGTGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGAt        <  2:222899/140‑1 (MQ=17)                                             ggCGTGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAg      <  1:1246538/140‑1 (MQ=17)                                              gCGTGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGt     >  1:728974/1‑140 (MQ=17)                                                 tGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGaa  >  2:1030943/1‑140 (MQ=11)                                                                                                                                           |                                                  CCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑GGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAA  >  NC_000913/4037786‑4037972

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: ‑