Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,423,878	T→C	F33L (TTT→CTT)	trkG →	Rac prophage; potassium transporter subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,423,878	0	T	C	100.0%	149.1 / NA	39	F33L (TTT→CTT)	trkG	Rac prophage; potassium transporter subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (19/20);  total (19/20)

GTATAATCATAATATTGCAGCAAGGTGGTTATAATTGAAAGAATATTTAGATATGAATACATCTCATGTAAGAGTTGTTACTCATATGTGTGGGTTCCTGGTTTGGCTCTATAGTCTTTCAATGTTGCCACCAATGGTTGTAGCATTGTTTTATAAAGAAAAAAGCCTGTTTGTTTTCTTTATAACTTTCGTTATATTTTTTTGCATTGGTGGCGGAGCGTGGTATACAACTAAGAAATCTGGCATTCAATTACGTACCCGTGATGGGTTTATTATAATTGTAATGTTTTGG  >  NC_000913/1423730‑1424021                                                                                                                                                     |                                                                                                                                                gTATAATCATAATATTGCAGCAAGGTGGTTATAATTGAAAGAATATTTAGATATGAATACATCTCATGTAAGAGTTGTTACTCATATGTGTGGGTTCCTGGTTTGGCTCTATAGTCTTTCAATGTTGCCACCAATGGTTGTAGCATTGctt                                                                                                                                               <  1:479081/151‑1 (MQ=255)        cATAATATTGCAGCAAGGTGGTTATAATTGAAAGAATATTTAGATATGAATACATCTCATGTAAGAGTTGTTACTCATATGTGTGGGTTCCTGGTTTGGCTCTATAGTCTTTCAATGTTGCCACCAATGGTTGTAGCATTGCTTTATaaa                                                                                                                                         <  1:174246/150‑1 (MQ=255)                          tggtTATAATTGAAAGAATATTTAGATATGAATACATCTCATGTAAGAGTTGTTACTCATATGTGTGGGTTCCTGGTTTGGCTCTATAGTCTTTCAATGTTGCCACCAATGGTTGTAGCATTGCTTTATAAAGAAAAAAGCCTgtttgttt                                                                                                                      >  1:277573/1‑151 (MQ=255)                               ataATTGAAAGAATATTTAGATATGAATACATCTCATGTAAGAGTTGTTACTCATATGTGTGGGTTCCTGGTTTGGCTCTATAGTCTTTCAATGTTGCCACCAATGGTTGTAGCATTGCTTTATAAAGAAAAAAGCCTGTTTGTTTTCtt                                                                                                                  <  1:260802/150‑1 (MQ=255)                                         gAATATTTAGATATGAATACATCTCATGTAAGAGTTGTTACTCATATGTGTGGGTTCCTGGTTTGGCTCTATAGTCTTTCAATGTTGCCACCAATGGTTGTAGCATTGCTTTATAAAGAAAAAAGCCTGTTTGTTTTCTTTATAACTTTCg                                                                                                       >  1:280275/1‑151 (MQ=255)                                             atTTAGATATGAATACATCTCATGTAAGAGTTGTTACTCATATGTGTGGGTTCCTGGTTTGGCTCTATAGTCTTTCAATGTTGCCACCAATGGTTGTAGCATTGCTTTATAAAGAAAAAAGCCTGTTTGTTTTCTTTATAACTTTCGTtat                                                                                                   >  1:169261/1‑151 (MQ=255)                                             atTTAGATATGAATACATCTCATGTAAGAGTTGTTACTCATATGTGTGGGTTCCTGGTTTGGCTCTATAGTCTTTCAATGTTGCCACCAATGGTTGTAGCATTGCTTTATAAAGAAAAAAGCCTGTTTGTTTTCTTTATAACTTTCGTtat                                                                                                   <  1:479259/151‑1 (MQ=255)                                              tttAGATATGAATACATCTCATGTAAGAGTTGTTACTCATATGTGTGGGTTCCTGGTTTGGCTCTATAGTCTTTCAATGTTGCCACCAATGGTTGTAGCATTGCTTTATAAAGAAAAAAGCCTGTTTGTTTTCTTTATAACTTTCGTtat                                                                                                   >  1:556190/1‑150 (MQ=255)                                                tAGATATGAATACATCTCATGTAAGAGTTGTTACTCATATGTGTGGGTTCCTGGTTTGGCTCTATAGTCTTTCAATGTTGCCACCAATGGTTGTAGCATTGCTTTATAAAGAAAAAAGCCTGTTTGTTTTCTTTATAACTTTCGTTatat                                                                                                 <  1:289987/150‑1 (MQ=255)                                                 aGATATGAATACATCTCATGTAAGAGTTGTTACTCATATGTGTGGGTTCCTGGTTTGGCTCTATAGTCTTTCAATGTTGCCACCAATGGTTGTAGCATTGCTTTATAAAGAAAAAAGCCTGTTTGTTTTCTTTATAACTTTCGTTATAtt                                                                                                <  1:44403/150‑1 (MQ=255)                                                   atatGAATACATCTCATGTAAGAGTTGTTACTCATATGTGTGGGTTCCTGGTTTGGCTCTATAGTCTTTCAATGTTGCCACCAATGGTTGTAGCATTGCTTTATAAAGAAAAAAGCCTGTTTGTTTTCTTTATAACTTTCGTTATAttttt                                                                                             >  2:137610/1‑151 (MQ=255)                                                                tcATGTAAGAGTTGTTACTCATATGTGTGGGTTCCTGGTTTGGCTCTATAGTCTTTCAATGTTGCCACCAATGGTTGTAGCATTGCTTTATAAAGAAAAAAGCCTGTTTGTTTTCTTTATAACTTTCGTTATATTTTTTTGCATtggtgg                                                                                 >  1:38297/1‑150 (MQ=255)                                                                 cATGTAAGAGTTGTTACTCATATGTGTGGGTTCCTGGTTTGGCTCTATAGTCTTTCAATGTTGCCACCAATGGTTGTAGCATTGCTTTATAAAGAAAAAAGCCTGTTTGTTTTCTTTATAACTTTCGTTATATTTTTTTGCATTGGTGGCg                                                                               <  2:168178/151‑1 (MQ=255)                                                                 cATGTAAGAGTTGTTACTCATATGTGTGGGTTCCTGGTTTGGCTCTATAGTCTTTCAATGTTGCCACCAATGGTTGTAGCATTGCTTTATAAAGAAAAAAGCCTGTTTGTTTTCTTTATAACTTTCGTTATATTTTTTTGCATTGGTGGCg                                                                               >  2:462445/1‑151 (MQ=255)                                                                             gttACTCATATGTGTGGGTTCCTGGTTTGGCTCTATAGTCTTTCAATGTTGCCACCAATGGTTGTAGCATTGCTTTATAAAGAAAAAAGCCTGTTTGTTTTCTTTATAACTTTCGTTATATTTTTTTGCATTGGTGGCGGAGCGTGGtat                                                                    >  2:92680/1‑150 (MQ=255)                                                                              ttACTCATATGTGTGGGTTCCTGGTTTGGCTCTATAGTCTTTCAATGTTGCCACCAATGGTTGTAGCATTGCTTTATAAAGAAAAAAGCCTGTTTGTTTTCTTTATAACTTTCGTTATATTTTTTTGCAt                                                                                       <  2:299894/130‑1 (MQ=255)                                                                              ttACTCATATGTGTGGGTTCCTGGTTTGGCTCTATAGTCTTTCAATGTTGCCACCAATGGTTGTAGCATTGCTTTATAAAGAAAAAAGCCTGTTTGTTTTCTTTATAACTTTCGTTATATTTTTTTGCAt                                                                                       >  1:299901/1‑130 (MQ=255)                                                                              ttACTCATATGTGTGGGTTCCTGGTTTGGCTCTATAGTCTTTCAATGTTGCCACCAATGGTTGTAGCATTGCTTTATAAAGAAAAAAGCCTGTTTGTTTTCTTTATAACTTTCGTTATATTTTTTTGCATTGGTGGCGGAGCGTGGTATAc                                                                  >  1:173713/1‑151 (MQ=255)                                                                                  tCATATGTGTGGGTTCCTGGTTTGGCTCTATAGTCTTTCAATGTTGCCACCAATGGTTGTAGCATTGCTTTATAAAGAAAAAAGCCTGTTTGTTTTCTTTATAACTTTCGTTATATTTTTTTGCATTGGTGGCGGAGCGTGGTATACAACt                                                              >  1:585761/1‑151 (MQ=255)                                                                                   cATATGTGTGGGTTCCTGGTTTGGCTCTATAGTCTTTCAATGTTGCCACCAATGGTTGTAGCATTGCTTTATAAAGAAAAAAGCCTGTTTGTTTTCTTTATAACTTTCGTTATATTTTTTTGCATTGGTGGCGGAGCGTGGTATACAACt                                                              >  1:98604/1‑150 (MQ=255)                                                                                    atatGTGTGGGTTCCTGGTTTGGCTCTATAGTCTTTCAATGTTGCCACCAATGGTTGTAGCATTGCTTTATAAAGAAAAAAGCCTGTTTGTTTTCTTTATAACTTTCGTTATATTTTTTTGCATTGGTGGCGGAGCGTGGTATACAACTaa                                                            <  1:554234/151‑1 (MQ=255)                                                                                      atGTGTGGGTTCCTGGTTTGGCTCTATAGTCTTTCAATGTTGCCACCAATGGTTGTAGCATTGCTTTATAAAGAAAAAAGCCTGTTTGTTTTCTTTATAACTTTCGTTATATTTTTTTGCATTGGTGGCGGAGCGTGGTATACAACTAAg                                                           <  2:528790/150‑1 (MQ=255)                                                                                           tgGGTTCCTGGTTTGGCTCTATAGTCTTTCAATGTTGCCACCAATGGTTGTAGCATTGCTTTATAAAGAAAAAAGCCTGTTTGTTTTCTTTATAACTTTCGTTATATTTTTTTGCATTGGTGGCGGAGCGTGGTATACAACTAAGAAATCt                                                     >  2:523018/1‑151 (MQ=255)                                                                                                  cTGGTTTGGCTCTATAGTCTTTCAATGTTGCCACCAATGGTTGTAGCATTGCTTTATAAAGAAAAAAGCCTGTTTGTTTTCTTTATAACTTTCGTTATATTTTTTTGCATTGGTGGCGGAGCGTGGTATACAACTAAGAAATCTGGCATTc                                              >  2:124566/1‑151 (MQ=255)                                                                                                  cTGGTTTGGCTCTATAGTCTTTCAATGTTGCCACCAATGGTTGTAGCATTGCTTTATAAAGAAAAAAGCCTGTTTGTTTTCTTTATAACTTTCGTTATATTTTTTTGCATTGGTGGCGGAGCGTGGTATACAACTAAGAAATCTGGCATTc                                              >  2:170776/1‑151 (MQ=255)                                                                                                   tGGTTTGGCTCTATAGTCTTTCAATGTTGCCACCAATGGTTGTAGCATTGCTTTATAAAGAAAAAAGCCTGTTTGTTTTCTTTATAACTTTCGTTATATTTTTTTGCATTGGTGGCGGAGCGTGGTATACAACTAAGAAATCTGGCATTc                                              >  1:410534/1‑150 (MQ=255)                                                                                                        tGGCTCTATAGTCTTTCAATGTTGCCACCAATGGTTGTAGCATTGCTTTATAAAGAAAAAAGCCTGTTTGTTTTCTTTATAACTTTCGTTATATTTTTTTGCATTGGTGGCGGAGCGTGGTATACAACTAAGAAATCTGGCATTCAATTa                                         <  1:250914/150‑1 (MQ=255)                                                                                                        tGGCTCTATAGTCTTTCAATGTTGCCACCAATGGTTGTAGCATTGCTTTATAAAGAAAAAAGCCTGTTTGTTTTCTTTATAACTTTCGTTATATTTTTTTGCATTGGTGGCGGAGCGTGGTATACAACTAAGAAATCTGGCATTCAATTa                                         <  2:357686/150‑1 (MQ=255)                                                                                                         ggCTCTATAGTCTTTCAATGTTGCCACCAATGGTTGTAGCATTGCTTTATAAAGAAAAAGGCCTGTTTGTTTTCTTTATAACTTTCGTTATATTTTTTTGCATTGGTGGCGGAGCGTGGTATACAACTAAGAAATCTGGCATTCAATTACg                                       >  1:571207/1‑151 (MQ=255)                                                                                                             ctATAGTCTTTCAATGTTGCCACCAATGGTTGTAGCATTGCTTTATAAAGAAAAAAGCCTGTTTGTTTTCTTTATAACTTTCGTTATATTTTTTTGCATTGGTGGCGGAGCGTGGTATACAACTAAGAAATCTGGCATTCAATTACGTAcc                                   >  1:224198/1‑151 (MQ=255)                                                                                                                      ttCAATGTTGCCACCAATGGTTGTAGCATTGCTTTATAAAGAAAAAAGCCTGTTTGTTTTCTTTATAACTTTCGTTATATTTTTTTGCATTGGTGGCGGAGCGTGGTATACAACTAAGAAATCTGGCATTCAATTACGTACCCGTGATgg                           <  2:77801/150‑1 (MQ=255)                                                                                                                            gTTGCCACCAATGGTTGTAGCATTGCTTTATAAAGAAAAAAGCCTGTTTGTTTTCTTTATAACTTTCGTTATATTTTTTTGCATTGGTGGCGGAGCGTGGTATACAACTAAGAAATCTGGCATTCAATTACGTACCCGTGATGGGTttat                     <  2:10529/150‑1 (MQ=255)                                                                                                                               gCCACCAATGGTTGTAGCATTGCTTTATAAAGAAAAAAGCCTGTTTGTTTTCTTTATAACTTTCGTTATATTTTTTTGCATTGGTGGCGGAGCGTGGTATACAACTAAGAAATCTGGCATTCAATTACGTACCCGTGATGGGTTtattat                  >  1:189966/1‑150 (MQ=255)                                                                                                                                ccaccaATGGTTGTAGCATTGCTTTATAAAGAAAAAAGCCTGTTTGTTTTCTCTATAACTTTCGTTATATTTTTTTGCATTGGTGGCGGAGCGTGGTATACAACTAAGAAATCTGGCATTCAATTACGTACCCGTGATGGGTTTATTATaa                <  2:563570/151‑1 (MQ=255)                                                                                                                                  accaATGGTTGTAGCATTGCTTTATAAAGAAAAAAGCCTGTTTGTTTTCTTTATAACTTTCGTTATATTTTTTTGCATTGGTGGCGGAGCGTGGTATACAACTAAGAAATCTGGCATTCAATTACGTACCCGTGATGGGTTTATTATAAtt              <  1:602377/151‑1 (MQ=255)                                                                                                                                      aTGGTTGTAGCATTGCTTTATAAAGAAAAAAGCCTGTTTGTTTTCTTTATAACTTTCGTTATATTTTTTTGCATTGGTGGCGGAGCGTGGTATACAACTAAGAAATCTGGCATTCAATTACGTACCCGTGATGGGTTTATTATAATTGTaa          <  2:302979/151‑1 (MQ=255)                                                                                                                                            gTAGCATTGCTTTATAAAGAAAAAAGCCTGTTTGTTTTCTTTATAACTTTCGTTATATTTTTTTGCATTGGTGGCGGAGCGTGGTATACAACTAAGAAATCTGGCATTCAATTACGTACCCGTGATGGGTTTATTATAATTGTAATGtttt    <  2:3422/151‑1 (MQ=255)                                                                                                                                            gTAGCATTGCTTTATAAAGAAAAAAGCCTGTTTGTTTTCTTTATAACTTTCGTTATATTTTTTTGCATTGGTGGCGGAGCGTGGTATACAACTAAGAAATCTGGCATTCAATTACGTACCCGTGATGGGTTTATTATAATTGTAATGtttt    <  1:244596/151‑1 (MQ=255)                                                                                                                                              aGCATTGCTTTATAAAGAAAAAAGCCTGTTTGTTTTCTTTATAACTTTCGTTATATTTTTTTGCATTGGTGGCGGAGCGTGGTATACAACTAAGAAATCTGGCATTCAATTACGTACCCGTGATGGGTTTATTATAATTGTAATGTTTTgg  <  2:73150/151‑1 (MQ=255)                                                                                                                                              aGCATTGCTTTATAAAGAAAAAAGCCTGTTTGTTTTCTTTATAACTTTCGTTATATTTTTTTGCATTGGTGGCGGAGCGTGGTATACAACTAAGAAATCTGGCATTCAATTACGTACCCGTGATGGGTTTATTATAATTGTAATGTTTTg   <  2:358122/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                GTATAATCATAATATTGCAGCAAGGTGGTTATAATTGAAAGAATATTTAGATATGAATACATCTCATGTAAGAGTTGTTACTCATATGTGTGGGTTCCTGGTTTGGCTCTATAGTCTTTCAATGTTGCCACCAATGGTTGTAGCATTGTTTTATAAAGAAAAAAGCCTGTTTGTTTTCTTTATAACTTTCGTTATATTTTTTTGCATTGGTGGCGGAGCGTGGTATACAACTAAGAAATCTGGCATTCAATTACGTACCCGTGATGGGTTTATTATAATTGTAATGTTTTGG  >  NC_000913/1423730‑1424021

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