Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,019,886	C→A	L56L (CTG→CTT)	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,886	0	C	A	100.0%	104.8 / NA	31	L56L (CTG→CTT)	ompA	outer membrane porin A
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (16/15);  total (16/15)

GAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTCGTAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTACTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTGCGCTACGGTAGCGAAAC  >  NC_000913/1019759‑1020013                                                                                                                                |                                                                                                                                gAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCTTTTGTACGGCATACGACCTAACCAGTCGTAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCAAGTTGGTTTTCa                                                                                                                     <  2:567387/140‑1 (MQ=255)    cGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTCGTAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCAAGTTGGTTTTCATgg                                                                                                                  >  2:308084/1‑140 (MQ=255)      cccTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTCGTAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCAAGTTGGTTTTCATGGGt                                                                                                                <  2:153412/140‑1 (MQ=255)               ttGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTCGTAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCAAGTTGGTTTTCATGGGTCGGGCCAtt                                                                                                       >  1:329342/1‑140 (MQ=255)                   aTGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTCGTAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCAAGTTGGTTTTCATGGGTCGGGCCAttgttg                                                                                                   <  2:409336/140‑1 (MQ=255)                       aCCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTCGTAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCAAGTTGGTTTTCATGGGTCgg                                                                                                             <  1:10352/126‑1 (MQ=255)                       aCCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTCGTAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCAAGTTGGTTTTCATGGGTCgg                                                                                                             >  2:10352/1‑126 (MQ=255)                          gTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTCGTAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCAAGTTGGTTTTCATGGGTCGGGCCATTGTTGTtgatga                                                                                            >  2:299493/1‑140 (MQ=255)                             ttCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTCGTAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCAAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAAc                                                                                         <  2:433403/140‑1 (MQ=255)                                aaCGCTGCCTTTGTACGGCATACGACCTAACCAGTCGTAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCAAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAg                                                                                      >  1:405943/1‑140 (MQ=255)                                     tGCCTTTGTACGGCATACGACCTAACCAGTCGTAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCAAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCa                                                                                 <  2:329342/140‑1 (MQ=255)                                       ccTTTGTACGGCATACGACCTAACCAGTCGTAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCAAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATg                                                                               <  1:677776/140‑1 (MQ=255)                                            gTACGGCATACGACCTAACCAGTCGTAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCAAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTACt                                                                          <  1:299493/140‑1 (MQ=255)                                                     aCGACCTAACCAGTCGTAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCAAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTACTGGGACCAGc                                                                 >  2:81956/1‑140 (MQ=255)                                                      cGACCTAACCAGTCGTAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCAAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTACTGGGACCAGcc                                                                <  2:405943/140‑1 (MQ=255)                                                       gACCTAACCAGTCGTAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCAAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTACTGGGACCAGccc                                                               >  2:643404/1‑140 (MQ=255)                                                            aaCCAGTCGTAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCAAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTACTGGGACCAGCCCAGttt                                                          >  1:649739/1‑140 (MQ=255)                                                            aaCCAGTCGTAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCAAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTACTGGGACCAGCCCAGttt                                                          >  2:420142/1‑140 (MQ=255)                                                            aaCCAGTCGTAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCAAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTACTGGGACCAGCCCAGttt                                                          >  2:10203/1‑140 (MQ=255)                                                              ccAGTCGTAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCAAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTACTGGGACCAGCCCAGTTTAg                                                        <  2:121283/140‑1 (MQ=255)                                                                 gTCGTAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCAAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTACTGGGACCAGCCCAGTTTAGCAc                                                     >  2:711944/1‑140 (MQ=255)                                                                   cGTAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCAAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTACTGGGACCAGCCCAGTTTAGCACCa                                                   >  2:415201/1‑140 (MQ=255)                                                                       aCCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCAAGTTGGTTTTCATGGGTAGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTACTGGGACCAGCCCAGTTTAg                                                        <  2:265322/131‑1 (MQ=255)                                                                       aCCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCAAGTTGGTTTTCATGGGTAGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTACTGGGACCAGCCCAGTTTAg                                                        >  1:265322/1‑131 (MQ=255)                                                                          cATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCAAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTACTGGGACCAGCCCAGTTTAGCACCAGTGTAcc                                            <  1:415201/140‑1 (MQ=255)                                                                            tttCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCAAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTACTGGGACCAGCCCAGTTTAGCAcc                                                    >  1:667542/1‑130 (MQ=255)                                                                            tttCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCAAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTACTGGGACCAGCCCAGTTTAGCAcc                                                    <  2:667542/130‑1 (MQ=255)                                                                            tttCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCAAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTACTGGGACCAGCCCAGTTTAGCACCAGTGTACCAg                                          <  1:10203/140‑1 (MQ=255)                                                                                                        gaaaCCAACAAAAGCACCAGCGCCAAGTTGGTTTTCAGGGGTCGGGCTTTTGTTGTTGATGAAACTGGTGTCATGGTACTGGGACCAGCCCAGTTTAGCACCAGTGTACCGGGTTTTATCTTTCGGGGCGGCCTGCGCTa              <  2:649739/138‑1 (MQ=255)                                                                                                          aaCCACCAAAAGCACCAGCGCCAAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTACTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTGCGCTACg            <  2:302594/140‑1 (MQ=255)                                                                                                           accaccAAAAGCACCAGCGCCAAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTACTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTGCGCTACgg           >  1:449279/1‑140 (MQ=255)                                                                                                             caccaAAAGCACCAGCGCCAAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTACTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTGCGCTACGGTa         >  2:817557/1‑140 (MQ=255)                                                                                                                   aaGCACCAGCGCCAAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTACTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTGCGCTACGGAAGCGaaa   >  2:367807/1‑140 (MQ=255)                                                                                                                    aGCACCAGCGCCAAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTACTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTGCGCTACGGTAGCGAAAc  <  1:367807/140‑1 (MQ=255)                                                                                                                                |                                                                                                                                GAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTCGTAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTACTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTGCGCTACGGTAGCGAAAC  >  NC_000913/1019759‑1020013

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