Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,101,553	G→A	pseudogene (192/349 nt)	wbbL ←	interrupted rhamnosyltransferase WbbL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,101,553	0	G	A	100.0%	110.3 / NA	30	pseudogene (192/349 nt)	wbbL	interrupted rhamnosyltransferase WbbL
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (16/14);  total (16/14)

ATAAGTATACAAAGTACTTATTAAAGTATAAATAGCTTATCCATGCTTATATGCTTACGGCTTTATATTACGGGTGAAAAACTGATGAAATTCGATCAAAGTTGCGATTTGATAAAATACGTTTTCTGGCTAAATATCTAAAAGTACTTTTTAAGTGCCATCTGAAGGCTTTTGAAAAAAAACTTCGATTGTCATGATGAGCATAATGTATCGCATGAAAAGCGGGAACATAATGAAGTCTGACACCAGCCAGGCTAAGCCTCAAGCACAGGTCAATATCTTCACAGTACATAAAGTAACCTTGATCGAAGCCATTTACACGCACAAAATCTGAAAAACGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCG  >  NC_000913/2101329‑2101781                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                     aTAAGTATACAAAGTACTTATTAAAGTATAAATAGCTTATCCATGCTTATATGCTTACGGCTTTATATTACGGGTGAAAAACTGATGAAATTCGATCAAAGTTGCGATTTGATAAAATACGTTTTCTGGCTAAATATCTAAAAGTACTTTTTAAGTGCCATCTGAAGGCTTTTGAAAAAAAACTTCGATTGTCATGATGAGCATAATGTATCGCATGAAAAGCGAGAACATAATGAAGTCTGACAccagcc                                                                                                                                                                                                            <  1:225648/251‑1 (MQ=255)    aGTATACAAAGTACTTATTAAAGTATAAATAGCTTATCCATGCTTATATGCTTACGGCTTTATATTACGGGTGAAAAACTGATGAAATTCGATCAAAGTTGCGATTTGATAAAATACGTTTTCTGGCTAAATATCTAAAAGTACTTTTTAAGTGCCATCTGAAGGCTTTTGAAAAAAAACTTCGATTGTCATGATGAGCATAATGTATCGCATGAAAAGCGAGAACATAATGAAGTCTGACACCAGCCAgg                                                                                                                                                                                                         <  2:4785/251‑1 (MQ=255)                                                  tatGCTTACGGCTTTATATTACGGGTGAAAAACTGATGAAATTCGATCAAAGTTGCGATTTGATAAAATACGTTTTCTGGCTAAATATCTAAAAGTACTTTTTAAGTGCCATCTGAAGGCTTTTGAAAAAAAACTTCGATTGTCATGATGAGCATAATGTATCGCATGAAAAGCGAGAACATAATGAAGTCTGACACCAGCCAGGCTAAGCCTCAAGCACAGGTCAATATCTTCACAGTACATAAAGTAAc                                                                                                                                                           >  2:208423/1‑251 (MQ=255)                                                          cGGCTTTATATTACTGGTGAAAAACTGATGAAATTCGATCAAAGTTGCGATTTTATAAAATACGTTTTCTGGCTAAATATCTAAAAGTACGTTTTAAGTGCCATCTGAAGGCTTTTGAAAAAAAACTTCGATTTGCATGATGAGCATAATGTATCGCATGAAAAGCGAGAACATAATGAAGTCTGACACCAGCCAGGCTAAGCCTCAAGCACAGGTCAATATCTTCACAGTACAGAAAGTAACCTTGATCg                                                                                                                                                   <  2:350295/251‑1 (MQ=255)                                                                 atatTACGGGTGAAAAACTGATGAAATTCGATCAAAGTTGCGATTTGATAAAATACGTTTTCTGGCTAAATATCTAAAAGTACTTTTTAAGTGCCATCTGAAGGCTTTTGAAAAAAAACTTCGATTGTCATGATGAGCATAATGTATCGCATGAAAAGCGAGAACATAATGAAGTCTGACACCAGCCAGTCTAAGCCTCAAGCACAGGTCAATATCTTCACAGTACATAAAGTCACCTTGATCGAAGCCAt                                                                                                                                            >  1:55321/1‑251 (MQ=255)                                                                        gggTGAAAAACTGATGAAATTCGATCAAAGTTGCGATTTGATAAAATACGTTTTCTGGCTAAATATCTAAAAGTACTTTTTAAGTGCCATCTGAAGGCTTTTGAAAAAAAACTTCGATTGTCATGATGAGCATAATGTATCGCATGAAAAGCGAGAACATAATGAAGTCTGACACCAGCCAGGCTAAGCCTCAAGCACAGGTCAATATCTTCACAGTACATAAAGTAACCTTGATCGAAGCCATTTACACg                                                                                                                                     <  1:188254/251‑1 (MQ=255)                                                                        gggTGAAAAACTGATGAAATTCGATAAAAGTTGCGATTTGATAAAATACGTTTTCTGGCTAAATATCTAAAAGTACTTTTTAAGTGCCATCTGAAGGCTTTTGAAAAAAAACTTCGATTGTCATGATGAGCATAATGTATCGCATGAAAAGCGAGAACATAATGAAGTCTGACACCAGCCAGGCTAAGCCTCAAGCACAGGTCAATATCTTCACAGTACATAAAGTAACCTTGATCGAAGCCATTTACACg                                                                                                                                     >  1:323204/1‑251 (MQ=255)                                                                          gTGAAAAACTGATGAAATTCGATCAAAGTTGCGATTTGATAAAATACGTTTTCTGGCTAAATATCTAAAAGTACTTTTTAAGTGCCATCTGAAGGCTTTTGAAAAAAAACTTCGATTGTCATGATGAGCATAATGTATCGCATGAAAAGCGAGAACATAATGAAGTCTGACACCAGCCAGGCTAAGCCTCAAGCACAGGTCAATATCTTCACAGTACATAAAGTAACCTTGATCGAAGCCATTTACACGca                                                                                                                                   <  1:332612/251‑1 (MQ=255)                                                                                 aCTGATGAAATTCGATCAAAGTTGCGATTTGATAAAATACGTTTTCTGGCTAAATATCTAAAAGTACTTTTTAAGTGCCATCTGAAGGCTTTTGAAAAAAAACTTCGATTGTCATGATGAGCATAATGTATCGCATGAAAAGCGAGAACATAATGAAGTCTGACACCAGCCAGGCTAAGCCTCAAGCACAGGTCAATATCTTCACAGTACATAAAGTAACCTTGATCGAAGCCATTTACACGCACAAAATc                                                                                                                            <  2:204111/251‑1 (MQ=255)                                                                                                      ttGCGATTTGATAAAATACGTTTTCTGGCTAAATCTCTAAAAGTACTTTTTAAGTGCCATCTGAAGGCTTTTGACAAAAAACTTCGATTGTCATGATGAGCATAATGTATCGCATGAAAAGCGAGAACATAATGAAGTCTGACACCAGCCAGGCTAAGCCTCAAGCACAGGTCAATATCTTCACAGTACATAAAGTAACCTTGATCGAAGCCATTTACACGCACCAAATCTGAAAAACGTACCAGCATaa                                                                                                        >  2:252898/1‑250 (MQ=255)                                                                                                       tGCGATTTGATAAAATACGTTTTCTGGCTAAATATCTAAAAGTACTTTTTAAGTGCCATCTGAAGGCTTTTGAAAAAAAACTTCGATTGTCATGATGAGCATAATGTATCGCATGAAAAGCGAGAACATAATGAAGTCTGACACCAGCCAGGCTAAGCCTCAAGCACAGGTCAATATCTTCACAGTACATAAAGTAACCTTGATCGAAGCCATTTACACGCACAAAATCTGAAAAACGTACCAGCATaaa                                                                                                       <  2:320035/250‑1 (MQ=255)                                                                                                            tttGATAAAATACGTTTTCTGGCTAAATATCTAAAAGTACTTTTTAAGTGCCATCTGAAGGCTTTTGAAAAAAAACTTCGATTGTCATGATGAGCATAATGTATCGCATGAAAAGCGAGAACATAATGAAGGCTGACACCAGCAAGGCTAAGCCTCAAGCACAGGTCAATATCTTCACAGTACATAAAGTAACCTTGATCGAAGCCATTTACACGCAGAAAATCTGAAAAACGTACCAGCATAAATGATcc                                                                                                 >  1:207435/1‑251 (MQ=255)                                                                                                                               tGGCTAAATATCTAAAAGTACTTTTTAAGTGCCATCTGAAGGCTTTTGAAAAAAAACTTCGATTGTCATGATGAGCATAATGTATCGCATGAAAAGCGAGAACATAATGAAGTCTGACACCAGCCAGGCTAAGCCTCAAGCACAGGTCAATATCTTCACAGTACATAAAGTAACCTTGATCGAAGCCATTTACACGcaca                                                                                                                                 <  1:328554/200‑1 (MQ=255)                                                                                                                               tGGCTAAATATCTAAAAGTACTTTTTAAGTGCCATCTGAAGGCTTTTGAAAAAAAACTTCGATTGTCATGATGAGCATAATGTATCGCATGAAAAGCGAGAACATAATGAAGTCTGACACCAGCCAGGCTAAGCCTCAAGCACAGGTCAATATCTTCACAGTACATAAAGTAACCTTGATCGAAGCCATTTACACGcaca                                                                                                                                 >  2:328554/1‑200 (MQ=255)                                                                                                                               tGGCTAAATATCTAAAAGTACTTTTTAAGTGCCATCTGAAGGCTTTTGAAAAAAAACTTCGATTGTCATGATGAGCATAATGTATCGCATGAAAAGCGAGAACATAATGAAGTCTGACACCAGCCAGGCTAAGCCTCAAGCACAGGTCAATATCTTCACAGTACATAAAGTAACCTTGATCGAAGCCATTTACACGCACAAAATCTGAAAAACGTACCAGCATAAATGATCCTGCGCACCAATCAACAAcc                                                                              >  1:78857/1‑251 (MQ=255)                                                                                                                                 gCTAAATATCTAAAAGTACTTTTTAAGTGCCATCTGAAGGCTTTTGAAAAAAAACTTCGATTGTCATGATGAGCATAATGTATCGCATGAAAAGCGAGAACATAATGAAGTCTGACACCAGCCAGGCTAAGCCTCAAGCACAGGTCAATATCTTCACAGTACATAAAGTAACCTTGATCGAAGCCATTTACACGCACAAAATCTGAAAAACGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGt                                                                            >  2:150098/1‑251 (MQ=255)                                                                                                                                        atCTAAAAGTACTTTTTAAGTGCCATCTGAAGGCTTTTGAAAAAAAACTTCGATTGTCATGATGAGCATAATGTATCGCATGAAAAGCGAGAACATAATGAAGTCTGACACCAGCCAGGCTAAGCCTCAAGCACAGGTCAATATCTTCACAGTACATAAAGTAACCTTGATCGAAGCCATTTACACGCACAAAATCTGAAAAACGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGa                                                                      >  1:290188/1‑250 (MQ=255)                                                                                                                                                gTACTTTTTAAGTGCCATCTGAAGGCTTTTGAAAAAAAACTTCGATTGTCATGATGAGCATAATGTATCGCATGAAAAGCGAGAACATAATGAAGTCTGACACCAGCCAGGCTAAGCCTCAAGCACAGGTCAATATCTTCACAGTACATAAAGTAACCTTGATCGAAGCCATTTACACGCACAAAATCTGAAAAACGTAACAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATAct                                                             <  2:79405/251‑1 (MQ=255)                                                                                                                                                       ttAAGTGCCATCTGAAGGCTTTTGAAAAAAAACTTCGATTGTCATGATGAGCATAATGTATCGCATGAAAAGCGAGAACATAATGAAGTCTGACACCAGCCAGGCTAAGCCTCAAGCACAGGTCAATATCTTCACAGTACATAAAGTAACCTTGATCGAAGCCATTTACACGCACAAAATCTGAAAAACGTACCAGCATAAATGATCCTGCGCACCAATAAACAACCGTATCAGAATAGATActttcttt                                                       >  2:201455/1‑250 (MQ=255)                                                                                                                                                         aaGTGCCATCTGAAGGCTTTTGAAAAAAAACTTCGATTGTCATGATGAGCATAATGTATCGCATGAAAAGCGAGAACATAATGAAGTCTGACACCAGCCAGGCTAAGCCTCAAGCACAGGTCAATATCTTCACAGTACATAAAGTAACCTTGATCGAAGCCATTTACACGCACAAAATCTGAAAAACGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAgg                                                    <  1:1875/251‑1 (MQ=255)                                                                                                                                                         aaGTGCCATCTGAAGGCTTTTGAAAAAAAACTTCGATTGTCATGATGAGCATAATGTATCGCATGAAAAGCGAGAACATAATGAAGTCTGACACCAGCCAGGCTAAGCCTCAAGCACAGGTCAATATCTTCACAGTACATAAAGTAACCTTGATCGAAGCCAGTTACACGCACAAAATCTGAAAAACGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGAGTCTTTCTTTAgg                                                    <  2:238426/251‑1 (MQ=255)                                                                                                                                                             gCCATCTGAAGGCTTTTGAAAAAAAACTTCGATTGTCATGATGAGCATAATGTATCGCATGAAAAGCGAGAACATAATGAAGTCTGACACCAGCCAGGCTAAGCCTCAAGCACAGGTCAATATCTTCACAGTACATAAAGTAACCTTGATCGAAGCCATTTACACGCACAAAATCTGAAAAACGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAAtt                                                >  1:234181/1‑251 (MQ=255)                                                                                                                                                                        gCTTTTGAAAAAAAACTTCGATTGTCATGATGAGCATAATGTATCGCATGAAAAGCGAGAACATAATGAAGTCTGACACCAGCCAGGCTAAGCCTCAAGCACAGGTCAATATCTTCACAGTACATAAAGTAACCTTGATCGAAGCCATTTACACGCACAAAATCTGAAAAACGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTttatt                                     >  1:267100/1‑251 (MQ=255)                                                                                                                                                                                    aaaCTTCGATTGTCATGATGAGCATAATGTATCGCATGAAAAGCGAGAACATAATGAAGTCTGACACCAGCCAGGCTAAGCCTCAAGCACAGGTCAATATCTTCACAGTACATAAAGTAACCTTGATCGAAGCCATTTACACGCACAAAATCTGAAAAACGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCg                         >  2:255403/1‑251 (MQ=255)                                                                                                                                                                                              tGTCATGATGAGCAGAATGTATCGCATGAAAAGCGAGAACATAATGAAGGCTGACACCAGACAGGCTAAGCCTCAAGCACAGGTCAATATCTTCACAGTACATAAAGTAACCTTGATCGAAGCCATTTACACGCACAAAATCTGAAAAACGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTGTTGTTTTATTAAGGAAGGTGCGAATAAGCgg                <  2:323204/250‑1 (MQ=255)                                                                                                                                                                                                  atgatgAGCATAATGTATCGCATGAAAAGCGAGAACATAATGAAGTCTGACACCAGCCAGGCTAAGCCTCAAGCACAGGTCAATATCTTCACAGTACATAAAGTAACCTTGATCGAAGCCATTTACACGCACAAAATCTGAAAAACGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGaa            >  1:25406/1‑250 (MQ=255)                                                                                                                                                                                                    gatgaGCATAATGTATCGCATGAAAAGCGAGAACATAATGAAGTCTGACACCAGCCAGGCTAAGCCTCAAGCACAGGTCAATATCTTCACAGTACATAAAGTAACCTTGATCGAAGCCATTTACACGCACAAAATCTGAAAAACGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGa                        <  2:107730/236‑1 (MQ=255)                                                                                                                                                                                                    gatgaGCATAATGTATCGCATGAAAAGCGAGAACATAATGAAGTCTGACACCAGCCAGGCTAAGCCTCAAGCACAGGTCAATATCTTCACAGTACATAAAGTAACCTTGATCGAAGCCATTTACACGCACAAAATCTGAAAAACGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGa                        >  1:107730/1‑236 (MQ=255)                                                                                                                                                                                                    gatgaGCATAATGTATCGCATGAAAAGCGAGAACATAATGAAGTCTGACACCAGCCAGGCTAAGCCTCAAGCACAGGTCAATATCTTCACAGTACATAAAGTAACCTTGATCGAAGCCATTTACACGCACAAAATCTGAAAAACGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAAtt         >  2:113502/1‑251 (MQ=255)                                                                                                                                                                                                           aTAATGTATCGCATGAAAAGCGAGAACATAATGAAGTCTGACACCAGCCAGGCTAAGCCTCAAGCACAGGTCAATATCTTCACAGTACATAAAGTAACCTTGATCGAAGCCATTTACACGCACAAAATCTGAAAAACGTACCAGCATAAATGATCCTGCGCACCAATCAACAGCCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGGTTTAGGAAGGAAGGTGCTAATGAGCGGGTAAATGCTTCTCg  <  2:217004/251‑1 (MQ=255)                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                     ATAAGTATACAAAGTACTTATTAAAGTATAAATAGCTTATCCATGCTTATATGCTTACGGCTTTATATTACGGGTGAAAAACTGATGAAATTCGATCAAAGTTGCGATTTGATAAAATACGTTTTCTGGCTAAATATCTAAAAGTACTTTTTAAGTGCCATCTGAAGGCTTTTGAAAAAAAACTTCGATTGTCATGATGAGCATAATGTATCGCATGAAAAGCGGGAACATAATGAAGTCTGACACCAGCCAGGCTAAGCCTCAAGCACAGGTCAATATCTTCACAGTACATAAAGTAACCTTGATCGAAGCCATTTACACGCACAAAATCTGAAAAACGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCG  >  NC_000913/2101329‑2101781

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