Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,382,820	T→C	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	98.2%	183.3 / ‑4.1	55	Y292H (TAT→CAT)	ompG	outer membrane porin G
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (28/25);  total (30/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.35e-01

ACGGTTATGATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTAGGCCGGATA  >  NC_000913/1382728‑1382916                                                                                             |                                                                                                 aCGGTTATGATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAgg                                                                                           <  2:280773/100‑1 (MQ=255)   ggTTATGATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGgtg                                                                                         >  1:1176345/1‑100 (MQ=255)     ttATGATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTc                                                                                       >  1:1468547/1‑100 (MQ=255)         gATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCg                                                                                   >  2:1192699/1‑100 (MQ=255)         gATTTCCAGAACGGACTTTCCGTTTCGCTGAAATACGCGTTTGAGTGGCAGGATGACGACGAAGGCGACAGTGATTAATTCCATCATGCAGGTGTCGGCg                                                                                   <  2:1377132/100‑1 (MQ=255)           tttCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTa                                                                                 >  2:138499/1‑100 (MQ=255)           tttCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTa                                                                                 >  2:124052/1‑100 (MQ=255)             tCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTaaa                                                                               >  1:783044/1‑100 (MQ=255)              ccAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAAt                                                                              <  1:1050971/100‑1 (MQ=255)              ccAGAACGGACTTTCCGTTTCGCGGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAAt                                                                              <  2:1229191/100‑1 (MQ=255)               cAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAAtt                                                                             >  1:813889/1‑100 (MQ=255)               cAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAAtt                                                                             >  2:19412/1‑100 (MQ=255)               cAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAAtt                                                                             <  2:1485292/100‑1 (MQ=255)               cAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGCGATAAATTCCATCATGCAGGTGTCGGCGTAAAtt                                                                             >  1:325769/1‑100 (MQ=255)                    cGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCg                                                                        <  2:275416/100‑1 (MQ=255)                     ggACTTTCCGTTTCGCTGGAATACGCGTTTTAGTGGCAGGATCACGACGAAGTCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCg                                                                        <  2:538650/99‑1 (MQ=255)                      gACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGtt                                                                      >  1:1064327/1‑100 (MQ=255)                         tttCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTg                                                                   <  2:856087/100‑1 (MQ=255)                         tttCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTg                                                                   <  2:1170060/100‑1 (MQ=255)                          ttCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTg                                                                   <  2:1542574/99‑1 (MQ=255)                          ttCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGa                                                                  >  1:497796/1‑100 (MQ=255)                           tCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGAt                                                                 <  2:1242689/100‑1 (MQ=255)                           tCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGAt                                                                 <  2:256331/100‑1 (MQ=255)                            ccGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGAt                                                                 <  2:506059/99‑1 (MQ=255)                                ttCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATgg                                                            <  2:1552289/100‑1 (MQ=255)                                  cGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGCCCGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGc                                                          >  2:418959/1‑100 (MQ=255)                                   gCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATGATGCAGGTGTCGGCGTAAAGTACTCGTTCTGCTAATGGACt                                                         >  2:1533354/1‑100 (MQ=255)                                   gCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCt                                                         <  1:417376/100‑1 (MQ=255)                                    cTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTACATCATACAGGTGTCGGCGTAAATTACTCGTTCTGCGAATGGGCTa                                                        >  1:686527/1‑100 (MQ=255)                                        aaTACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAAtt                                                    >  1:992868/1‑100 (MQ=255)                                              cgTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGa                                              >  1:697872/1‑100 (MQ=255)                                               gTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGAt                                             >  2:554191/1‑100 (MQ=255)                                                tttGAGTGGCAGTTTCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCg                                                <  1:59320/96‑1 (MQ=255)                                                 ttGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATgc                                           >  1:655996/1‑100 (MQ=255)                                                   gAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATgcgg                                         <  2:783040/100‑1 (MQ=255)                                                     gTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATgcggcg                                       <  2:1312637/100‑1 (MQ=255)                                                     gTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATgcggcg                                       >  1:1196549/1‑100 (MQ=255)                                                     gTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATgcggcg                                       >  1:385645/1‑100 (MQ=255)                                                       ggCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGcgcg                                     >  1:1133063/1‑100 (MQ=255)                                                        gCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGa                                    >  1:406893/1‑100 (MQ=255)                                                         cAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGan                                   <  2:29159/100‑2 (MQ=255)                                                              tCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACtt                              <  2:242079/100‑1 (MQ=255)                                                                acgacgAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTa                            >  1:1440457/1‑100 (MQ=255)                                                                   acgaAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATcc                         <  2:1468539/100‑1 (MQ=255)                                                                    caaaGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCg                        >  1:1073304/3‑100 (MQ=255)                                                                       aGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATc                     <  2:300861/100‑1 (MQ=255)                                                                        ggCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCt                    <  2:287566/100‑1 (MQ=255)                                                                        ggCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCt                    <  1:124052/100‑1 (MQ=255)                                                                        ggCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCt                    >  2:1228648/1‑100 (MQ=255)                                                                        ggCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCt                    <  2:1176340/100‑1 (MQ=255)                                                                        ggCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGGACTTTATCCGATCt                    >  2:80052/1‑100 (MQ=255)                                                                               gTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATg             >  1:470152/1‑100 (MQ=255)                                                                                    aaaTTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCTATCTATAAATGTAGGc        >  2:463354/1‑100 (MQ=255)                                                                                        tCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATACGAAATATAAATGTAGACCGGa    >  2:174381/1‑100 (MQ=255)                                                                                         ccATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTAGGCCGGAt   >  2:1139752/1‑100 (MQ=255)                                                                                          catcatGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTAGGCCGGATa  >  2:289455/1‑100 (MQ=255)                                                                                          catcatGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTAGGCCGGATa  >  2:286809/1‑100 (MQ=255)                                                                                          catcatGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTAGGCCGGATa  <  2:614697/100‑1 (MQ=255)                                                                                             |                                                                                                 ACGGTTATGATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTAGGCCGGATA  >  NC_000913/1382728‑1382916

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: ‑