Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,191,696	G→A	*355* (TAG→TAA)	yqiI →	fimbrial protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,191,696	0	G	A	100.0%	149.0 / NA	47	*355* (TAG→TAA)	yqiI	fimbrial protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (24/23);  total (24/23)

GCTTCTTTTTCATCGTCTAGTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAGTATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTGTTGACAATACCAGGCATAACGTGCGCAAAACCAGGTTCTATGTTCTTCA  >  NC_000913/3191565‑3191826                                                                                                                                    |                                                                                                                                   gCTTCTTTTTCATCGTCTAGTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTaa                                                                                                                             >  2:1000183/1‑139 (MQ=255) gCTTCTTTTTCATCGTCTAGTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTaa                                                                                                                             >  1:490449/1‑139 (MQ=255)     cttTTTCATCGTCTAGTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACa                                                                                                                         <  1:989195/139‑1 (MQ=255)     cttTTTCATCGTCTAGTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACa                                                                                                                         <  1:860203/139‑1 (MQ=255)            aTCGTCTAGTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTa                                                                                                                  >  1:328609/1‑139 (MQ=255)                       cgcCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCAt                                                                                                       <  2:457841/139‑1 (MQ=255)                        gcCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATc                                                                                                      >  2:147407/1‑139 (MQ=255)                          cTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAg                                                                                                    <  2:490449/139‑1 (MQ=255)                           ttCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAgg                                                                                                   <  1:1000183/139‑1 (MQ=255)                                tCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGt                                                                                              <  2:535129/139‑1 (MQ=255)                                     gCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCAt                                                                                                       >  2:471626/1‑125 (MQ=255)                                     gCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCAt                                                                                                       <  1:471626/125‑1 (MQ=255)                                      cACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATc                                                                                        >  1:1004043/1‑139 (MQ=255)                                      cACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATc                                                                                        >  2:211535/1‑139 (MQ=255)                                            tCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGct                                                                                  <  2:81631/139‑1 (MQ=255)                                            tCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGct                                                                                  >  1:184121/1‑139 (MQ=255)                                              aGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGctct                                                                                <  1:536993/139‑1 (MQ=255)                                                     aaCTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATAc                                                                                                                          <  2:1030212/90‑1 (MQ=255)                                                     aaCTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATAc                                                                                                                          >  1:1030212/1‑90 (MQ=255)                                                     aaCTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCa                                                                         <  1:739181/139‑1 (MQ=255)                                                     aaCTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCa                                                                         >  2:739181/1‑139 (MQ=255)                                                     aaCTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCa                                                                         <  1:211535/139‑1 (MQ=255)                                                     aaCTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGATATATCAGTACTCTAACTCCa                                                                         >  1:511711/1‑139 (MQ=255)                                                          aGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGctct                                                                    >  1:915086/1‑139 (MQ=255)                                                                      gCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGcatc                                                        >  2:421101/1‑139 (MQ=255)                                                                      gCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAAACCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGcatc                                                        <  2:664364/139‑1 (MQ=255)                                                                          aaaCGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCt                                                    <  2:184121/139‑1 (MQ=255)                                                                                    tgAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTGTTGACAATa                                          >  2:654776/1‑139 (MQ=255)                                                                                     gAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTGTTGACAATAc                                         <  1:174664/139‑1 (MQ=255)                                                                                      aGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTGTTGACAATAcc                                        >  1:704002/1‑139 (MQ=255)                                                                                          aCCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTGTTGACAATACCAGGc                                    >  2:990257/1‑139 (MQ=255)                                                                                            cGTTTAGTAAAGACGTACTCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTGTTGACAATACCAGGCAt                                  >  1:660727/1‑139 (MQ=255)                                                                                              tttAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTg                                                                                    <  2:460630/87‑1 (MQ=255)                                                                                              tttAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTg                                                                                    >  1:460630/1‑87 (MQ=255)                                                                                                        aCGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTGTTGACAATAc                                         >  2:242545/1‑120 (MQ=255)                                                                                                        aCGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTGTTGACAATAc                                         <  1:242545/120‑1 (MQ=255)                                                                                                        aCGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTGTTGACAATACCAGGCATAACGTGCGCaaa                      >  1:1004836/1‑139 (MQ=255)                                                                                                                  gAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGcatcat                                                      >  1:188522/1‑97 (MQ=255)                                                                                                                  gAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGcatcat                                                      <  2:188522/97‑1 (MQ=255)                                                                                                                   aaaaTCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTGTTGACAATACCAGGCATAACGTGCGCAAAACCAGGTTCTa           <  1:771242/139‑1 (MQ=255)                                                                                                                    agaTCACCTATAACTAATATCTAATACAAACACTAACACGGGCCATCACGCCCGTTATATCAGTGCTCTAACTCCTGCTCTCTTGCCTGCATCATCTGTTGACAATACCAGGCATAACGTGCGCAAAACCAGGTTCTAt          <  2:511711/137‑1 (MQ=255)                                                                                                                        cACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTGTTGACAATACCAGGCATAACGTGCGCAAAACCAGGTTCTATGttc      >  2:273998/1‑139 (MQ=255)                                                                                                                         aCCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCt                                                    <  2:43263/92‑1 (MQ=255)                                                                                                                         aCCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCt                                                    >  1:43263/1‑92 (MQ=255)                                                                                                                            tataACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTGTTGACAATACCAGGCATAACGTGCGCAAAACCAGGTTCTATGTTCTTCa  <  2:1004836/139‑1 (MQ=255)                                                                                                                             ataACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTGTTGACAATAc                                         >  2:220827/1‑99 (MQ=255)                                                                                                                             ataACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTGTTGACAATAc                                         <  1:220827/99‑1 (MQ=255)                                                                                                                                    |                                                                                                                                   GCTTCTTTTTCATCGTCTAGTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAGTATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTGTTGACAATACCAGGCATAACGTGCGCAAAACCAGGTTCTATGTTCTTCA  >  NC_000913/3191565‑3191826

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: ‑