Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,019,827	T→C	Y76C (TAC→TGC)	ompA ←	outer membrane porin A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,019,827	0	T	C	100.0%	146.0 / NA	42	Y76C (TAC→TGC)	ompA	outer membrane porin A
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (20/22);  total (20/22)

CACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTCGTAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTACTGGGACCAGCCCAGTTTAGCAC  >  NC_000913/1019696‑1019962                                                                                                                                    |                                                                                                                                        cACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTCGCAACCCAtt                                                                                                                                 >  2:650369/1‑140 (MQ=255)   cccAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTCGCAACCCATTTc                                                                                                                               >  2:69280/1‑140 (MQ=255)   cccAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTCGCAACCCATTTc                                                                                                                               <  1:441685/140‑1 (MQ=255)         cGAGTGTAGATGTCCAGGTCGTAAGTGATTGAGTAGCCCAGTTTAGCGTTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTCGCAACCCATTTCAAAGcc                                                                                                                         <  2:304520/140‑1 (MQ=255)              gtAGATGTCCAGGTCGTCAGTGATTTGGTAACCCAGTCTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTCGCAACCCATTTCAAAGCCAACAt                                                                                                                    <  1:277248/140‑1 (MQ=255)              gtAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTCGCAACCCATTTCAAAGCCAACAt                                                                                                                    <  2:1227/140‑1 (MQ=255)                   tGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTCGCAACCCATTTCAAAGCCAACATACggg                                                                                                               <  2:523188/140‑1 (MQ=255)                   tGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTCGCAACCCATTTCAAAGCCAACATACggg                                                                                                               <  1:445703/140‑1 (MQ=255)                    gTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTCGCAACCCATTTCAAAGCCAACATACGGGt                                                                                                              >  1:453841/1‑140 (MQ=255)                     tCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTAGTACGGCATACGACCTAAGCAGACGCAACCCATTTCTAAGCCAACATACGGGtt                                                                                                             >  2:263531/1‑140 (MQ=255)                      ccAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTATCCAGTCGCAACACATTTCAAATCCAACATACGGGTTa                                                                                                            >  1:352548/1‑140 (MQ=255)                      ccAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTCGCAACCCATTTCAAAGCCAACATACGGGTTa                                                                                                            >  1:681646/1‑140 (MQ=255)                       cAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTAAGCAGTCGCAACCGATTTCAAAGGCAACATACGGGTTaa                                                                                                           >  1:193849/1‑140 (MQ=255)                       cAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTCGCAACCCATTTCAAAGCCAACATACGGGTTaa                                                                                                           >  2:74729/1‑140 (MQ=255)                         ggTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATCCGACCTAAACAGCCGCAACCCATTTCAAAGCCAACATACGGGTTAAcc                                                                                                         >  2:268674/1‑140 (MQ=255)                               cAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTCGCAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTaa                                                                                                   >  2:538501/1‑140 (MQ=255)                                            aCCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTCGCAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCAc                                                                                      <  1:74729/140‑1 (MQ=255)                                                aGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTCGCAACCCAATTCAAAGCTAACGTACAGGTTAACCAGGTAACCACCAAAAGCACCAgc                                                                                  >  2:237417/1‑140 (MQ=255)                                                 gTTTAGCGTTCAGTTGAGCGCCATGAGCTTTTTTTGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTCGCAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAgcg                                                                                 <  2:352548/140‑1 (MQ=255)                                                       cGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTCGCAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGt                                                                           >  2:43478/1‑140 (MQ=255)                                                            aGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTCGCAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGtt                                                                      >  2:639625/1‑140 (MQ=255)                                                                      ccTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTCGCAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAg                                                                            <  2:475481/124‑1 (MQ=255)                                                                      ccTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTCGCAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAg                                                                            >  1:475481/1‑124 (MQ=255)                                                                       cTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTCGCAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCg                                                           <  1:69280/140‑1 (MQ=255)                                                                          aGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTCGCAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGc                                                        <  2:453841/140‑1 (MQ=255)                                                                                gTATGCACCGTTTTCAACGCTGCCTTTGTATGGCATACGACCTAACCAGTCGCAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAAACACCAAAAACACGAGCGCCCAGTTGGTTTTGATGGGTCGGGCCAttgt                                                  >  2:211653/1‑140 (MQ=255)                                                                                  aTGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTCGCAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCAttgttg                                                <  1:43478/140‑1 (MQ=255)                                                                                  aTGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTCGCAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCAttgttg                                                <  1:268674/140‑1 (MQ=255)                                                                                               aaCGCTGCCTTTGTACGGCATACGACCTAACCAGTCGCAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAg                                   <  1:169734/140‑1 (MQ=255)                                                                                                 cGCTGCCTTTGTACGTCATACGTCCGAACAAGTCGCAGCCCATTTCAAAGCCAACATACGGGTTAACCTGGTAGCCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAgtg                                 <  1:263531/140‑1 (MQ=255)                                                                                                    tGCCTTTGTACGGCATACGACCTAACCAGTCGCAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCa                              <  1:639625/140‑1 (MQ=255)                                                                                                     gCCTTTGTACGGCATACGACCTAACCAGTCGCAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCAt                             <  1:538501/140‑1 (MQ=255)                                                                                                                gCATACGACCTAACCAGTCGCAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTACTGGGAc                  >  2:19509/1‑140 (MQ=255)                                                                                                                    aCGACCTAACCAGTCGCAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTACTGGGACCAGc              >  2:174605/1‑140 (MQ=255)                                                                                                                      gACCTAACCAGTCGCAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTACTGGGACCAGccc            <  1:174605/140‑1 (MQ=255)                                                                                                                      gACCTAACCAGTCGCAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTACTGGGACCAGccc            <  2:60547/140‑1 (MQ=255)                                                                                                                      gACCTAACCAGTCGCAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTACTGGGACCAGccc            <  1:650369/140‑1 (MQ=255)                                                                                                                      gACCTAACCAGTCGCAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTACTGGGACCAGccc            <  2:131840/140‑1 (MQ=255)                                                                                                                           aaCCAGTCGCAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTACTGGGACCAGCCCAGttt       >  1:47818/1‑140 (MQ=255)                                                                                                                           aaCCAGTCGCAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTACTGGGACCAGCCCAGttt       >  1:716497/1‑140 (MQ=255)                                                                                                                           aaCCAGTCGCAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTACTGGGACCAGCCCAGttt       >  1:650193/1‑140 (MQ=255)                                                                                                                             ccAGTCGCAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTACTGGGACCAGCCCAGTTTAg     <  2:650193/140‑1 (MQ=255)                                                                                                                                gTCGCAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTACTGGGACCAGCCCAGTTTAGCAc  <  1:362083/140‑1 (MQ=255)                                                                                                                                    |                                                                                                                                        CACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTCGTAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTACTGGGACCAGCCCAGTTTAGCAC  >  NC_000913/1019696‑1019962

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