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%	151.4 / NA	48	*355* (TAG→TAA)	yqiI	fimbrial protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (23/25);  total (23/25)

TCTTTTTCATCGTCTAGTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAGTATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTGTTGACAATACCAGGCATAACGTGCGCAAAACCAGGTTCTATGTTCT  >  NC_000913/3191568‑3191823                                                                                                                                 |                                                                                                                                tcttTTTCATCGTCTAGTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATAc                                                                                                                       >  2:1640153/1‑139 (MQ=255)  cttTTTCATCGTCTAGTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACa                                                                                                                      <  2:1306816/139‑1 (MQ=255)     tttCATCGTCTAGTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCAGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAac                                                                                                                   >  2:229694/1‑139 (MQ=255)             tCTAGTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACa                                                                                                                      >  2:1235610/1‑128 (MQ=255)             tCTAGTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACa                                                                                                                      <  1:1235610/128‑1 (MQ=255)              cTAGTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTaaagcg                                                                                                          >  2:1236096/1‑136 (MQ=255)                       cTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGGAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAg                                                                                                 <  1:1526520/139‑1 (MQ=255)                        ttCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAgg                                                                                                <  1:1640153/139‑1 (MQ=255)                        ttCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAgg                                                                                                <  1:1252801/139‑1 (MQ=255)                          cGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAGCGGGCCATCAGGcc                                                                                              <  1:1236096/139‑1 (MQ=255)                                  gCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTatat                                                                                      <  1:229694/139‑1 (MQ=255)                                    aCCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCa                                                                                    <  2:384073/139‑1 (MQ=255)                                       aTTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATGACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTg                                                                                 >  1:1531862/1‑139 (MQ=255)                                             ggCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAg                                                                                   <  2:490003/131‑1 (MQ=255)                                             ggCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAg                                                                                   >  1:490003/1‑131 (MQ=255)                                                  aaCTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCa                                                                      >  2:62561/1‑139 (MQ=255)                                                  aaCTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTATAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCa                                                                      >  1:1079932/1‑139 (MQ=255)                                                       aGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAgg                                                                                                <  1:159955/108‑1 (MQ=255)                                                       aGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAgg                                                                                                >  2:159955/1‑108 (MQ=255)                                                       aGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGctct                                                                 >  1:1474244/1‑139 (MQ=255)                                                       aGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACATATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGctct                                                                 >  2:1129716/1‑139 (MQ=255)                                                              aCCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCt                                                          <  1:1129716/139‑1 (MQ=255)                                                                  ggCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTa                                                                                                                           >  1:1568303/1‑70 (MQ=255)                                                                  ggCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTa                                                                                                                           <  2:1568303/70‑1 (MQ=255)                                                                  ggCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGcat                                                      <  1:1610099/139‑1 (MQ=255)                                                                       aaaCGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCt                                                 <  2:1079932/139‑1 (MQ=255)                                                                             ttgttgAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGcat                                                      <  2:846144/128‑1 (MQ=255)                                                                             ttgttgAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGcat                                                      >  1:846144/1‑128 (MQ=255)                                                                              tgttgAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTGTTGACa                                          <  2:1474244/139‑1 (MQ=255)                                                                                ttgAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAg                                                                     <  2:496852/110‑1 (MQ=255)                                                                                ttgAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAg                                                                     >  1:496852/1‑110 (MQ=255)                                                                                 tgAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTGTTGACAATa                                       >  1:12487/1‑139 (MQ=255)                                                                                       aCCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTGTTGACa                                          <  2:1466762/130‑1 (MQ=255)                                                                                       aCCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTGTTGACa                                          >  1:1466762/1‑130 (MQ=255)                                                                                                     aCGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGcc                                                                                              >  2:1378383/1‑64 (MQ=255)                                                                                                     aCGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGcc                                                                                              <  1:1378383/64‑1 (MQ=255)                                                                                                     aCGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTGTTGACAATACCAGGCATAACGTGCGCaaa                   >  2:677752/1‑139 (MQ=255)                                                                                                     aCGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTGTTGACAATACCAGGCATAACGTGCGCaaa                   >  1:830313/1‑139 (MQ=255)                                                                                                     aCGTAATCCTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTGTTGACAATACCAGGCATAACGTGCGCaaa                   >  2:35600/1‑139 (MQ=255)                                                                                                             cTGAAAATCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTGTTGACAATACCAGGCATAACGTGCGCAAAACCAGGtt           <  2:12487/139‑1 (MQ=255)                                                                                                                aaaaTCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAg                                                                     >  2:1226401/1‑78 (MQ=255)                                                                                                                aaaaTCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAg                                                                     <  1:1226401/78‑1 (MQ=255)                                                                                                                aaaaTCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCt                                                 >  1:1275090/1‑98 (MQ=255)                                                                                                                aaaaTCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCt                                                 <  2:1275090/98‑1 (MQ=255)                                                                                                                  aaTCACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTGTTGACAATACCAGGCATAACGTGCGCAAAACCAGGTTCTATg      <  2:1531862/139‑1 (MQ=255)                                                                                                                     cACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTGTTGACa                                          <  2:1014374/100‑1 (MQ=255)                                                                                                                     cACCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTCATATCAGTGCTCTAACTCCAGGTCACTTGCCTGCATCATCTGTTGACa                                          >  1:1014374/1‑100 (MQ=255)                                                                                                                      aCCTATAACTAATATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTGTTGACAATACCAGGCATAACGTGCGCAAAACCAGGTTCTATGttct  <  1:677752/139‑1 (MQ=255)                                                                                                                                 |                                                                                                                                TCTTTTTCATCGTCTAGTGCGCCTTCGGTCACCGCACCATTCAGGGCAGAACTGAGTGCGAACCCGGCAGAAACGGTTGTTGAGGGACCGTTTAGTAAAGACGTAATCCTGAAAATCACCTATAACTAGTATCTAATACAAACACTAAAACGGGCCATCAGGCCCGTTATATCAGTGCTCTAACTCCAGCTCTCTTGCCTGCATCATCTGTTGACAATACCAGGCATAACGTGCGCAAAACCAGGTTCTATGTTCT  >  NC_000913/3191568‑3191823

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