Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,769,648	(C)6→5	intergenic (+162/‑55)	yfjQ → / → yfjR	CP4‑57 prophage, uncharacterized protein/CP4‑57 prophage, putative DNA‑binding transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,769,643	0	C	.	100.0%	173.0 / NA	43	intergenic (+157/‑60)	yfjQ/yfjR	CP4‑57 prophage, uncharacterized protein/CP4‑57 prophage, putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (22/21);  total (22/21)

TTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCAGCCGTCTGGTTGCGGGT  >  NC_000913/2769508‑2769777                                                                                                                                        |                                                                                                                                       ttAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGcccc                                                                                                                                     <  2:256796/139‑1 (MQ=255)        cGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCG‑CCCCCTCTTTa                                                                                                                             >  2:1524101/1‑139 (MQ=255)           acacCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGccc                                                                                                                                      <  2:284595/128‑1 (MQ=255)           acacCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGccc                                                                                                                                      >  1:284595/1‑128 (MQ=255)                 ccTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCG‑CCCCCTCTTTAAAGTAATca                                                                                                                    >  2:1351464/1‑139 (MQ=255)                    gTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCG‑CCCCCTCTTTAAAGTAATCAcat                                                                                                                 <  2:699752/139‑1 (MQ=255)                     tCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCG‑CCCCCTCTTTAAAGTAATCAcatc                                                                                                                >  1:1255243/1‑139 (MQ=255)                            cGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCG‑CCCCCTCTTTa                                                                                                                             >  1:498799/1‑119 (MQ=255)                            cGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCG‑CCCCCTCTTTa                                                                                                                             <  2:498799/119‑1 (MQ=255)                            cGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCG‑CCCCCTCTTTAAAGTAATCACATCATTTtca                                                                                                         <  2:1510025/139‑1 (MQ=255)                                     cTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCG‑CCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTaa                                                                                                <  1:1584120/139‑1 (MQ=255)                                      tCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCG‑CCCCCTCTTTAAAGTAATCAcat                                                                                                                 <  1:456960/121‑1 (MQ=255)                                      tCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCG‑CCCCCTCTTTAAAGTAATCAcat                                                                                                                 >  2:456960/1‑121 (MQ=255)                                        gtgtGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGcccc                                                                                                                                     <  1:611224/100‑1 (MQ=255)                                        gtgtGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGcccc                                                                                                                                     >  2:611224/1‑100 (MQ=255)                                         tgtgCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCG‑CCCCCTCTTTAAAGTAATc                                                                                                                     >  1:843374/1‑114 (MQ=255)                                         tgtgCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCG‑CCCCCTCTTTAAAGTAATc                                                                                                                     <  2:843374/114‑1 (MQ=255)                                         tgtgCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCG‑CCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACttt                                                                                            >  2:395742/1‑139 (MQ=255)                                                cAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCG‑CCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGgag                                                                                     <  1:451147/139‑1 (MQ=255)                                                  aTCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGcccc                                                                                                                                     <  1:523118/90‑1 (MQ=255)                                                  aTCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGcccc                                                                                                                                     >  2:523118/1‑90 (MQ=255)                                                  aTCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCG‑CCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGaa                                                                                   <  2:1259412/139‑1 (MQ=255)                                                  aTCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCG‑CCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGaa                                                                                   <  2:1255243/139‑1 (MQ=255)                                                   tCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCG‑CCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAAt                                                                                  <  1:1524101/139‑1 (MQ=255)                                                    cccccGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCG‑CCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAAtc                                                                                 <  1:84732/139‑1 (MQ=255)                                                      cccGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGcccc                                                                                                                                     >  2:756982/1‑86 (MQ=255)                                                      cccGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGcccc                                                                                                                                     <  1:756982/86‑1 (MQ=255)                                                                   tttAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGccc                                                                                                                                      >  2:992695/1‑72 (MQ=255)                                                                   tttAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGccc                                                                                                                                      <  1:992695/72‑1 (MQ=255)                                                                    ttAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGcccc                                                                                                                                     <  1:964583/72‑1 (MQ=255)                                                                    ttAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGcccc                                                                                                                                     >  2:964583/1‑72 (MQ=255)                                                                    ttAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCG‑CCCCCTCTTTAAAGTAATCACATCATTTtcag                                                                                                        <  1:1116385/100‑1 (MQ=255)                                                                    ttAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCG‑CCCCCTCTTTAAAGTAATCACATCATTTtcag                                                                                                        >  2:1116385/1‑100 (MQ=255)                                                                    ttAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCG‑CCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCa                                                                 >  2:830419/1‑139 (MQ=255)                                                                    ttAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCG‑CCCCCTCTTTAAAGTAATAACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCa                                                                 >  1:1568087/1‑139 (MQ=255)                                                                     tAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCG‑CCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAg                                                                >  2:974448/1‑139 (MQ=255)                                                                            cGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCG‑CCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAAcgccg                                                         >  2:556712/1‑139 (MQ=255)                                                                                      ggCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGcccc                                                                                                                                     >  2:858330/1‑54 (MQ=255)                                                                                      ggCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGcccc                                                                                                                                     <  1:858330/54‑1 (MQ=255)                                                                                              cAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCG‑CCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTggctgt                                       <  1:1351464/139‑1 (MQ=255)                                                                                                 aaaTAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCG‑CCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTggctgtcag                                    >  2:634727/1‑139 (MQ=255)                                                                                                  aaTAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCG‑CCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTggctgtcagg                                   <  1:634727/139‑1 (MQ=255)                                                                                                            aGTTTCTGCGTTGTCCATACCCTGTCCG‑CCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGa                         >  2:807509/1‑139 (MQ=255)                                                                                                                   gCGTTGTCCATACCCTGTCCG‑CCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTGACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCAg                  <  1:974448/139‑1 (MQ=255)                                                                                                                        gTCCATACCCTGTCCG‑CCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAAtc                                                                                 >  1:680709/1‑71 (MQ=255)                                                                                                                        gTCCATACCCTGTCCG‑CCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAAtc                                                                                 <  2:680709/71‑1 (MQ=255)                                                                                                                                   gTCCG‑CCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGCGAATCTCTCATGACACAGGCAGACCGCCGCCATGACAGGCTGGCTGTCAGGCTGTTACTGATAATCAGCCGTCTGGTTGCGGGt  >  1:92525/1‑139 (MQ=14)                                                                                                                                   gTCCG‑CCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAAtc                                                                                 >  2:631132/1‑60 (MQ=255)                                                                                                                                   gTCCG‑CCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAAtc                                                                                 <  1:631132/60‑1 (MQ=255)                                                                                                                                   gTCCG‑CCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCggc                                              >  2:544927/1‑95 (MQ=255)                                                                                                                                   gTCCG‑CCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCggc                                              <  1:544927/95‑1 (MQ=255)                                                                                                                                   gTCCG‑CCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTggctgtcaggctgt                               >  2:797507/1‑110 (MQ=255)                                                                                                                                   gTCCG‑CCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTggctgtcaggctgt                               <  1:797507/110‑1 (MQ=255)                                                                                                                                   gTCCG‑CCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTggctgtca                                     <  1:333433/104‑1 (MQ=255)                                                                                                                                   gTCCG‑CCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTggctgtca                                     >  2:333433/1‑104 (MQ=255)                                                                                                                                   gTCCG‑CCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTggctgt                                       <  1:1592868/102‑1 (MQ=255)                                                                                                                                   gTCCG‑CCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTggctgt                                       >  2:1592868/1‑102 (MQ=255)                                                                                                                                   gTCCG‑CCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCAGCCGTCTGGTTGCGGGt  >  2:1429780/1‑139 (MQ=255)                                                                                                                                        |                                                                                                                                       TTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGTTTCTGCGTTGTCCATACCCTGTCCGCCCCCCTCTTTAAAGTAATCACATCATTTTCAGTCAGTTAACTTTCCTGGAGAATCTCTCATGACACAGGCAGAACGCCGCCATGACCGGCTGGCTGTCAGGCTGTCACTGATAATCAGCCGTCTGGTTGCGGGT  >  NC_000913/2769508‑2769777

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