Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,632,654	T→C	intergenic (‑369/‑236)	ydfJ ← / → ynfT	putative transporter YdfJ/Qin prophage; protein YnfT

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,632,654	0	T	C	100.0%	125.8 / NA	47	intergenic (‑369/‑236)	ydfJ/ynfT	putative transporter YdfJ/Qin prophage; protein YnfT
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (27/20);  total (27/20)

TTTACCAGCCAATCGCAGCACGTTCTTGCATAAGGTGTGTCTGCGGTTTTTCAACTATTCAGATACATCACTCCCATCACATTCATTCCTCCGCATCAAAGGCATATAGGCTATATCACCTTGATATTTTTCTTTTTCAGATAAAAACTGTTATCTATGTATACTTTTAAACCCAATCCGTGTAGAGTCTCTACATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAGCATTCCTG  >  NC_000913/1632520‑1632790                                                                                                                                       |                                                                                                                                         tttACCAGCCAATCGCAGCACGTTCTTGCATAAGGTGTGTCTGCGGTTTTTCAACTATTCAGATACATCACTCCCATCACATTCATTCCTCCGCATCAAAGGCATATAGGCTATATCACCTTGATATTTTTCTTCTTCAg                                                                                                                                     >  2:1125636/1‑140 (MQ=255)   tACCAGCCAATCGCAGCACGTTCTTGCATAAGGTGTGTCTGCGGTTTTTCAACTATTCAGATACATCACTCCCATCACATTCATTCCTCCGCATCAAAGGCATATAGGCTATATCACCTTGATATTTTTCTTCTTCAGAt                                                                                                                                   >  1:181795/1‑140 (MQ=255)        gCCAATCGCAGCACGTTCTTGCATAAGGTGTGTCTGCGGTTTTTCAACTATTCAGATACATCACTCCCATCACATTCATTCCTCCGCATCAAAGGCATATAGGCTATATCACCTTGATATTTTTCTTCTTCAGATaaaaa                                                                                                                              <  2:1214695/140‑1 (MQ=255)          cAATCGCAGCACGTTCTTGCATAAGGTGTGTCTGCGGTTTTTCAACTATTCAGATACATCACTCCCATCACATTCATTCCTCCGCATCAAAGGCATATAGGCTATATCACCTTGATATTTTTCTTCTTCAGATAAAAACt                                                                                                                            <  2:807892/140‑1 (MQ=255)          cAATCGCAGCACGTTCTTGCATAAGGTGTGTCTGCGGTTTTTCAACTATTCAGATACATCACTCCCATCACATTCATTCCTCCGCATCAAAGGCATATAGGCTATATCACCTTGATATTTTTCTTCTTCAGATAAAAACt                                                                                                                            <  2:740808/140‑1 (MQ=255)           aaTCGCAGCACGTTCTTGCATAAGGTGTGTCTGCGGTTTTTCAACTATTCAGATACATCACTCCCATCACATTCATTCCTCCGCATCAAAGGCATATAGGCTATATCACCTTGATATTTTTCTTCTTCAGATAAAAACTg                                                                                                                           >  2:737759/1‑140 (MQ=255)              cGCAGCACGTTCTTGCATAAGGTGTGTCTGCGGTTTTTCAACTATTCAGATACATCACTCCCATCACATTCATTCCTCCGCATCAAAGGCATATAGGCTATATCACCTTGATATTTTTCTTCTTCAGATAAAAACTGTTa                                                                                                                        <  2:737086/140‑1 (MQ=37)              cGCAGCACGTTCTTGCATAAGGTGTGTCTGCGGTTTTTCAACTATTCAGATACATCACTCCCATCACATTCATTCCTCCGCATCAAAGGCATATAGGCTATATCACCTTGATATTTTTCTTCTTCAGATAAAAACTGTTa                                                                                                                        <  1:1251656/140‑1 (MQ=37)                 agcaCGTTCTTGCATAAGGTGTGTCTGCGGTTTTTCAACTATTCAGATACATCACTCCCATCACATTCATTCCTCCGCATCAAAGGCATATAGGCTATATCACCTTGATATTTTTCTTCTTCAGATAAAAACTGTTATCt                                                                                                                     <  1:10761/140‑1 (MQ=37)                      gTTCTTGCATAAGGTGTGTCTGCGGTTTTTCAACTATTCAGATACATCACTCCCATCACATTCATTCCTCCGCATCAAAGGCATATAGGCTATATCACCTTGATATTTTTCTTCTTCAGATAAAAACTGTTATCTATGta                                                                                                                <  1:737759/140‑1 (MQ=37)                       ttCTTGCATAAGGTGTGTCTGCGGTTTTTCAACTATTCAGATACATCACTCCCATCACATTCATTCCTCCGCATCAAAGGCATATAGGCTATATCACCTTGATATTTTTCTTCTTCAGATAAAAACTGTTATCTATGtat                                                                                                               <  1:1127632/140‑1 (MQ=37)                        tCTTGCATAAGGTGTGTCTGCGGTTTTTCAACTATTCAGATACATCACTCCCATCACATTCATTCCTCCGCATCAAAGGCATATAGGCTATATCACCTTGATATTTTTCTTCTTCAGATAAAAACTGTTATCTATGtata                                                                                                              >  1:288641/1‑140 (MQ=255)                                 aGGTGTGTCTGCGGTTTTTCAACTATTCAGATACATCACTCCCATCACATTCATTCCTCCGCATCAAAGGCACATAGGCTATATCACCTTGATATTTTTCTTCTTCAGATAAAAACTGTTATCTATGTATACTTTTAAAc                                                                                                     <  1:110529/140‑1 (MQ=12)                                  ggTGTGTCTGCGGTTTTTCAACTATTCAGATACATCACTCCCATCACATTCATTCCTCCGCATCAAAGGCATATAGGCTATATCACCTTGATATTTTTCTTCTTCAGATAAAAACTGTTATCTATGTATACTTTTAAAcc                                                                                                    >  2:517074/1‑140 (MQ=17)                                  ggTGTGTCTGCGGTTTTTCAACTATTCAGATACATCACTCCCATCACATTCATTCCTCCGCATCAAAGGCATATAGGCTATATCACCTTGATATTTTTCTTCTTCAGATAAAAACTGTTATCTATGTATACTTTTAAAcc                                                                                                    >  2:187543/1‑140 (MQ=18)                                  ggTGTGTCTGCGGTTTTTCAACTATTCAGATACATCACTCCCATCACATTCATTCCTCCGCATCAAAGGCATATAGGCTATATCACCTTGATATTTTTCTTCTTCAGATAAAAACTGTTATCTATGTATACTTTTAAAcc                                                                                                    >  1:671312/1‑140 (MQ=18)                                  ggTGTGTCTGCGGTTTTTCAACTATTCAGATACATCACTCCCATCACATTCATTCCTCCGCATCAAAGGCATATAGGCTATATCACCTTGATATTTTTCTTCTTCAGATAAAAACTGTTATCTATGTATACTTTTAAAcc                                                                                                    >  2:1189725/1‑140 (MQ=18)                                                                                 attcattcCTCCGCATCAAAGGCATATAGGCTATATCACCTTGATATTTTTCTTCTTCAGATAAAAACTGTTATCTATGTATACTTTTAAACCCAATCCGTGTAGAGTCTCTACATAAGATAGTTTGCAGTTGCCGCTTc                                                     >  1:10889/1‑140 (MQ=17)                                                                                  ttcattcCTCCGCATCAAAGGCATATAGGCTATATCACCTTGATATTTTTCTTCTTCAGATAAAAACTGTTATCTATGTATACTTTTAAACCCAATCCGTGTAGAGTCTCTACATAAGATAGTTTGCAGTTGCCGCTTCa                                                    <  2:1232196/140‑1 (MQ=18)                                                                                     attcCTCCGCATCAAAGGCATATAGGCTATATCACCTTGATATTTTTCTTCTTCAGATAAAAACTGTTATCTATGTATACTTTTAAACCCAATCCGTGTAGAGTCTCTACATAAGATAGTTTGCAGTTGCCGCTTCAGCt                                                 >  2:930730/1‑140 (MQ=14)                                                                                      ttcCTCCGCATCAAAGGTATATAGGCTATATCACCTTGATATTTTTCTTCTTCAGATAAAAACTGTTATCTATGTATACTTTTAAACCCAATCCGTGTAGAGTCTCTACATAAGATAGTTTGCAGTTGCCGCTTCAGCtt                                                <  1:1186673/140‑1 (MQ=18)                                                                                      ttcCTCCGCATCAAAGGCATATAGGCTATATCACCTTGATATTTTTCTTCTTCAGATAAAAACTGTTATCTATGTATACTTTTAAACCCAATCCGTGTAGAGTCTCTACATAAGATAGTTTGCAGTTGCCGCTTCAGCtt                                                <  1:1125636/140‑1 (MQ=18)                                                                                      ttcCTCCGCATCAAAGGCATATAGGCTATATCACCTTGATATTTTTCTTCTTCAGATAAAAACTGTTATCTATGTATACTTTTAAACCCAATCCGTGTAGAGTCTCTACATAAGATAGTTTGCAGTTGCCGCTTCAGCtt                                                <  1:295047/140‑1 (MQ=17)                                                                                         ctccGCATCAAAGGCATATAGGCTATATCACCTTGATATTTTTCTTCTTCAGATAAAAACTGTTATCTATGTATACTTTTAAACCCAATCCGTGTAGAGTCTCTACATAAGATAGTTTGCAGTTGCCGCTTCAGCTTgcg                                             >  2:1200693/1‑140 (MQ=14)                                                                                          tccGCATCAAAGGCATATAGGCTATATCACCTTGATATTTTTCTTCTTCAGATAAAAACTGTTATCTATGTATACTTTTAAACCCAATCCGTGTAGAGTCTCTACATAAGATAGTTTGCAGTTGCCGCTTCAGCTTgcgc                                            <  1:233195/140‑1 (MQ=18)                                                                                             gCATCAAAGGCATATAGGCTATATCACCTTGATATTTTTCTTCTTCAGATAAAAACTGTTATCTATGTATACTTTTAAACCCAATCCGTGTAGAGTCTCTACATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCAt                                         >  1:960766/1‑140 (MQ=14)                                                                                              cATCAAAGGCATATAGGCTATATCACCTTGATATTTTTCTTCTTCAGATAAAAACTGTTATCTATGTATACTTTTAAACCCAATCCGTGTAGAGTCTCTACATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATa                                        >  2:102834/1‑140 (MQ=14)                                                                                                                tataTCACCTTGATATTTTTCTTCTTCAGATAAAAACTGTTATCTATGTATACTTTTAAACCCAATCCGTGTAGAGTCTCTACATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTtgct                      >  2:1258567/1‑140 (MQ=17)                                                                                                                 atatCACCTTGATATTTTTCTTCTTCAGATAAAAACTGTTATCTATGTATACTTTTAAACCCAATCCGTGTAGAGTCTCTACATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTtgctg                     >  2:114786/1‑140 (MQ=18)                                                                                                                 atatCACCTTGATATTTTTCTTCTTCAGATAAAAACTGTTATCTATGTATACTTTTAAACCCAATCCGTGTAGAGTCTCTACATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTtgctg                     >  1:1129042/1‑140 (MQ=18)                                                                                                                  tatCACCTTGATATTTTTCTTCTTCAGATAAAAACTGTTATCTATGTATACTTTTAAACCCAATCCGTGTAGAGTCTCTACATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTtgctgc                    <  1:1200693/140‑1 (MQ=21)                                                                                                                  tatCACCTTGATATTTTTCTTCTTCAGATAAAAACTGTTATCTATGTATACTTTTAAACCCAATCCGTGTAGAGTCTCTACATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTtgctgc                    >  2:32982/1‑140 (MQ=18)                                                                                                                   atCACCTTGATATTTTTCTTCTTCAGATAAAAACTGTTATCTATGTATACTTTTAAACCCAATCCGTGTAGAGTCTCTACATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGcc                   <  2:181795/140‑1 (MQ=21)                                                                                                                   atCACCTTGATATTTTTCTTCTTCAGATAAAAACTGTTATCTATGTATACTTTTAAACCCAATCCGTGTAGAGTCTCTACATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGcc                   >  2:582525/1‑140 (MQ=21)                                                                                                                       ccTTGATATTTTTCTTCTTCAGATAAAAACTGTTATCTATGTATACTTTTAAACCCAATCCGTGTAGAGTCTCTACATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCt               <  1:163385/140‑1 (MQ=25)                                                                                                                       ccTTGATATTTTTCTTCTTCAGATAAAAACTGTTATCTATGTATACTTTTAAACCCAATCCGTGTAGAGTCTCTACATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCt               >  2:1029015/1‑140 (MQ=18)                                                                                                                        cTTGATATTTTTCTTCTTCAGATAAAAACTGTTATCTATGTATACTTTTAAACCCAATCCGTGTAGAGTCTCTACATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTg              >  1:697575/1‑140 (MQ=18)                                                                                                                        cTTGATATTTTTCTTCTTCAGATAAAAACTGTTATCTATGTATACTTTTAAACCCAATCCGTGTAGAGTCTCTACATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTg              >  2:1236581/1‑140 (MQ=21)                                                                                                                        cTTGATATTTTTCTTCTTCAGATAAAAACTGTTATCTATGTATACTTTTAAACCCAATCCGTGTAGAGTCTCTACATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTg              >  2:289580/1‑140 (MQ=18)                                                                                                                        cTTGATATTTTTCTTCTTCAGATAAAAACTGTTATCTATGTATACTTTTAAACCCAATCCGTGTAGAGTCTCTACATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTg              >  2:361193/1‑140 (MQ=21)                                                                                                                        cTTGATATTTTTCTTCTTCAGATAAAAACTGTTATCTATGTATACTTTTAAACCCAATCCGTGTAGAGTCTCTACATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTg              >  2:420136/1‑140 (MQ=18)                                                                                                                        cTTGATATTTTTCTTCTTCAGATAAAAACTGTTATCTATGTATACTTTTAAACCCAATCCGTGTAGAGTCTCTACATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTACTGCCACCTg              >  2:789645/1‑140 (MQ=16)                                                                                                                            atatTTTTCTTCTTCAGATAAAAACTGTTATCTATGTATACTTTTAAACCCAATCCGTGTAGAGTCTCTACATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAg          >  2:98413/1‑140 (MQ=18)                                                                                                                            atatTTTTCTTCTTCAGATAAAAACTGTTATCTATGTATACTTTTAAACCCAATCCGTGTAGAGTCTCTACATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAg          >  2:244244/1‑140 (MQ=18)                                                                                                                            atatTTTTCTTCTTCAGATAAAAACTGTTATCTATGTATACTTTTAAACCCAATCCGTGTAGAGTCTCTACATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAg          <  1:789645/140‑1 (MQ=25)                                                                                                                              atTTTTCTTCTTCAGATAAAAACTGTTATCTATGTATACTTTTAAACCCAATCCGTGTAGAGTCTCTACATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAGCa        <  1:1007960/140‑1 (MQ=25)                                                                                                                              atTTTTCTTCTTCAGATAAAAACTGTTATCTATGTATACTTTTAAACCCAATCCGTGTAGAGTCTCTACATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAGCa        <  2:1129042/140‑1 (MQ=21)                                                                                                                                 tttcttcttcAGATAAAAACTGTTATCTATGTATACTTTTAAACCCAATCCGTGTAGAGTCTCTACATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAGCATTc     >  1:1130670/1‑140 (MQ=25)                                                                                                                                    cttcttcAGATAAAAACTGTTATCTATGTATACTTTTAAACCCAATCCGTGTAGAGTCTCTACATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAGCATTCCTg  <  1:305528/140‑1 (MQ=25)                                                                                                                                       |                                                                                                                                         TTTACCAGCCAATCGCAGCACGTTCTTGCATAAGGTGTGTCTGCGGTTTTTCAACTATTCAGATACATCACTCCCATCACATTCATTCCTCCGCATCAAAGGCATATAGGCTATATCACCTTGATATTTTTCTTTTTCAGATAAAAACTGTTATCTATGTATACTTTTAAACCCAATCCGTGTAGAGTCTCTACATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAGCATTCCTG  >  NC_000913/1632520‑1632790

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