Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	407,573	A→G	30.3%	Y49C (TAC→TGC)	aroM →	AroM family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	407,573	0	A	G	30.3%	67.2 / 37.6	56	Y49C (TAC→TGC)	aroM	AroM family protein
Reads supporting (aligned to +/- strand):  ref base A (18/21);  new base G (6/11);  total (24/32)
Fisher's exact test for biased strand distribution p-value = 5.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAGTTTTGCCGCTCCTGACGGAATACATTGACGAAGATAATATTTCCCATCATAGCCTGCTGGGGAAGTTAAGTCGTGAAGAAGTGATGGCGGAGTACGCGCCAGAAGCAGGCGAAGACACCATTCTCACATTATTAAATGACAACCAGCTGGCCCATGTTTCGCGTCGCAAAGTGGAGCGTGACCTGCAAG  >  NC_000913/407476‑407668                                                                                                  |                                                                                                gAAGTTTTGCCGCTCCTGACGGAATACATTGACGAAGATAATATTTCCCATCATAGCCTGCTGGGGAAGTTAAGTCGTGAAGAAGTGATGGCGGAGTAcg                                                                                               <  1:158894/100‑1 (MQ=255)     ttttGCCGCTCCTGAGGGAATACATTGACGAAGATAATATGTCCCATCATAGCCTGCTGGGGAAGTTAAGTCGTGAAGAAGTGATGGCGGAGTACGCGcc                                                                                           <  1:1879857/100‑1 (MQ=255)     ttttGCCGCTCCTGACGGAATACATTGACGAAGATAATATTTCCCATCATAGCCTGCTGGGGAAGTTAAGTCGTGAAGAAGTGATGGCGGAGTGCGCGcc                                                                                           <  1:426305/100‑1 (MQ=255)       ttGCCGCTCCTGACGGAATACATTGACGAAGATAATATTTCCCATCATAGCCTGCTGGGGAAGTTAAGTCGTGAAGAAGTGATGGCGGAGTGCGCGCCAg                                                                                         <  2:1958841/100‑1 (MQ=255)        tGCCGCTCCTGACGGAATACATTGACGAAGATAATATTTCCCATCATAGCCTGCTGGGGAAGTTAAGTCGTGAAGAAGTGATGGCGGAGTACGCGCCAGa                                                                                        <  1:125222/100‑1 (MQ=255)           cGCTCCTGACGGAATACATTGACGAAGATAATATTTCCCATCATAGCCTGCTGGGGAAGTTAAGTCGTGAAGAAGTGATGGCGGAGTGCGCGCCAGAAg                                                                                      <  1:761786/99‑1 (MQ=255)           cGCTCCTGACGGAATACATTGACGAAGATAATATTTCCCATCATAGCCTGCTGGGGAAGTTAAGTCGTGAAGAAGTGATGGCGGAGTACGCGCCAGAAGc                                                                                     <  1:1733878/100‑1 (MQ=255)           cGCTCCTGACGGAATACATTGACGAAGATAATATTTCCCATCATAGCCTGCTGGGGAAGTTAAGTCGTGAAGAAGTGATGGCGGAGTACGCGCCAGAAGc                                                                                     <  1:2246723/100‑1 (MQ=255)            gCTCCTGACGGAATACATTGACGAAGATAATATTTCCCATCATAGCCTGCTGGGGAAGTTAAGTCGTGAAGAAGTGATGGCGGAGTACGCGCCAGAAGCa                                                                                    <  2:908015/100‑1 (MQ=255)            gCTCCTGACGGAATACATTGACGAAGATAATATTTCCCATCATAGCCTGCTGGGGAAGTTAAGTCGTGAAGAAGTGATGGCGGAGTACGCGCCAGAAGCa                                                                                    <  1:927386/100‑1 (MQ=255)                 tGACGGAATACATTGACGAAGATAATATTTCCCATCATAGCCTGCTGGGGAAGTTAAGTCGTGAAGAAGTGATGGCGGAGTGCGCGCCAGAAGCAGGCGa                                                                               >  1:2063190/1‑100 (MQ=255)                 tGACGGAATACATTGACGAAGATAATATTTCCCATCATAGCCTGCTGGGGAAGTTAAGTCGTGAAGAAGTGATGGCGGAGTACGCGCCAGAAGCAGGCGa                                                                               <  1:497172/100‑1 (MQ=255)                  gACGGAATACATTGACGAAGATAATATTTCCCATCATAGCCTGCTGGGGAAGTTAAGTCGTGAAGAAGTGATGGCGGAGTGCGCGCCAGAAGCAGGCGaa                                                                              <  2:1138096/100‑1 (MQ=255)                  gACGGAATACATTGACGAAGATAATATTTCCCATCATAGCCTGCTGGGGAAGTTAAGTCGTGAAGAAGTGATGGCGGAGTACGCGCCAGAAGCAGGCGaa                                                                              >  1:1882859/1‑100 (MQ=255)                    cGGTATACATTGACGAAGATAATATTGTGGAACATAGTCTGCTGGGGAAGTTAAGTGGGGAAGAAGTGATGGCTGAGTACGCGCCAGAAGCAgg                                                                                  <  2:1083680/94‑1 (MQ=255)                    cGGAATACATTGACGAAGATAATATTTCCCATCATAGCCTGCTGGGGAAGTTAAGTCGTGAAGAAGTGATGGCGGAGTGCGCGCCAGAAGCAGGCGAAGa                                                                            <  1:135592/100‑1 (MQ=255)                        aTACATTGACGAAGATAATATTTCCCATCATAGCCTGCTGGGGAAGTTAAGTCGTGAAGAAGTGATGGCGGAGTGCGCGCCAGAAGCAGGCGAAGACAcc                                                                        <  2:699542/100‑1 (MQ=255)                        aTACATTGACGAAGATAATATTTCCCATCATAGCCTACTGGGGAAGTTAAGTCGTGAAGAAGTGATGGCGGAGTAGGCGACAGAAGCAGGCGAAGACAcc                                                                        >  2:541750/1‑100 (MQ=255)                         tACATTGACGAAGATAATATTTCCCATCATAGCCTGCTGGGGAAGTTAAGTCGTGAAGAAGTGATGGCGGAGTGCGCGCCAGAAGCAGGCGAAGACACCa                                                                       >  2:454388/1‑100 (MQ=255)                               gACGAAGATAATATTTCCCATCATAGCCTGCTGGGGAAGTTAAGTCGTGAAGAAGTGATGGCGGAGTACGCGCCAGAAGCAGGCGAAGACACCATTCTca                                                                 >  2:798998/1‑100 (MQ=255)                                 cGAAGATAATATTTCCCATCATAGCCTGCTGGGGAAGTTAAGTCGTGAAGAAGTGATGGCGGAGTGCGCGCCAGAAGCAGGCGAAGACACCATTCTcaca                                                               <  2:457916/100‑1 (MQ=255)                                  gAAGATAATATTTCCCATCATAGCCTGCTGGGGAAGTTAAGTCGTGAAGAAGTGATGGCGGAGTACGCGCCAGAAGCAGGCGAAGACACCATTCTCACat                                                              <  1:1100638/100‑1 (MQ=255)                                   aaGATAATATTTCCCATCATAGCCTGCTGGGGAAGTTAAGTCGTGAAGAAGTGATGGCGGAGTGCGCGCCAGAAGCAGGCGAAGACACCATTCTCACatt                                                             <  2:1372531/100‑1 (MQ=255)                                   aaGATAATATTTACCATCATAGCCTGCTGGGGAAGTTAAGTCGTGAAGAAGTGATGGCGGAGTACGCGCCAGAAGCAGGCGAAGACACCATTCTCACatt                                                             <  2:527976/100‑1 (MQ=255)                                       taataTTTCCCATCATAGCCTGCTGGGGAAGTTAAGTCGTGAAGAAGTGATGGCGGAGTACGCGCCAGAAGCAGGCGAAGACACCATTCTCACAttatta                                                         >  2:2248061/1‑100 (MQ=255)                                        aataTTTCCCATCATAGCCTGCTGGGGAAGTTAAGTCGTGAAGAAGTGATGGCGGAGTACGCGCCAGAAGCAGGCGAAGACACCATTCTCACATTATTaa                                                        <  1:928929/100‑1 (MQ=255)                                          tatTTCCCATCATAGCCTGCTGGGGAAGTTAAGTCGTGAAGAAGTGATGGCGGAGTGCGCGCCAGAAGCAGGCGAAGACACCATTCTCACATTATTAAAt                                                      <  1:2042702/100‑1 (MQ=255)                                          tatTTCCCATCATAGCCTGCTGGGGAAGTTAAGTCGTGAAGAAGTGATGGCGGAGTACGCGCCAGAAGCAGGCGAAGACACCATTCTCACATTATTAAAt                                                      <  2:316360/100‑1 (MQ=255)                                           atTTCCCATCATAGCCTGCTGGGGAAGTTAAGTCGTGAAGAAGTGATGGCGGAGTACGCGCCAGAAGCAGGCGAAGACACCATTCTCACATTATTAAATg                                                     >  1:1371668/1‑100 (MQ=255)                                           atTTCCCAGCATAGCCTGCTGGGGAAGTTCAGTCGTGAAGAAGTGATGGCGGAGTACGCGCCAGAAACAGGCGAAGACACCATTCTCACATTATTAAATg                                                     <  2:1349435/100‑1 (MQ=255)                                            tttCCCATCATAGCCTGCTGGGGAAGTTAAGTCGTGAAGAAGTGATGGCGGAGTGCGCGCCAGAAGCAGGCGAAGACACCATTCTCACATTATTAAATGa                                                    <  2:639635/100‑1 (MQ=255)                                               cccATCATAGCCTGCTGGGGAAGTTAAGTCGTGAAGAAGTGATGGCGGAGTGCGCGCCAGAAGCAGGCGAAGACACCATTCTCACATTATTAAATGACaa                                                 <  2:567481/100‑1 (MQ=255)                                               cccATCATAGCCTGCTGGGGAAGTTAAGTCGTGAAGAAGTGATGGCGGAGTGCGCGCCAGAAGCAGGCGAAGACACCATTCTCACATTATTAAATGACaa                                                 >  1:2222592/1‑100 (MQ=255)                                                 catcatAGCCTGCTGGGGAAGTTAAGTCGTGAAGAAGTGATGGCGGAGTACGCGCCAGAAGCAGGCGAAGACACCATTCTCACATTATTAAATGACAAcc                                               >  2:911219/1‑100 (MQ=255)                                                  atcatAGCCTGCTGGGGAAGTTAAGTCGTGAAGAAGTGATGGCGGAGTACGCGCCAGAAGCAGGCGAAGACACCATTCTCACATTATTAAATGACAACCa                                              <  1:651511/100‑1 (MQ=255)                                                       aGCCTGCTGGGGACGTTAAGTCGTGAAGAAGTGATGGCGGAGTACGCGCCAGAAGCAGGCGAAGACACCATTCTCACATTATTAAATGACACCCAGc                                            >  1:743158/1‑97 (MQ=255)                                                        gCCTGCTGGGGAAGTTAAGTCGTGAAGAAGTGATGGCGGAGTACGCGCCAGAAGCAGGCGAAGACACCATTCTCACATTATTAAATGACAACCAGCTGGc                                        <  2:1992352/100‑1 (MQ=255)                                                          ctgctgGGGAAGTTAAGTCGTGAAGAAGTGATGGCGGAGTACGCGCCAGAAGCAGGCGAAGACACCATTCTCACATTATTAAATGACAACCAGCTGGc                                        >  1:97751/1‑98 (MQ=255)                                                              tgGGGAAGTTAAGTCGTGAAGAAGTGATGGCGGAGTACGCGCCAGAAGCAGGCGAAGACACCATTCTCACATTATTAAATGACAACCAGCTGGCCCATGt                                  >  2:920074/1‑100 (MQ=255)                                                                 ggAAGTTAAGTCGTGAAGAAGTGATGGCGGAGTACGCGCCAGAAGCAGGCGAAGACACCATTCTCACATTATTAAATGACAACCAGCTGGCCCATGtt                                 >  2:2176167/1‑98 (MQ=255)                                                                  gAAGTTAAGTCGTGAAGAAGTGATGGCGGAGTACGCGCCAGAAGCAGGCGAAGACACCATTCTCACATTATTAAATGACAACCAGCTGGCCCATGTTTcg                              >  2:832262/1‑100 (MQ=255)                                                                     gTTAAGTCGTGAAGAAGTGATGGCTGAGTACGCGCCAGAAGCAGGCGAAGACACCATTCTCACATTATTAAATGACAACCAGCTGGCCCATGTTTagc                             >  1:875818/1‑98 (MQ=255)                                                                     gTTAAGTCGTGAAGAAGTGATGGCGGAGTACGCGCCAGAAGCAGGCGAAGACACCATTCTCACATTATTAAATGACAACCAGCTGGCCCATGTTTCGCGt                           >  2:636702/1‑100 (MQ=255)                                                                       tAAGTCGTGAAGAAGTGATGGCGGAGTACGCGCCAGAAGCAGGCGAAGACACCATTCTCACATTATTAAATGACAACCAGCTGGCCCATGTTTCGCGTCg                         >  2:2051769/1‑100 (MQ=255)                                                                       tAAGTCGTGAAGAAGTGATGGCGGAGTACGCGCCAGAAGCAGGCGAAGACACCATTCTCACATTATTAAATGACAACCAGCTGGCCCATGTTTCGCGTCg                         <  1:1594382/100‑1 (MQ=255)                                                                          gTCGTGAAGAAGTGATGGCGGAGTGCGCGCCAGAAGCAGGCGAACACACCACTCTCACGTTATTAAATGACAACCCGGTGGCCCAAGTTTCGCGTCGCa                       >  2:2138936/1‑99 (MQ=255)                                                                          gTCGTGAAGAAGTGATGGCGGAGTACGCGCCAGAAGCAGGCGAAGACACCATTCTCACATTATTAAATGACAACCAGCTGGCCCATGTTTCGCGTCGCaa                      >  1:1122898/1‑100 (MQ=255)                                                                           tCGTGAAGAAGTGATGGCGGAGTACGCGCAAGAAGCAGGCGCAGAAACCATTCTCACATTATTAAATGACAACCAGTTTGCCCATGTTTCGCGTCGCaaa                     <  2:742886/100‑1 (MQ=255)                                                                            cGTGAAGAAGTGATGGCGGAGTACGCGCCAGAAGCAGGCGAAGACACCATTCTCACATTATTAAATGACAACAAGCTGGCCCATGTTTCGCGTCGCaaa                     >  2:1440578/1‑99 (MQ=255)                                                                             gTGAAGAAGTGATCGCGGAGTGCGCGCCAGAAGCAGGCGAAGACACCATTCTCACATTATTAAATGACAACCAGCTGGCCCATGTTTCGCGTCGCAAAGt                   >  1:1506489/1‑100 (MQ=255)                                                                              tGAAGAAGTGATGGCGGAATACGCGCCAGAAGCAGGCGAAGACACCATTCTCACATTATTAAATGACAACCAGCTGGCCCATGTTTCGCGTCGCAAAGTg                  <  1:572657/100‑1 (MQ=255)                                                                                   aagTGATGGCGGAGTACGCGCCAGAAGCAGGCGAAGACACCATTCTCACATTATTAAATGACAACCAGCTGGCCCATGTTTCGCGTCGCAAAGTGGAGCg             >  1:1001693/1‑100 (MQ=255)                                                                                    agTGATGGCGGAGTGCGCGCCAGAAGCAGGCGAAGACACCATTCTCACATTATTAAATGACAACCAGCTGGCCCATGTTTCGCGTCGCAAAGTGGAGCGt            >  1:1373770/1‑100 (MQ=255)                                                                                        aTGGCGGAGTACGCGCCAGAAGCAGGCGAAGACACCATTCTCACATTATTAAATGACAACCAGCTGGCCCATGTTTCGCGTCGCAAAGTGGAGCGTGAcc        <  1:1939701/100‑1 (MQ=255)                                                                                           gCGGAGTACGCGCCAGAAGCAGGCGAAGACACCATTCTCACATTATTAAATGACAACCAGCTGGCCCATGTTTCGCGTCGCAAAGTGGAGCGTGACCt       >  2:612132/1‑98 (MQ=255)                                                                                              gAGTACGCGCCAGAAGCAGGCGAAGACACCATTCTCACATTATTAAATGACAACCAGCTGGCCCATGTTTCGCGTCGCAAAGTGGAGCGTGACCTGCAAg  <  2:1622952/100‑1 (MQ=255)                                                                                                  |                                                                                                GAAGTTTTGCCGCTCCTGACGGAATACATTGACGAAGATAATATTTCCCATCATAGCCTGCTGGGGAAGTTAAGTCGTGAAGAAGTGATGGCGGAGTACGCGCCAGAAGCAGGCGAAGACACCATTCTCACATTATTAAATGACAACCAGCTGGCCCATGTTTCGCGTCGCAAAGTGGAGCGTGACCTGCAAG  >  NC_000913/407476‑407668

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