Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,068,622	A→G	intergenic (‑1301/‑13)	insH‑6 ← / → yoeA	CP4‑44 prophage; IS5 transposase and trans‑activator/CP4‑44 prophage; TonB‑dependent receptor plug domain‑containing protein YoeA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,068,622	0	A	G	100.0%	121.8 / NA	34	intergenic (‑1301/‑13)	insH‑6/yoeA	CP4‑44 prophage; IS5 transposase and trans‑activator/CP4‑44 prophage; TonB‑dependent receptor plug domain‑containing protein YoeA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (19/15);  total (19/15)

AAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATGAAGATGACAATAATTATAGATATATGAGAACATAAATGAAAATAATTATCATTACAGTAATCATTTGTACTTTGTATTAATGAGGGATGAAATGTTATATAATATACCTTGTCGAATTTATATCCTTTCCACTCTGTCATTATGCATTTCTGGGATAGTTTCTACTGCAACCGCAACTTCTTCAGAAACAAAAATCAGCAACGAAGAGACGC  >  NC_000913/2068486‑2068755                                                                                                                                         |                                                                                                                                      aaaaTATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATGAAGATGACAATAATTATAGATATATGAGAACATAAATGAAAATAATTATCATTACAGTAATCATTTGTACTTTGTATTGATg                                                                                                                                    >  2:1295297/1‑140 (MQ=255)        ctatcAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATGAAGATGACAATAATTATAGATATATGAGAACATAAATGAAAATAATTATCATTACAGTAATCATTTGCACTTTGTATTGATGAGGGATg                                                                                                                             <  2:1161835/140‑1 (MQ=255)          atcAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATGAAGATGACAATAATTATAGATATATGAGAACATAAATGAAAATAATTATCATTACAGTAATCATTTGTACTTTGTATTGATGAGGGATGaa                                                                                                                           <  2:943198/140‑1 (MQ=255)          atcAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATGAAGATGACAATAATTATAGATATATGAGAACATAAATGAAAATAATTATCATTACAGTAATCATTTGTACTTTGTATTGATGAGGGATGaa                                                                                                                           <  1:1904790/140‑1 (MQ=255)           tcAAAAAACATGATTCAGGTCAAAACTGTCTGATATTTACAATTTATGAAGATGACAATAATTATAGATATATGAGAACATAAATGAAAATAATTATCATTACAGTAATCATTTGTACTTTGTATTGATGAGGGATGaaa                                                                                                                          <  1:795362/140‑1 (MQ=255)               aaaaCATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATGAAGATGACAATAATTATAGATATATGAGAACATAAATGAAAATAATTATCATTACAGTAATCATTTGTACTTTGTATTGATGAGGGATGAAATGtt                                                                                                                      <  1:1685421/140‑1 (MQ=255)                aaaCATGATTAAGGTGAAAAATGTTTGATGTTTACAACTTATGAAGATGACAATAATTATAGATATATGAGAACATAAATGAAAATAATTATCATTACAGTAATCATTTGTACTTTGTATTGATGAGGGATGAAATGTta                                                                                                                     <  1:899489/140‑1 (MQ=255)                aaaCATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATGAAGATGACAATAATTATAGATATATGAGAACATAAATGAAAATAATTATCATTACAGTAATCATTTGTACTTTGTATTGATGAGGGATGAAATGTta                                                                                                                     <  1:80649/140‑1 (MQ=255)                         tAAGGTCAAAAATGTTTGATATTTACAATTTATGAAGATGACAATAATTATAGATATATGAGAACATAAATGAAAATAATTATCATTACAGTAATCATTTGTACTTTGTATTGATGAGGGATGAAATGTTatataatata                                                                                                            <  2:1216101/140‑1 (MQ=255)                         tAAGGTCAAAAATGTTTGATATTTACAATTTATGAAGATGACAATAATTATAGATATATGAGAACATAAATGAAAATAATTATCATTACAGCAATCATTTCTACTTGGCATTGATGATGGATCAAATGTTatataatata                                                                                                            >  2:899489/1‑140 (MQ=255)                             gTCAAAAATGTTTGATATTTACAATTTATGAAGATGACAATAATTATAGATATATGAGAACATAAATGAAAATAATTATCATTACAGTAATCATTTGTACTTTGTATTGATGAGGGATGAAATGTTATATAATATACCtt                                                                                                        >  2:989262/1‑140 (MQ=255)                                 aaaaTGTTTGATATTCACAATTGATCAAGTTGACAATGATTATAGATATATGAGAACATAAATGAAAATAATTATCATTACAGTAATCATTTGTACTTTGTATTGATGAGGGATGAAATGTTATATAATATACCTTGTCg                                                                                                    <  1:896558/140‑1 (MQ=255)                                               ttACAATTTATGAAGCTGACAACAATTATCGATATATCAGAACATAAATGAAAATAATTATCATTACAGTAATCATTTGTACTTTGTATTGATGAGGGATGAAATGTTATATAATATACCTTGTCGAATTTATATCCttt                                                                                      <  1:497022/140‑1 (MQ=255)                                                  cAATTTATGAAGATGACAATAATTATAGATATATGAGAACATAAATGAAAATAATTATCATTACAGTAATCATTTGTACTTTGTATTGATGAGGGATGAAATGTTATATAATATACCTTGTCGAATTTATATCCTTTCCa                                                                                   <  1:989262/140‑1 (MQ=255)                                                               tGACAATAATTATAGATATATGAGAACATAAATGAAAATAATTATCATTACAGTAATCATTTGTACTTTGTATTGATGAGGGATGAAATGTTATATAATATACCTTGTCGAATTTATATCCTTTCCACTCTGTCATTATg                                                                      >  1:1218288/1‑140 (MQ=255)                                                                           taGATATATGAGAACATAAATGAAAATAATTATCATTACAGTAATCATTTGTACTTTGTATTGATGAGGGATGAAATGTTATATAATATACCTTGTCGAATTTATATCCTTTCCACTCTGTCATTATGCATTTCTGGGAt                                                          >  2:1851075/1‑140 (MQ=255)                                                                            aGATATATGAGAACATACATGAAAATAATTATCATCACAATAATCATTTGTACTTTGTATTGATGAGGGATGAAATGTTATATAATATACCTTGTCGAATTTATATCCTTTCCACTCTGTCATTATGCATTTCTGGGATa                                                         <  2:782572/140‑1 (MQ=255)                                                                             gATATATGAGAACATAAATGAAAATAATTATCATTACAGTAATCATTTGTACTTTGTATTGATGAGGGATGAAATGTTATATAATATACCTTGTCTATTTTATATCGTTTCCACTCTGGCATTATGCATTTCTGGGATAg                                                        >  2:1796772/1‑140 (MQ=255)                                                                             gATATATGAGAACATAAATGAAAATAATTATCATTACAGTAATCATTTGTACTTTGTATTGATGAGGGATGAAATGTTATATAATATACCTTGTCGAATTTATATCCTTTCCACTCTGTCATTATGCATTTCTGGGATAg                                                        >  2:1710202/1‑140 (MQ=255)                                                                             gATATATGAGAACATAAATGAAAATAATTATCATTACAGTAATCATTTGTACTTTGTATTGATGAGGGATGAAATGTTATATAATATACCTTGTCGAATTTAAATCCTTTCCACACTCTCATTACGCATTTCTGGGATAg                                                        >  1:782572/1‑140 (MQ=255)                                                                               tatatGAGAACATAAATGAAAATAATTATCATTACAGTAATCATTTGTACTTTGTATTGATGAGGGATGAAATGTTATATAATATACCTTGTCGAATTTATATCCTTTCCACTCTGTCATTATGCATTTCTGGGATAGtt                                                      >  1:979192/1‑140 (MQ=255)                                                                                atatGAGAACATAAATGAAAATAATTATCATTACAGTAATCATTTGTACTTTGTATTGATGAGGGATGAAATGTTATATAATATACCTTGTTGAATTTATATCCCTTCAACTCTGTCCTTATGCATTTCTGGGATAGttt                                                     >  2:743311/1‑140 (MQ=255)                                                                                atatGAGAACATAAATGAAAATAATTATCATTACAGTAATCATTTGTACTTTGTATTGATGAGGGATGAAATGTTATATAATATACCTTGTCGAATTTATATCCTTTCCACTCTGTCATTATGCATTTCTggg                                                            >  1:1312136/1‑133 (MQ=255)                                                                                atatGAGAACATAAATGAAAATAATTATCATTACAGTAATCATTTGTACTTTGTATTGATGAGGGATGAAATGTTATATAATATACCTTGTCGAATTTATATCCTTTCCACTCTGTCATTATGCATTTCTggg                                                            <  2:1312136/133‑1 (MQ=255)                                                                                     agaACATAAATGAAAATAATTATCATTACAGTAATCATTTGTACTTTGTATTGATGAGGGATGAAATGTTATATAATATACCTTGTCGAATTTATATCCTTTCCACTCTGTCATTATGCATTTCTGGGATAGTTTCTACt                                                >  2:1287948/1‑140 (MQ=255)                                                                                         cATAAATGAAAATAATTATCATTACAGTAATCATTTGTACTTTGTATTGATGAGGGATGGAATGTTATATAATATACCTTGTCGAATTTATGTCCTTTCCTCTCTGTCTTTATGCATGTCTGAGATAGTTTCTACTGCaa                                            >  2:528771/1‑140 (MQ=255)                                                                                           tAAATGAAAATAATTATCATTACAGTAATCATTTGTACTTTGTATTGATGAGGGATGAAATGTTATATAATATACCTTGTCGAATTTATATCCTTTCCACTCTGTCATTATGCATTTCTGGGATAGTTTCTACTGCAAcc                                          >  2:1520130/1‑140 (MQ=255)                                                                                                               tACAGTAATCATTTGTACTTTGTATTGATGAGGGATGAAATGTTATATAATATACCTTGTCGAATTTATATCCTTTCCACTCTGTCATTATGCATTTCTGGGATAGTTTCTACTGCAACCGCAACTTCTTCAGAAACaaa                      >  2:1334098/1‑140 (MQ=255)                                                                                                                    tAATCATTTGTACTTTGTATTGATGAGGGATGAAATGTTATATAATATACCTTGTCGAATTTATATCCTTTCCACTCTGTCATTATGCATTTCTGGGATAGTTTCTACTGCAACCGCAACTTCTTCAGAAACAAAAATCa                 >  2:330317/1‑140 (MQ=255)                                                                                                                     aaTCATTTGTACTTTGTATTGATGAGGGATGAAATGTTATATAATATACCTTGTCGAATTTATATCCTTTCCACTCTGTCATTATGCATTTCTGGGATAGTTTCTACTGCAACCGCAACTTCTTCAGAAACAAAAATCAg                <  1:330317/140‑1 (MQ=255)                                                                                                                      aTCATTTGTACTTTGTATTGATGAGGGATGAAATGTTATATAATATACCTTGTCGAATTTATATCCTTTCCACTCTGTCATTATGCATTTCTGGGATAGTTTGTACTGCAGCCGTAAGTTCTTCGGAAACAAAAATCAGc               >  1:532099/1‑140 (MQ=255)                                                                                                                              gaCTTTGTATTGATGTGGGATGCAATGTTCTATAATAGACCTTGTCGAATTTATATCCTTTCCACTCTGTCATTATGCATTTCTGGGATAGTTTCTACTGCAACCGCAACTTCTTCAGAAACAAAAATCAGCAACGAaga       <  1:576361/139‑1 (MQ=255)                                                                                                                                  ttGTATTGATGAGGGATGAAATGTTATATAATATACCTTGTCGAATTTATATCCTTTCCACTCTGTCATTATGCATTTCTGGGATAGTTTCTACTGCAACCGCAACTTCTTCAGAAACAAAAATCAGCAACGAAGAGACg   >  1:1523444/1‑140 (MQ=255)                                                                                                                                   tGTATTGATGAGGGATGAAATGTTATATAATATACCTTGTCGAATTTATATCCTTTCCACTCTGTCATTATGCATTTCTGGGATAGTTTCTACTGCAACCGCGACTTCTTGAGAGACAAAAATCAGCAACGAAGAGACGc  >  1:568303/1‑140 (MQ=255)                                                                                                                                         |                                                                                                                                      AAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATGAAGATGACAATAATTATAGATATATGAGAACATAAATGAAAATAATTATCATTACAGTAATCATTTGTACTTTGTATTAATGAGGGATGAAATGTTATATAATATACCTTGTCGAATTTATATCCTTTCCACTCTGTCATTATGCATTTCTGGGATAGTTTCTACTGCAACCGCAACTTCTTCAGAAACAAAAATCAGCAACGAAGAGACGC  >  NC_000913/2068486‑2068755

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