Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,760,492	T→C	intergenic (+98/‑55)	yfjI → / → yfjJ	CP4‑57 prophage; protein YfjI/CP4‑57 prophage; protein YfjJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,760,492	0	T	C	100.0%	149.2 / NA	40	intergenic (+98/‑55)	yfjI/yfjJ	CP4‑57 prophage; protein YfjI/CP4‑57 prophage; protein YfjJ
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (19/21);  total (19/21)

TAAAACAAAAATAATAAACATAACTCATTGTGTTTTTTGATTTTAATGACAATCTAAAAGGACTTAATTTAAGAGTAAAAACAAACTCTAAACCCAGTTCGGATAAAAGCTTGTCATCTTGATAGGTAATATATAGGTACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACG  >  NC_000913/2760355‑2760628                                                                                                                                          |                                                                                                                                         tAAAACAAAAATAATAAACATAACTCATTGTGTTTTTTGATTTTAATGACAATCTAAAAGGACTTAATTTAAGAGTAAAAACAAACTCTAAACCCAGTTCGGATAAAAGCTTGTCATCTTGATAGGTAATATATAGGCAc                                                                                                                                        >  2:328759/1‑140 (MQ=255)      cAAAAATAATAAACATAACTCATTGTGTTTTTTGATTTTAATGACAATCTAAAAGGACTTAATTTAAGAGTAAAAACAAACTCTAAACCCAGTTCGGATAAAAGCTTGTCATCTTGATAGGTAATATATAGGCACATAGa                                                                                                                                   <  2:1026455/140‑1 (MQ=255)           ataataAACATAACTCATTGTGTTTTTTGATTTTAATGACAATCTAAAAGGACTTAATTTAAGAGTAAAAACAAACTCTAAACCCAGTTCGGATAAAAGCTTGTCATCTTGATAGGTAATATATAGGCACATAGACACAt                                                                                                                              >  2:1281643/1‑140 (MQ=255)             aataaACATAACTCATTGTGTTTTTTGATTTTAATGACAATCTAAAAGGACTTAATTTAAGAGTAAAAACAAACTCTAAACCCAGTTCGGATAAAAGCTTGTCATCTTGATAGGTAATATATAGGCACATAGACACATCa                                                                                                                            <  2:1070107/140‑1 (MQ=255)               taaACATAACTCATTGTGTTTTTTGATTTTAATGACAATCTAAAAGGACTTAATTTAAGAGTAAAAACAAACTCTAAACCCAGTTCGGATAAAAGCTTGTCATCTTGATAGGTAATATATAGGCACATAGACACATCAGa                                                                                                                          >  1:1439659/1‑140 (MQ=255)                        cTCATTGTGTTTTTTGATTTTAATGACAATCTAAAAGGACTTAATTTAAGAGTAAAAACAAACTCTAAACCCAGTTCGGATAAAAGCTTGTCATCTTGATAGGTAATATATAGGCACATAGACACATCAGAACGATCACt                                                                                                                 <  2:938035/140‑1 (MQ=255)                           aTTGTGTTTTTTGATTTTAATGACAATCTAAAAGGACTTAATTTAAGAGTAAAAACAAACTCTAAACCCAGTTCGGATAAAAGCTTGTCATCTTGATAGGTAATATATAGGCACATAGACACATCAGAACGATCACTTAc                                                                                                              <  2:1012208/140‑1 (MQ=255)                             tgtgTTTTTTGATTTTAATGACAATCTAAAAGGACTTAATTTAAGAGTAAAAACAAACTCTAAACCCAGTTCGGATAAAAGCTTGTCATCTTGATAGGTAATATATAGGCACATAGACACATCAGAACGATCACt                                                                                                                 <  1:1278784/135‑1 (MQ=255)                             tgtgTTTTTTGATTTTAATGACAATCTAAAAGGACTTAATTTAAGAGTAAAAACAAACTCTAAACCCAGTTCGGATAAAAGCTTGTCATCTTGATAGGTAATATATAGGCACATAGACACATCAGAACGATCACt                                                                                                                 >  2:1278784/1‑135 (MQ=255)                             tgtgTTTTTTGATTTTAATGACAATCTAAAAGGACTTAATTTAAGAGTAAAAACAAACTCTAAACCCAGTTCGGATAAAAGCTTGTCATCTTGATAGGTAATATATAGGCACATAGACACATCAGAACGATCACTTACTc                                                                                                            >  2:3189/1‑140 (MQ=255)                             tgtgTTTTTTGATTTTAATGACAATCGAAAAGGACTCAATTTAAGAGTAAAAACAAACTCTAAACCCAATTCGAACAAAAACTTGTCATCTTGATAGAGAATATATAGGCACATAGACACATCAGAACGATCACTTACTc                                                                                                            >  2:1311794/1‑140 (MQ=255)                                       gATTTTAATGACAATCTAAAAGGACTTAATTTAAGAGTAAAAACAAACTCTAAACCCAGTTCGGATAAAAGCTTGTCATCTTGATAGGTAATATATAGGCACATAGACACATCAGAACGATCACTTACTCCAATATCAtt                                                                                                  <  2:768050/140‑1 (MQ=255)                                          tttAATGACAATCTAAAAGGACTTAATTTAAGAGTAAAAACAAACTCTAAACCCAGTTCGGATAAAAGCTTGTCATCTTGATAGGTAATATATAGGCACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGc                                                                                               <  2:1246905/140‑1 (MQ=255)                                            tAATGACAATCTAAAAGGACTTAATTTAAGAGTAAAAACAAACTCTAAACCCAGTTCGGATAAAAGCTTGTCATCTTGATAGGTAATATATAGGCACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCa                                                                                             <  1:1226289/140‑1 (MQ=255)                                                  cAATCTAAAAGGACTTAATTTAAGAGTAAAAACAAACTCTAAACCCAGTTCGGATAAAAGCTTGTCATCTTGATAGGTAATATATAGGCACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATAc                                                                                       <  1:817985/140‑1 (MQ=255)                                                      cTAAAAGGACTTAATTTAAGAGTAAAAACAAACTCTAAACCCAGTTCGGTTAAAAGCTTGTCATCTTGATAGGTAATATATAGGCACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACtata                                                                                   <  2:964343/140‑1 (MQ=255)                                                               cTTAATTTAAGAGTAAAAACAAACTCTAAACCCAGTTCGGATAAAAGCTTGTCATCTTGATAGGTAATATATAGGCACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACg                                                                          >  1:686451/1‑140 (MQ=255)                                                                ttAATTTAAGAGTAAAAACAAACTCTAAACCCAGTTCGGATAAAAGCTTGTCATCTTGATAGGTAATATATAGGCACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGa                                                                         >  1:574994/1‑140 (MQ=255)                                                                 tAATTTAAGAGTAAAAACAAACTCTAAACCCAGTTCGGATAAAAGCTTGTCATCTTGATAGGTAATATATAGGCACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGAt                                                                        >  2:422199/1‑140 (MQ=255)                                                                       aaGAGTAAAAACAAACTCTAAACCCAGTTCGGATAAAAGCTTGTCATCTTGATAGGTAATATATAGGCACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGaa                                                                  >  1:1134667/1‑140 (MQ=255)                                                                         gagTAAAAACAAACTCTAAACCCAGTTCGGATAAAAGCTTGTCATCTTGATAGGTAATATATAGGCACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAAtt                                                                <  1:328759/140‑1 (MQ=255)                                                                              aaaaCAAACTCTAAACCCAGTTCGGATAAAAGCTTGTCATCTTGATAGGTAATATATAGGCACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGAt                                                           >  1:78479/1‑140 (MQ=255)                                                                                    aacTCTAAACCCAGTTCGGATAAAAGCTTGTCATCTTGATAGGTAATATATAGGCACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATg                                                     <  1:1281643/140‑1 (MQ=255)                                                                                    aacTCTAAACCCAGTTCGGATAAAAGCTTGTCATCTTGATAGGTAATATATAGGCACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATg                                                     >  1:1291099/1‑140 (MQ=255)                                                                                        ctAAACCCAGTTCGGATAAAAGCTTGTCATCTTGATAGGTAATATATAGGCACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATc                                                 <  1:402879/140‑1 (MQ=255)                                                                                               cAGTTCGGATAAAAGCTTGTCATCTTGATAGGTAATATATAGGCACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTatat                                          <  2:1124240/140‑1 (MQ=255)                                                                                                 gTTCGGATAAAAGCTTGTCATCTTGATAGGTAATATATAGGCACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCg                                        >  2:1180071/1‑140 (MQ=255)                                                                                                 gTTCGGATAAAAGCTTGTCATCTTGATAGGTAATATATAGGCACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCg                                        >  2:1488053/1‑140 (MQ=255)                                                                                                 gTTCGGATAAAAGCTTGTCATCTTGATAGGTAATATATAGGCACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCg                                        >  2:301628/1‑140 (MQ=255)                                                                                                     ggATAAAAGCTTGTCATCTTGATAGGTAATATATAGGCACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGAt                                    >  1:68524/1‑140 (MQ=255)                                                                                                       aTAAAAGCTTGTCATCTTGATAGGTAATATATAGGCACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAg                                  >  1:1389923/1‑140 (MQ=255)                                                                                                        tAAAAGCTTGTCATCTTGATAGGTAATATATAGGCACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGa                                 >  2:1113714/1‑140 (MQ=255)                                                                                                              cTTGTCATCTTGATAGGTAATATATAGGCACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCaaa                           <  1:723548/140‑1 (MQ=255)                                                                                                              cTTGTCATCTTGATAGGTAATATATAGGCACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCaaa                           <  1:855722/140‑1 (MQ=255)                                                                                                              cTTGTCATCTTGATAGGTAATATATAGGCACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCaaa                           <  2:602473/140‑1 (MQ=255)                                                                                                              cTTGTCATCTTGATAGGTAATATATAGGCACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCaaa                           <  2:1134667/140‑1 (MQ=255)                                                                                                                  tCATCTTGATAGGTAATATATAGGCACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAgg                       <  2:574994/140‑1 (MQ=255)                                                                                                                          aTAGGTAATATATAGGCACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATtgat               <  1:3189/140‑1 (MQ=255)                                                                                                                            aGGTAATATATAGGCACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATtgatga             <  2:1291099/140‑1 (MQ=255)                                                                                                                                aaTATATAGGCACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTcc         >  1:1453947/1‑140 (MQ=255)                                                                                                                                       aGGCACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACg  >  2:53593/1‑140 (MQ=255)                                                                                                                                          |                                                                                                                                         TAAAACAAAAATAATAAACATAACTCATTGTGTTTTTTGATTTTAATGACAATCTAAAAGGACTTAATTTAAGAGTAAAAACAAACTCTAAACCCAGTTCGGATAAAAGCTTGTCATCTTGATAGGTAATATATAGGTACATAGACACATCAGAACGATCACTTACTCCAATATCATTGGCCAGGATACTATATGACATACGATAGTGAATTCGGATCACATGTATCCCTATATCGGGATAGAATCAAACAGGTTATTGATGACTCCCTAAACG  >  NC_000913/2760355‑2760628

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