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	92.8%	262.2 / 7.6	83	Y292H (TAT→CAT)	ompG	outer membrane porin G
Reads supporting (aligned to +/- strand):  ref base T (3/3);  new base C (42/35);  total (45/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

TTACGGTTATGATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTAGGCCGGA  >  NC_000913/1382726‑1382914                                                                                               |                                                                                               ttACGGTTATGATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCa                                                                                           <  1:325678/100‑1 (MQ=255)  tACGGTTATGATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAg                                                                                          <  1:1889955/100‑1 (MQ=255)  tACGGTTATGATTTCCAGAACGGACTTTCCGTTGCGCTGGAATACGCGTTTGAGTTGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAg                                                                                          <  1:1915068/100‑1 (MQ=255)      gTTATGATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGgtgt                                                                                      >  2:953053/1‑100 (MQ=255)      gTTATGATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGgtgt                                                                                      >  2:999227/1‑100 (MQ=255)        tATGATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCg                                                                                    >  1:37765/1‑100 (MQ=255)        tATGATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCg                                                                                    >  1:37376/1‑100 (MQ=255)           gATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCg                                                                                 >  1:208496/1‑100 (MQ=255)           gATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCg                                                                                 <  2:1679513/100‑1 (MQ=255)            aTTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGt                                                                                <  1:288559/100‑1 (MQ=255)            aTTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGt                                                                                >  2:142376/1‑100 (MQ=255)            aTTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGt                                                                                >  2:230153/1‑100 (MQ=255)             tttCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTa                                                                               <  1:895189/100‑1 (MQ=255)              ttCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTaa                                                                              >  2:1801809/1‑100 (MQ=255)              ttCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTaa                                                                              >  1:1557777/1‑100 (MQ=255)              ttCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTaa                                                                              <  1:1284229/100‑1 (MQ=255)                ccAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAAt                                                                            <  2:208497/100‑1 (MQ=255)                 cAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAAtt                                                                           <  1:686979/100‑1 (MQ=255)                   gAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTAc                                                                         >  2:1729676/1‑100 (MQ=255)                   gAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTAc                                                                         <  1:897556/100‑1 (MQ=255)                    aaCGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACt                                                                        <  2:337765/100‑1 (MQ=255)                     aCGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTc                                                                       >  2:1192953/1‑100 (MQ=255)                      cGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCg                                                                      <  1:1262109/100‑1 (MQ=255)                       ggACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGTTGTCGGCGTACATTACTCGt                                                                     >  1:768538/1‑100 (MQ=255)                          cTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCt                                                                  <  1:500454/100‑1 (MQ=255)                               cGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATaa                                                             >  1:1093090/1‑100 (MQ=255)                               cGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATaa                                                             >  2:834876/1‑100 (MQ=255)                               cGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATaa                                                             >  1:1193488/1‑100 (MQ=255)                                  ttCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATgg                                                          <  1:1617176/100‑1 (MQ=255)                                  ttCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATgg                                                          <  2:539405/100‑1 (MQ=255)                                  ttCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATgg                                                          <  1:552380/100‑1 (MQ=255)                                   tCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATggg                                                         >  2:833027/1‑100 (MQ=255)                                    cGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGc                                                        >  2:1171524/1‑100 (MQ=255)                                    cGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGc                                                        >  2:1339512/1‑100 (MQ=255)                                     gCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCt                                                       <  2:37376/100‑1 (MQ=255)                                     gCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCt                                                       <  2:37765/100‑1 (MQ=255)                                        ggAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTaaa                                                    >  1:742594/1‑100 (MQ=255)                                        ggAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATGACTCGTTCTGATAATGGGCTaaa                                                    >  2:476252/1‑100 (MQ=255)                                             aCGCGTTTGAGTGGCAGTATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGcc                                               <  1:1750125/100‑1 (MQ=255)                                              cgcgTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCg                                              <  2:1063747/100‑1 (MQ=255)                                                cgTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGa                                            <  1:863748/100‑1 (MQ=255)                                                 gTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGAt                                           >  1:1374793/1‑100 (MQ=255)                                                 gTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGAt                                           >  1:947454/1‑100 (MQ=255)                                                  tttGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATg                                          <  2:902361/100‑1 (MQ=255)                                                   ttGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATgc                                         <  2:1624509/100‑1 (MQ=255)                                                    tGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATgcg                                        >  1:1407102/1‑100 (MQ=255)                                                     gAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATgcgg                                       >  1:414449/1‑100 (MQ=255)                                                         ggCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGcgcg                                   >  2:1833609/1‑100 (MQ=255)                                                          gCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGa                                  >  2:1783230/1‑100 (MQ=255)                                                           cAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAg                                 >  2:1219325/1‑100 (MQ=255)                                                           cAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAg                                 >  1:943747/1‑100 (MQ=255)                                                             ggATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGt                                <  1:702973/99‑1 (MQ=255)                                                             ggATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTa                               >  2:247220/1‑100 (MQ=255)                                                              gATCACGACGAAGGCGACAGTGATAAATTCCAGCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTAc                              <  2:578931/100‑1 (MQ=255)                                                                tCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACtt                            <  1:632700/100‑1 (MQ=255)                                                                 cACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACttt                           >  2:439192/1‑100 (MQ=255)                                                                  acgacgAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTa                          <  1:1783219/100‑1 (MQ=255)                                                                  acgacgAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTa                          <  2:1148993/100‑1 (MQ=255)                                                                   cgacgaAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTAt                         >  1:1778212/1‑100 (MQ=255)                                                                   cgacgaAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTAt                         >  2:1491860/1‑100 (MQ=255)                                                                   cgacgaAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTAt                         >  2:289912/1‑100 (MQ=255)                                                                       gaAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGa                     >  1:971389/1‑100 (MQ=255)                                                                         aGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATc                   <  2:1072510/100‑1 (MQ=255)                                                                          ggCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCt                  >  1:660703/1‑100 (MQ=255)                                                                          ggCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCt                  <  1:1687754/100‑1 (MQ=255)                                                                               cAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATaaa             >  2:2748/1‑100 (MQ=255)                                                                               cAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATaaa             <  2:1407110/100‑1 (MQ=255)                                                                                aGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAAt            <  2:1433037/100‑1 (MQ=255)                                                                                 gTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATg           <  1:759427/100‑1 (MQ=255)                                                                                 gTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATg           <  2:723976/100‑1 (MQ=255)                                                                                  tGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGt          >  1:651349/1‑100 (MQ=255)                                                                                   gATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTa         >  2:1305029/1‑100 (MQ=255)                                                                                   gATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTa         <  2:691658/100‑1 (MQ=255)                                                                                   gATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTa         >  1:1885893/1‑100 (MQ=255)                                                                                    aTAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAGTGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTAg        >  1:1085082/1‑100 (MQ=255)                                                                                    aTAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTAg        <  1:1833597/100‑1 (MQ=255)                                                                                    aTAAATTCCATCACGCAGGTGTCGGCGTAAATTACTCGATCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTAg        >  1:1198844/1‑100 (MQ=255)                                                                                      aaaTTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTAGGc      >  2:1119609/1‑100 (MQ=255)                                                                                      aaaTTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGAGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTAGGc      >  2:466769/1‑100 (MQ=255)                                                                                        aTTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTAGGCCg    >  2:888010/1‑100 (MQ=255)                                                                                         ttCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTAGGCCgg   <  1:833026/100‑1 (MQ=255)                                                                                          tCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTAGGCCGGa  <  2:1093091/100‑1 (MQ=255)                                                                                          tCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTAGGCCGGa  >  2:1203461/1‑100 (MQ=255)                                                                                               |                                                                                               TTACGGTTATGATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTAGGCCGGA  >  NC_000913/1382726‑1382914

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