Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,194,916	T→G	Q12H (CAA→CAC)	rluE ←	23S rRNA pseudouridine(2457) synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,194,916	0	T	G	100.0%	193.6 / NA	50	Q12H (CAA→CAC)	rluE	23S rRNA pseudouridine(2457) synthase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (27/23);  total (27/23)

TTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTTTGCATGGTATTTTCAGAGATTATGAATTGCCGCATTATAGCCTAATAACGCGCATCTTTCATGACGGCAAACAATAGGGTAGTATTGACAAGCCAATTACAAATCATTAACAAAAAATTGCTCTAAAGCATCCGTATCGCAGGACGC  >  NC_000913/1194768‑1195063                                                                                                                                                     |                                                                                                                                                    ttttGCGTCCGGCTTCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTtgtg                                                                                                                                                   >  1:298210/1‑151 (MQ=255)     gCGTCCGGCTTCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTGTGCAt                                                                                                                                                <  1:768881/150‑1 (MQ=255)     gCGTCCGGCTTCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTGTGCAt                                                                                                                                                <  1:707550/150‑1 (MQ=255)      cGTCCGGCTTCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTGTGCATgg                                                                                                                                              >  1:422211/1‑151 (MQ=255)      cGTCCGGCTTCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTGTGCATgg                                                                                                                                              >  1:772564/1‑151 (MQ=255)      cGTCCGGCTTCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTGTGCATgg                                                                                                                                              <  1:386629/151‑1 (MQ=255)          cGGCTTCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTGTGCATGGTAtt                                                                                                                                          <  1:46096/151‑1 (MQ=255)            gCTTCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTGTGCATGGTAtttt                                                                                                                                        <  1:864885/151‑1 (MQ=255)             cTTCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTGTGCATGGTATTTTc                                                                                                                                       >  2:392899/1‑151 (MQ=255)             cTTCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTGTGCATGGTATTTTc                                                                                                                                       >  1:397302/1‑151 (MQ=255)              ttCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTGTGCATGGTATTTTc                                                                                                                                       >  2:503520/1‑150 (MQ=255)               tCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTGTGCATGGTATTTTc                                                                                                                                       >  1:504027/1‑149 (MQ=255)                 aTCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTGTGCATGGTATTTTCag                                                                                                                                     >  2:548101/1‑149 (MQ=255)                   cGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTGTGCATGGTATTTTCag                                                                                                                                     >  2:80101/1‑147 (MQ=255)                       gAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTGTGCATGGTATTTTCAGAGATTATg                                                                                                                             <  2:720622/151‑1 (MQ=255)                         aCTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTGTGCATGGTATTTTCAGAGATTATGa                                                                                                                            >  2:143880/1‑150 (MQ=255)                            gCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTGTGCATGGTATTTTCAGAGATTATGAAtt                                                                                                                         <  1:826628/150‑1 (MQ=255)                            gCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTGTGCATGGTATTTTCAGAGATTATGAATTg                                                                                                                        <  2:891850/151‑1 (MQ=255)                              ggAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTGTGCATGGTATTTTCAGAGATTATGAATTGcc                                                                                                                      >  2:799736/1‑151 (MQ=255)                                  gaaCATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTGTGCATGGTATTTTCAGAGATTATGAATTGCCGCAt                                                                                                                  >  1:447280/1‑151 (MQ=255)                                    aCATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTGTGCATGGTATTTTCAGAGATTATGAATTGCCGCATta                                                                                                                <  1:32322/151‑1 (MQ=255)                                     cATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTGTGCATGGTATTTTCAGAGATTATGAATTGCCGCATtat                                                                                                               >  2:828344/1‑151 (MQ=255)                                             ggTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTGTGCATGGTATTTTCAGAGATTATGAATTGCCGCATTATAGCCtaat                                                                                                       >  2:402735/1‑151 (MQ=255)                                                    tGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTGTGCATGGTATTTTCAGAGATTATGAATTGCCGCATTATAGCCTAATAAcgcg                                                                                                 <  1:417867/150‑1 (MQ=255)                                                      aaCAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTGTGCATGGTATTTTCAGAGATTATGAATTGCCGCATTATAGCCTAATAACgcgctt                                                                                              <  1:258857/151‑3 (MQ=255)                                                      aaCAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTGTGCATGGTATTTTCAGAGATTATGAATTGCCGCATTATAGCCTAATAACGCGCAt                                                                                              >  1:26075/1‑151 (MQ=255)                                                         aGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTGTGCATGGTATTTTCAGAGATTATGAATTGCCGCATTATAGCCTAATAACGCGCATCtt                                                                                           >  1:770236/1‑151 (MQ=255)                                                           gATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTGTGCATGGTATTTTCAGAGATTATGAATTGCCGCATTATAGCCTAATAACGCGCATCTTTc                                                                                         >  2:219503/1‑151 (MQ=255)                                                              cacaCGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTGTGCATGGTATTTTCAGAGATTATGAATTGCCGCATTATAGCCTAATAACGCGCATCTTTCATg                                                                                      >  1:386717/1‑151 (MQ=255)                                                                    cgTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTTAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTGTGCATGGTATTTTCAGAGATTATGAATTGCCGCATTATAGCCTAATAACGCGCATCTTTCATGACGGCa                                                                                <  2:671414/151‑1 (MQ=255)                                                                      tGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTGTGCATGGTATTTTCAGAGATTATGAATTGCCGCATTATAGCCTAATAACGCGCATCTTTCATGACGGCaa                                                                               >  1:104841/1‑150 (MQ=255)                                                                           tGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTGTGCATGGTATTTTCAGAGATTATGAATTGCCGCATTATAGCCTAATAACGCGCATCTTTCATGACGGCAAACAAt                                                                          >  1:430035/1‑150 (MQ=255)                                                                                ttCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTGTGCATGGTATTTTCAGAGATTATGAATTGCCGCATTATAGCCTAATAACGCGCATCTTTCATGACGGCAAACAATAGGgt                                                                     <  1:566330/150‑1 (MQ=255)                                                                                 tCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTGTGCATGGTATTTTCAGAGATTATGAATTGCCGCATTATAGCCTAATAACGCGCATCTTTCATGACGGCAAACAATAGGgt                                                                     >  2:144943/1‑149 (MQ=255)                                                                                 tCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTGTGCATGGTATTTTCAGAGATTATGAATTGCCGCATTATAGCCTAATAACGCGCATCTTTCATGACGGCAAACAATAGGgt                                                                     <  1:479197/149‑1 (MQ=255)                                                                                      tttACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTGTGCATGGTATTTTCAGAGATTATGAATTGCCGCATTATAGCCTAATAACGCGCATCTTTCATGACGGCAAACAATAGGGTAGTATTg                                                              >  1:386082/1‑151 (MQ=255)                                                                                         aCGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTGTGCATGGTATTTTCAGAGATTATGAATTGCCGCATTATAGCCTAATAACGCGCATCTTTCATGACGGCAAACAATAGGGTAGTATTGAc                                                            <  1:809833/150‑1 (MQ=255)                                                                                                       gTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTGTGCATGGTATTTTCAGAGATTATGAATTGCCGCATTATAGCCTAATAACGCGCATCTTTCATGACGGCAAACAATAGGGTAGTATTGACAAGCCAATTACAAAt                                             >  1:235820/1‑151 (MQ=255)                                                                                                             aGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTGTGCATGGTATTTTCAGAGATTATGAATTGCCGCATTATAGCCTAATAACGCGCATCTTTCATGACGGCAAACAATAGGGTAGTATTGACAAGCCAATTACAAATCATTa                                        <  2:142792/150‑1 (MQ=255)                                                                                                                         tAACCTGGTGATTTCTAAAAGAAGTTTTGTGCATGGTATTTTCAGAGATTATGAATTGCCGCATTATAGCCTAATAACGCGCATCTTTCATGACGGCAAACAATAGGGTAGTATTGACAAGCCAATTACAAATCATTAACAAAAAATTg                             <  2:266639/149‑1 (MQ=255)                                                                                                                         tAACCTGGTGATTTCTAAAAGAAGTTTTGTGCATGGTATTTTCAGAGATTATGAATTGCCGCATTATAGCCTAATAACGCGCATCTTTCATGACGGCAAACAATAGGGTAGTATTGACAAGCCAATTACAAATCATTAACAAAAAATTg                             <  1:625287/149‑1 (MQ=255)                                                                                                                            ccTGGTGATTTCTAAAAGAAGTTTTGTGCATGGTATTTTCAGAGATTATGAATTGCCGCATTATAGCCTAATAACGCGCATCTTTCATGACGGCAAACAATAGGGTAGTATTGACAAGCCAATTACAAATCATTAACAAAAAATTGctct                         >  2:747738/1‑150 (MQ=255)                                                                                                                            ccTGGTGATTTCTAAAAGAAGTTTTGTGCATGGTATTTTCAGAGATTATGAATTGCCGCATTATAGCCTAATAACGCGCATCTTTCATGACGGCAAACAATAGGGTAGTATTGACAAGCCAATTACAAATCATTAACAAAAAATTGCTCTa                        <  1:828362/151‑1 (MQ=255)                                                                                                                               ggTGATTTCTAAAAGAAGTTTTGTGCATGGTATTTTCAGAGATTATGAATTGCCGCATTATAGCCTAATAACGCGCATCTTTCATGACGGCAAACAATAGGGTAGTATTGACAAGCCAATTACAAATCATTAACAAAAAATTGCTCTaaa                      <  2:59241/150‑1 (MQ=255)                                                                                                                               ggTGATTTCTAAAAGAAGTTTTGTGCATGGTATTTTCAGAGATTATGAATTGCCGCATTATAGCCTAATAACGCGCATCTTTCATGACGGCAAACAATAGGGTAGTATTGACAAGCCAATTACAAATCATTAACAAAAAATTGCTCTaa                       <  1:548116/149‑1 (MQ=255)                                                                                                                                     ttCTAAAAGAAGTTTTGTGCATGGTATTTTCAGAGATTATGAATTGCCGCATTATAGCCTAATAACGCGCATCTTTCATGACGGCAAACAATAGGGTAGTATTGACAAGCCAATTACAAATCATTAACAAAAAATTGCTCTAAAGCATCCg               >  2:618476/1‑151 (MQ=255)                                                                                                                                       cTAAAAGAAGTTTTGTGCATGGTATTTTCAGAGATTATGAATTGCCGCATTATAGCCTAATAACGCGCATCTTTCATGACGGCAAACAATAGGGTAGTATTGACAAGCCAATTACAAATCATTAACAAAAAATTGCTCTAAAGCATCCGt              >  1:261669/1‑150 (MQ=255)                                                                                                                                        tAAAAGAAGTTTTGTGCATGGTATTTTCAGAGATTATGAATTGCCGCATTATAGCCTAATAACGCGCATCTTTCATGACGGCAAACAATAGGGTAGTATTGACAAGCCAATTACAAATCATTAACAAAAAATTGCTCTAAAGCATCCGTAt            >  2:164813/1‑151 (MQ=255)                                                                                                                                         aaaaGAAGTTTTGTGCATGGTATTTTCAGAGATTATGAATTGCCGCATTATAGCCTAATAACGCGCATCTTTCATGACGGCAAACAATAGGGTAGTATTGACAAGCCAATTACAAATCATTAACAAAAAATTGCTCTAAAGCATCCGt              <  1:520129/148‑1 (MQ=255)                                                                                                                                          aaaGAAGTTTTGTGCATGGTATTTTCAGAGATTATGAATTGCCGCATTATAGCCTAATAACGCGCATCTTTCATGACGGCAAACAATAGGGTAGTATTGACAAGCCAATTACAAATCATTAACAAAAAATTGCTCTAAAGCATCCGTATCg          >  2:264202/1‑151 (MQ=255)                                                                                                                                               agTTTTGTGCATGGTATTTTCAGAGATTATGGATTGCCGCATTATAGCCTAATAACGCGCATCTTTCATGACGGCAAACAATAGGGTAGTATTGACAAGCCAATTACAAATCATTAACAAAAAATTGCTCTAAAGCATCCGTATCGCAGGa     <  1:393882/151‑1 (MQ=255)                                                                                                                                                  tttgtGCATGGTATTTTCAGAGATTATGAATTGCCGCATTATAGCCTAATAACGCGCATCTTTCATGACGGCAAACAATAGGGTAGTATTGACAAGCCAATTACAAATCATTAACAAAAAATTGCTCTAAAGCATCCGTATCGCAGGACGc  >  1:764047/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                    TTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTTTGCATGGTATTTTCAGAGATTATGAATTGCCGCATTATAGCCTAATAACGCGCATCTTTCATGACGGCAAACAATAGGGTAGTATTGACAAGCCAATTACAAATCATTAACAAAAAATTGCTCTAAAGCATCCGTATCGCAGGACGC  >  NC_000913/1194768‑1195063

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