Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	986,671	A→G	R104R (CGT→CGC)	ompF ←	outer membrane porin 1a (Ia,b,F)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	986,671	0	A	G	100.0%	173.2 / NA	48	R104R (CGT→CGC)	ompF	outer membrane porin 1a (Ia,b,F)
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (24/24);  total (24/24)

AGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGACGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAG  >  NC_000913/986535‑986790                                                                                                                                         |                                                                                                                        agTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGcg                                                                                                                       >  1:782044/1‑139 (MQ=255)   tCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGcgcg                                                                                                                     >  1:18637/1‑139 (MQ=255)                  gTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGttt                                                                                                      <  1:1046891/139‑1 (MQ=255)                        ccaccaAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGt                                                                                                >  1:606939/1‑139 (MQ=255)                           ccaAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAg                                                                                             >  2:1038479/1‑139 (MQ=255)                           ccaAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAg                                                                                             <  2:18637/139‑1 (MQ=255)                             aaaTTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAgcg                                                                                           >  1:399069/1‑139 (MQ=255)                                  cTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTc                                                                                      <  1:703749/139‑1 (MQ=255)                                   tGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAgc                                                                                            >  2:1029918/1‑132 (MQ=255)                                   tGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAgc                                                                                            <  1:1029918/132‑1 (MQ=255)                                          atatCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAgttgt                                                                              <  2:782044/139‑1 (MQ=255)                                                  gtAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTgg                                                                      <  1:363712/139‑1 (MQ=255)                                                   tAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGa                                                                     >  1:933628/1‑139 (MQ=255)                                                      cccAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCAAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCa                                                                                     <  2:443663/120‑1 (MQ=255)                                                      cccAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCAAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCa                                                                                     >  1:443663/1‑120 (MQ=255)                                                           tGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCa                                                                                     <  2:978669/115‑1 (MQ=255)                                                           tGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCa                                                                                     >  1:978669/1‑115 (MQ=255)                                                           tGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAgttgt                                                                              >  2:456200/1‑122 (MQ=255)                                                           tGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAgttgt                                                                              <  1:456200/122‑1 (MQ=255)                                                           tGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTgg                                                                      >  2:394942/1‑130 (MQ=255)                                                           tGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTgg                                                                      <  1:394942/130‑1 (MQ=255)                                                           tGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTatat                                                             >  1:720193/1‑139 (MQ=255)                                                                catAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAgc                                                                                            >  2:330637/1‑103 (MQ=255)                                                                catAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAgc                                                                                            <  1:330637/103‑1 (MQ=255)                                                                     aCCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGgcgc                                                                                                                      <  1:441480/72‑1 (MQ=255)                                                                     aCCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGgcgc                                                                                                                      >  2:441480/1‑72 (MQ=255)                                                                          acCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGcg                                                                                                                       <  1:21829/66‑1 (MQ=255)                                                                          acCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAgttgt                                                                              <  1:720071/107‑1 (MQ=255)                                                                          acCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAgttgt                                                                              >  2:720071/1‑107 (MQ=255)                                                                          acCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCAt                                                <  2:522557/137‑1 (MQ=255)                                                                          acCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCAt                                                >  1:522557/1‑137 (MQ=255)                                                                          acCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATaa                                              >  2:92406/1‑139 (MQ=255)                                                                           ccGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAgc                                                                                            >  2:533360/1‑92 (MQ=255)                                                                           ccGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAgc                                                                                            <  1:533360/92‑1 (MQ=255)                                                                           ccGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCa                                                        <  1:460439/128‑1 (MQ=255)                                                                           ccGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCa                                                        >  2:460439/1‑128 (MQ=255)                                                                               aGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATAACCGGt                                         <  1:1038479/139‑1 (MQ=255)                                                                                      gCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATAACCGGTCAGATCg                                  >  2:293253/1‑139 (MQ=255)                                                                                                     aCCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGAc                                                   >  2:527084/1‑107 (MQ=255)                                                                                                     aCCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGAc                                                   <  1:527084/107‑1 (MQ=255)                                                                                                        aaCGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATAACCGGTCAGATCGGAATTGATTTGAGTTTcc                <  2:606939/139‑1 (MQ=255)                                                                                                                      tAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAg  <  2:720193/139‑1 (MQ=255)                                                                                                                               cgAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCa                                                        >  1:9048/1‑76 (MQ=255)                                                                                                                               cgAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCa                                                        >  2:542432/1‑76 (MQ=255)                                                                                                                               cgAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCa                                                        <  2:9048/76‑1 (MQ=255)                                                                                                                               cgAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCa                                                        <  1:542432/76‑1 (MQ=255)                                                                                                                               cgAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATaa                                              >  2:600510/1‑86 (MQ=255)                                                                                                                               cgAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATaa                                              <  1:600510/86‑1 (MQ=255)                                                                                                                               cgAATGCCAGGCGCGTTTTGTTA‑CAGTTTGAGCGTCAGCGCCTTCaga                                                                                   <  2:119580/48‑1 (MQ=255)                                                                                                                               cgAATGCCAGGCGCGTTTTGTTA‑CAGTTTGAGCGTCAGCGCCTTCaga                                                                                   >  1:119580/1‑48 (MQ=38)                                                                                                                                         |                                                                                                                        AGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGACGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAG  >  NC_000913/986535‑986790

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