Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	273,607	A→G	pseudogene (761/1108 nt)	ykfC →	CP4‑6 prophage, conserved protein,Phage or Prophage Related

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	273,607	0	A	G	100.0%	154.7 / NA	44	pseudogene (761/1108 nt)	ykfC	CP4‑6 prophage, conserved protein,Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (24/20);  total (24/20)

GAAGGTGTGCCACAGGGCGGTGTTATATCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTACTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCGGTCAGAGAAAACTGGCAGTGGAAACCCGCGGTGGCGTACTGCCGCTATGCCGATGATTTTGTCCTCATCGTCAAAGGCACCAAAGCA  >  NC_000913/273489‑273731                                                                                                                       |                                                                                                                             gAAGGTGTGCCACAGGGCGGTGTTATATCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTGCTGGAATAACAGTATCCAAc                                                                                                          <  1:645215/139‑1 (MQ=255)     gtgtGCCACAGGGCGGTGTTATATCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTGCTGGAATAACAGTATCCAACgggg                                                                                                      >  1:1068236/1‑139 (MQ=255)                ggCGGTGTTATATCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTGCTGGAATAACAGTATCCAACGGGGCCGAAGTAcgg                                                                                           >  1:724245/1‑139 (MQ=255)                ggCGGTGTTATATCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTGCTGGAATAACAGTATCCAACGGGGCCGAAGTAcgg                                                                                           >  1:321629/1‑139 (MQ=255)                    gtgtTATATCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTGCTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCGGt                                                                                       >  2:362242/1‑139 (MQ=255)                               ccgcTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTGCTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCGGTCAGAGaaaa                                                                              <  1:780382/137‑1 (MQ=255)                               ccgcTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTGCTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCGGTCAGAGaaaa                                                                              >  2:780382/1‑137 (MQ=255)                               ccgcTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTGCTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCGGTCAGAGAAAACt                                                                            >  1:1010593/1‑139 (MQ=255)                               ccgcTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTGCTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCGGTCAGAGAAAACt                                                                            >  1:432142/1‑139 (MQ=255)                                  ctactaTCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTGCTGGAATGACAGTATCCAACGGGGCCGAAGTACGGCGGTCAGAGAAAACTGGc                                                                         >  2:728866/1‑139 (MQ=255)                                  cTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTGCTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCGGTCAGAGAAAACTGGc                                                                         >  1:431610/1‑139 (MQ=255)                                                catGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTGCTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCGGTCAGAGAAAACTGGCAGTGGAAACCcgcg                                                           <  2:946405/139‑1 (MQ=255)                                                 atGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTGCTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCGGTCAGAGAAAACTGGCAGTGGAAACCCGCgg                                                          <  2:321629/139‑1 (MQ=255)                                                  tGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTGCTGGAATAACAGTATCCAACGGGGc                                                                                                     <  1:640456/95‑1 (MQ=255)                                                  tGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTGCTGGAATAACAGTATCCAACGGGGc                                                                                                     >  2:640456/1‑95 (MQ=255)                                                  tGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTGCTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCGGTCAGAGAAAACTGGCAGTGGAAACCCGCGGt                                                         >  2:700672/1‑139 (MQ=255)                                                       aaTGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTGCTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCGGTCAGAGAAAACTGGCAGTGGAAACCCGCGGTGGCg                                                     <  2:850433/138‑1 (MQ=255)                                                       aaTGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTGCTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCGGTCAGAGAAAACTGGCAGTGGAAACCCGCGGTGGCg                                                     >  1:850433/1‑138 (MQ=255)                                                         tGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTGCTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCGGTCAGAGAAAACTGGCAGTGGAAACCCGCGGTGGCGTAc                                                  <  1:700672/139‑1 (MQ=255)                                                            gTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTGCTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCGGTCAGAGAAAACTGGCAGTGGAAACCCGCGGTGGCGTACTGc                                               <  1:362242/139‑1 (MQ=255)                                                                    aaTACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGTAAAGATCGGTGGTGCTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCGGTCAGAGAAAACTGGCAGTGGAAACCCGCGGTGGCGTACTGCCGCTATGc                                       <  2:724245/139‑1 (MQ=255)                                                                    aaTACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTGCTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCGGTCAGAGAAAACTGGCAGTGGAAACCCGCGGTGGCGTACTGCCGCTATGc                                       <  2:431610/139‑1 (MQ=255)                                                                           gCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTGCTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCGGTCAGAGAAAACTGGCAGTGGAAACCCGCGGTGGCGTACTGCCGCTATGCCgatgat                                >  1:691751/1‑139 (MQ=255)                                                                            cATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTGCTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCGGTCAGAGAAAACTGGCAGTGGAAACCCGCGGTGGCGTAc                                                  <  1:166544/120‑1 (MQ=255)                                                                            cATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTGCTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCGGTCAGAGAAAACTGGCAGTGGAAACCCGCGGTGGCGTAc                                                  >  2:166544/1‑120 (MQ=255)                                                                            cATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTGCTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCGGTCAGAGAAAACTGGCAGTGGAAACCCGCGGTGGCGTACTGCCGCTATGCCGATGAtt                               >  2:89619/1‑139 (MQ=255)                                                                            cATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTGCTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCGGTCAGAGAAAACTGGCAGTGGAAACCCGCGGTGGCGTACTGCCGCTATGCCGATGAtt                               >  2:401629/1‑139 (MQ=255)                                                                                aGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTGCTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCGGTCAGAGAAAACTGGCAGTGGAAACCCGCGGTGGCGTACTGCCGCt                                           >  2:304287/1‑123 (MQ=255)                                                                                aGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTGCTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCGGTCAGAGAAAACTGGCAGTGGAAACCCGCGGTGGCGTACTGCCGCt                                           <  1:304287/123‑1 (MQ=255)                                                                                aGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTGCTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCGGTCAGAGAAAACTGGCAGTGGAAACCCGCGGTGGCGTACTGCCGCTATGCCGATGATTTTGt                           >  2:84538/1‑139 (MQ=255)                                                                                aGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTGCTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCGGTCAGAGAAAACTGGCAGTGGAAACCCGCGGTGGCGTACTGCCGCTATGCCGATGATTTTGt                           >  2:1026988/1‑139 (MQ=255)                                                                                aGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTGCTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCGGTCAGAGAAAACTGGCAGTGGAAACCCGCGGTGGCGTACTGCCGCTATGCCGATGATTTTGt                           >  2:1098397/1‑139 (MQ=255)                                                                                 gcgcTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTGCTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCGGTCAGAGAAAACTGGCAGTGGAAAc                                                                >  2:266153/1‑101 (MQ=255)                                                                                 gcgcTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTGCTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCGGTCAGAGAAAACTGGCAGTGGAAAc                                                                <  1:266153/101‑1 (MQ=255)                                                                                 gcgcTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTGCTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCGGTCAGAGAAAACTGGCAGTGGAAACCCGCGGTGGCGTACTGCCGCTATGCCGATGATTTTGTc                          <  1:1098397/139‑1 (MQ=255)                                                                                 gcgcTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTGCTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCGGTCAGAGAAAACTGGCAGTGGAAACCCGCGGTGGCGTACTGCCGCTATGCCGATGATTTTGTc                          <  1:1026988/139‑1 (MQ=255)                                                                                  cgcTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTGCTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCGGTCAGAGAAAACTGGCAGTGGAAACCCGCGGTGGCGTACTGCCGCTAt                                         >  1:825826/1‑123 (MQ=255)                                                                                  cgcTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTGCTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCGGTCAGAGAAAACTGGCAGTGGAAACCCGCGGTGGCGTACTGCCGCTAt                                         <  2:825826/123‑1 (MQ=255)                                                                                      aCCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTGCTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCGGTCAGAGAAAACTGGCAGTGGAAACCCGCGGTGGCGTACTGCCGCTATGCCgatg                                  >  2:1115150/1‑126 (MQ=255)                                                                                      aCCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTGCTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCGGTCAGAGAAAACTGGCAGTGGAAACCCGCGGTGGCGTACTGCCGCTATGCCgatg                                  <  1:1115150/126‑1 (MQ=255)                                                                                          gAGCGGGAAAGCCAGAAAAGATCGGTGGTGCTGGAATAACa                                                                                                                   <  2:350744/41‑1 (MQ=255)                                                                                          gAGCGGGAAAGCCAGAAAAGATCGGTGGTGCTGGAATAACa                                                                                                                   >  1:350744/1‑41 (MQ=255)                                                                                                       aGAAAAGATCGGTGGTGCTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCGGTCAGAGAAAACTGGCAGTGGAAACCCGCGGTGGCGTACTGCCGCTATGCCGATGATTTTGTCCTCATCGTCAAAGGCACCAAAg    <  2:1010593/139‑1 (MQ=255)                                                                                                         aaaaGATCGGTGGTGCTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCGGTCAGAGAAAACTGGCAGTGGAAACCCGCGGTGGCGTACTGCCGCTATGCCGATGATTTTGTCCTCATCGTCAAAGGCACCAAAGca  <  2:1068236/139‑1 (MQ=255)                                                                                                                       |                                                                                                                             GAAGGTGTGCCACAGGGCGGTGTTATATCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTACTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCGGTCAGAGAAAACTGGCAGTGGAAACCCGCGGTGGCGTACTGCCGCTATGCCGATGATTTTGTCCTCATCGTCAAAGGCACCAAAGCA  >  NC_000913/273489‑273731

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