Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,382,820	T→C	22.0%	Y292H (TAT→CAT)	ompG →	outer membrane porin G

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,382,820	0	T	C	22.0%	123.6 / 33.3	68	Y292H (TAT→CAT)	ompG	outer membrane porin G
Reads supporting (aligned to +/- strand):  ref base T (21/32);  new base C (7/8);  total (28/40)
Fisher's exact test for biased strand distribution p-value = 7.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.17e-01

ATTTTACGGTTATGATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTAGGCCGGA  >  NC_000913/1382723‑1382914                                                                                                  |                                                                                               aTTTTACGGTTATGATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTa                                                                                               <  1:2389276/99‑1 (MQ=255)     tACGGTTATGATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAg                                                                                          <  1:2064677/100‑1 (MQ=255)     tACGGTTATGATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAg                                                                                          <  1:2064001/100‑1 (MQ=255)       cGGTTATGATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGgt                                                                                        >  2:1344246/1‑100 (MQ=255)        ggTTATGATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGgtg                                                                                       <  2:355148/100‑1 (MQ=255)         gTTATGATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGCTCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGgtgt                                                                                      >  1:586117/1‑100 (MQ=255)            aTGATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCgg                                                                                   <  2:1675775/100‑1 (MQ=255)             tGATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGc                                                                                  >  1:1199214/1‑100 (MQ=255)             tGATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGc                                                                                  <  1:1589858/100‑1 (MQ=255)             tGATTTCCAGAACGGACTATCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGGGTCgg                                                                                   >  2:771307/1‑99 (MQ=255)              gATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCg                                                                                 <  2:1837360/100‑1 (MQ=255)               aTTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGt                                                                                >  2:1486571/1‑100 (MQ=255)                  tCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTaaa                                                                             >  2:799782/1‑100 (MQ=255)                   ccAGAACTGAATCTCCCCATCGCTAGAATACGCGTTCAAATTCCAGTATCACGACGAAGGCGACAATGATAAATTCCATCTTTCATCTGTCGGCGTAAAt                                                                            <  2:1191440/100‑1 (MQ=255)                   ccAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTaaa                                                                             >  1:1074104/1‑99 (MQ=255)                   ccAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAAt                                                                            <  2:1273407/100‑1 (MQ=255)                   ccAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAAt                                                                            >  2:2353807/1‑100 (MQ=255)                   ccAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAAt                                                                            <  2:2359847/100‑1 (MQ=255)                    cAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAAtt                                                                           <  1:1413433/100‑1 (MQ=255)                      gAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTAc                                                                         <  1:1890093/100‑1 (MQ=255)                      gAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTAc                                                                         <  1:1753243/100‑1 (MQ=255)                          ggACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGt                                                                     >  1:2322003/1‑100 (MQ=255)                           gACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGtt                                                                    >  1:339172/1‑100 (MQ=255)                           gACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGt                                                                     <  2:1099626/99‑1 (MQ=255)                             cTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTc                                                                   >  1:791283/1‑99 (MQ=255)                              tttCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTg                                                                 <  1:1593327/100‑1 (MQ=255)                                    tttCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATg                                                           <  2:612831/100‑1 (MQ=255)                                    tctCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTGCATCATGCAGTTGTCGGCGTAAATTACTCGTTCTGATAATg                                                           <  2:2064318/98‑1 (MQ=255)                                     ttCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTAGTTCTGATAATgg                                                          >  1:1763523/1‑100 (MQ=255)                                      tCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATggg                                                         <  1:2297188/100‑1 (MQ=255)                                       cGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGc                                                        >  1:1703113/1‑100 (MQ=255)                                                aCGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGc                                                >  1:1213227/1‑99 (MQ=255)                                                 cgcgTTTGAGTGGCTGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCg                                              <  1:63957/100‑1 (MQ=255)                                                 cgcgTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGTCTAAAtt                                                  >  1:1196976/1‑96 (MQ=255)                                                 cgcgTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCg                                              <  2:1749874/100‑1 (MQ=255)                                                 cgcgTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCg                                              <  1:193833/100‑1 (MQ=255)                                                    gTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGAt                                           >  1:386177/1‑100 (MQ=255)                                                      ttGAGTGTCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATgc                                         <  2:1670534/100‑1 (MQ=255)                                                        gAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCGTTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATgcgg                                       >  1:438699/1‑100 (MQ=255)                                                        gAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATgcgg                                       <  1:537012/100‑1 (MQ=255)                                                          gTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATgcggcg                                     >  2:1459669/1‑100 (MQ=255)                                                          gTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCAGATAATGGGCTAAATTGCAGGATgcggcg                                     >  2:226799/1‑100 (MQ=255)                                                              cAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAg                                 <  1:1837056/100‑1 (MQ=255)                                                              cAGGATCAAGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTACTTATAATGGTCTAAATTGCCGGATGCGGcgcg                                   <  2:1574572/98‑1 (MQ=255)                                                               tttattacgacgAATGTGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAg                                 <  1:2053896/93‑1 (MQ=255)                                                               aGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGt                                <  2:123625/100‑1 (MQ=255)                                                                  aTCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATCATGGGCTAAATTGCCGGATACGACGCAAGTAc                              >  1:296679/1‑99 (MQ=255)                                                                  aTCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACt                             <  2:160060/100‑1 (MQ=255)                                                                  aTCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACt                             >  2:38404/1‑100 (MQ=255)                                                                   tCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACtt                            <  2:1964557/100‑1 (MQ=255)                                                                    cACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACttt                           <  2:1199203/100‑1 (MQ=255)                                                                       gacgaAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATc                        <  1:209273/100‑1 (MQ=255)                                                                        acgaAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATcc                       <  2:586114/100‑1 (MQ=255)                                                                        acgaAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATcc                       <  1:564504/100‑1 (MQ=255)                                                                         cgaAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCg                      <  2:1344955/100‑1 (MQ=255)                                                                         cgaAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCg                      <  1:1006271/100‑1 (MQ=255)                                                                           aaGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGAt                    >  1:2224079/1‑100 (MQ=255)                                                                            aGGCGACAGTGATAGATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATc                   <  2:2122567/100‑1 (MQ=255)                                                                            aGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATc                   <  1:1855924/100‑1 (MQ=255)                                                                            aGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATc                   <  1:925931/100‑1 (MQ=255)                                                                               cGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCtat                >  2:964294/1‑100 (MQ=255)                                                                                   aGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAAt            <  1:242282/100‑1 (MQ=255)                                                                                      gATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGt          >  1:1132714/1‑99 (MQ=255)                                                                                      gATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAAATGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTa         >  2:1682734/1‑100 (MQ=255)                                                                                        tAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTAgg       >  2:357016/1‑100 (MQ=255)                                                                                         aaaTTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTAGGc      >  1:730210/1‑100 (MQ=255)                                                                                            ttCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTAGGCCg    >  2:312748/1‑99 (MQ=255)                                                                                             tCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTAATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGTTCTATAAATGTAGGCCGGa  <  2:556027/100‑1 (MQ=255)                                                                                                  |                                                                                               ATTTTACGGTTATGATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTAGGCCGGA  >  NC_000913/1382723‑1382914

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