Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,944,122	T→C	noncoding (419/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,122	0	T	C	100.0%	97.7 / NA	38	noncoding (419/2904 nt)	rrlC	23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (22/16);  total (22/16)

GTATGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCA‑‑GTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGG  >  NC_000913/3943986‑3944198                                                                                                                                         |                                                                               gTATGTGTGTTAGTGGAAGCGTGTGGAAAGGCGCGCGGTACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCCCCa                                                                             <  1:353749/140‑1 (MQ=255) gTATGTGTGTTAGTGCAAGCGTTTGTAAACGCGCGCGCTACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCCCCa                                                                             <  1:218423/140‑1 (MQ=11)    tgtgtgTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCCCCAAgg                                                                          <  2:35206/140‑1 (MQ=21)    tgtgtgTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCCCCAAgg                                                                          <  1:93609/140‑1 (MQ=21)    tgtgtgTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCCCCAAgg                                                                          <  2:21490/140‑1 (MQ=21)    tgtgtgTTAGTGGAAGCCTCTGTAAAGGCGCCCGAAACATGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCCCCAAgg                                                                          <  2:192442/140‑1 (MQ=11)      tgtgtTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCCGGACACGTAGTATCCTTTCTGAATATGGGGGGACCATCCCCCAAGGCt                                                                        >  2:271014/1‑140 (MQ=255)        tgtTAGTGGAAGCGTCTGGAAATGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCCCCAAGGCTaa                                                                      <  1:290133/140‑1 (MQ=19)        tgtTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCCCCAAGGCTaa                                                                      <  1:533906/140‑1 (MQ=21)        tgtTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCCCCAAGGCTaa                                                                      <  1:74757/140‑1 (MQ=21)        tgtTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCCCCAAGGCTaa                                                                      >  1:382434/1‑140 (MQ=21)            aGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCAtccccc                                                                              >  1:622037/1‑128 (MQ=25)            aGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTCTCTGAATCTGGAGGGACCATCGCCCAAGGCTAAATAc                                                                  >  1:293465/1‑140 (MQ=11)            aGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGCAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCCCCAAGGCTAAATAc                                                                  >  1:660543/1‑140 (MQ=21)              tGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCCCCAAGGCTAAATACTc                                                                >  2:558194/1‑140 (MQ=21)                gAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCCCCAAGGCTAAATACTCct                                                              >  1:166247/1‑140 (MQ=21)                gAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCCCCAAGGCTAAATACTCct                                                              <  1:644792/140‑1 (MQ=21)                  aGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑GGGGGACCATCCCCCAAGGCTAAATACTCctgac                                                           >  1:171749/1‑140 (MQ=255)                       cTGGAAAGGCGCGCGCTACAGGCTGACAGCCCTGTACCCAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCCCCAAGGCTAAATACTCCtgactgac                                                       <  1:196931/140‑1 (MQ=11)                        tGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCCCCAAGGCTAAATACTCCTGACTGAcc                                                      <  1:663517/140‑1 (MQ=21)                             aGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCCCCAAGGCTAAATACTCCTGACTGACCGATAg                                                 >  1:1751/1‑140 (MQ=21)                             aGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCCCCAAGGCTAAATACTCCTGACTGACCGATAg                                                 >  1:31706/1‑140 (MQ=21)                               gcgcgcGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCCCCAAGGCTAAATACTCCTGACTGACCGATAGTg                                               <  2:427966/140‑1 (MQ=21)                                cgcgcgATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCCCCAAGGCTAAATACTCCGGGCGGACCGATAGTGa                                              >  1:502215/1‑140 (MQ=255)                                  cgcgATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCCGCAAGGCTAAATAGTCCTGACTGACCGATAGTGAAc                                            >  2:253388/1‑140 (MQ=255)                                            ggTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCCCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCA‑‑GTACCGTg                                >  2:183194/1‑140 (MQ=37)                                               gACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCCCCAAGGCTAAAAACTCCTGCCTGTCCCATAGTGTACCA‑‑GTACCGTGAgg                             >  1:310316/1‑140 (MQ=11)                                                 cAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCCCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCA‑‑GTACCGTGAGGGa                           >  1:393916/1‑140 (MQ=37)                                                 cAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCCCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCA‑‑GTACCGTGAGGGa                           >  2:579854/1‑140 (MQ=21)                                                 cAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑GGGGGACCATCCCCCAAGGCTAAATTCTCCTGACTGACCGATAGTGAACCA‑‑GTACCGTGAGGGaa                          >  2:276531/1‑140 (MQ=255)                                                   gCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCCCCAAGGCTAAATACTGCTGACTGACCGATAGTGAACCT‑‑GTACCGTGAGGGaaa                         >  1:225536/1‑140 (MQ=255)                                                    ccccGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCCCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCA‑‑GTACCGTGAGGGAAAg                        >  1:59400/1‑140 (MQ=37)                                                                   tGCACATATTGTGAGCTCGATGAGGAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCCCCAAGGCTAAATACTCCTGATTGACCGATGGTGAACCT‑‑GTACCGTGTGGGACAGGCGAAAAGAACCCCg         >  2:221715/1‑140 (MQ=11)                                                                        atatTGTGAGCTCGATGCGTAGTGCGGGACAGGTGCTATTCTGTCTGACTATGGGGGGACCATCCCCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCA‑‑GTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAg    <  2:264433/140‑1 (MQ=11)                                                                        atatTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCCCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCA‑‑GTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAg    <  2:560917/140‑1 (MQ=18)                                                                        atatTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCCCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCA‑‑GTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAg    >  2:90040/1‑140 (MQ=18)                                                                        atatTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCCCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCA‑‑GTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAg    >  1:560917/1‑140 (MQ=18)                                                                          atTGTGAGCTCGATGAGTAGGGTGGGACACGAGGTGTCCTGTCTGAATGTGGGGGGACCATCCCCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCA‑‑GTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAggg  <  2:293465/140‑1 (MQ=11)                                                                          atTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCCCCAAGACTAAATACTCCTGACTGACCGATAGTGAACCAGTGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAg    >  2:144447/1‑140 (MQ=11)                                                                                                                                         |                                                                               GTATGTGTGTTAGTGGAAGCGTCTGGAAAGGCGCGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCA‑‑GTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGG  >  NC_000913/3943986‑3944198

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