Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	1,047,955	G→A	8.2%	D324N (GAT→AAT)	menF →	isochorismate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	1,047,955	0	G	A	8.2%	142.8 / 4.2	49	D324N (GAT→AAT)	menF	isochorismate synthase
Reads supporting (aligned to +/- strand):  ref base G (20/25);  new base A (2/2);  total (22/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.96e-01

AACGTTCACAAGATGAGGACGAAGATACAAAAAATGTTGAAGCATTTTTAAAAGATAATAAAAACTTAATCGAACATCGATTTGTTGTTGACAGTATTTTACATGATATTAAACCTTATATCACTGAATTACATTATGATAAGACGCCTAAAATTCTAAAAAATGATCATTTATATCACTTGTACACTGAAATAAAGGCGCCACTGAAGGATGATTCGTATATTAGTTTAATTGATAATTTACATCCAACACCTGCTTTAGGTGGCTATCCA  >  CP000730/1047818‑1048089                                                                                                                                          |                                                                                                                                       aaCGTTCACAAGATGAGGACGAAGATACAAAAAATGTTGAAGCATTTTTAAAAGATAATAAAAACTTAATCGAACATCGATTTGTTGTTGACAGTATTTTACATGATATTAAACCTTATATCACTGAATTACATTATGATa                                                                                                                                     >  2:394337/1‑141 (MQ=255)     ttCACAAGATGAGGACGAAGATACAAAAAATGTTGAAGCATTTTTAAAAGATAATAAAAACTTAATCGAACATCGATTTGTTGTTGACAGTATTTTACATGATATTAAACCTTATATCACTGAATTACATTATGATAAGAc                                                                                                                                 >  1:230960/1‑141 (MQ=255)       cacaAGATGAGGACGAAGATACAAAAAATGTTGAAGCATTTTTAAAAGATAATAAAAACTTAATCGAACATCGATTTGTTGTTGACAGTATTTTACATGATATTAAACCTTATATCACTGAATTACATTATGATAAGACGc                                                                                                                               >  2:56299/1‑141 (MQ=255)       cacaAGATGAGGACGAAGATACAAAAAATGTTGAAGCATTTTTAAAAGATAATAAAAACTTAATCGAACATCGATTTGTTGTTGACAGTATTTTACATGATATTAAACCTTATATCACTGAATTACATTATGATAAGACGc                                                                                                                               >  2:412012/1‑141 (MQ=255)             aTGAGGACGAAGATACAAAAAATGTTGAAGCATTTTTAAAAGATAATAAAAACTTAATCGAACATCGATTTGTTGTTGACAGTATTTTACATGATATTAAACCTTATATCACTGAATTACATTATGATAAGACGCCTaaaa                                                                                                                         >  2:502100/1‑141 (MQ=255)               gAGGTCGAAGATACAAAAAATGTTGAAGCATTTTTAAAAGATAATAAAAACTTAATCGAACATCGATTTGTTGTTGACAGTATTTTACATGATATTAAACCTTATATCACTGAATTACATTATGATAAGACGCCTAAAAtt                                                                                                                       <  2:443434/141‑1 (MQ=255)                 ggACGAAGATACAAACAATGTAGAAGCAGTTTTAAAAGATAATAAAAACTTAATCGAATAGCGATTTGTTGTTGACAGTATTTTGCATGATATAAAACAATATATTACTTAATGACATTATGTGAAGCGGCCTGAAATTCt                                                                                                                     >  2:94316/1‑141 (MQ=255)                  gACGAAGATACAAAAAATGTTGAAGCATTTTTAAAAGATAATAAAAACTTAATCGAACATCGATTTGTTGTTGACAGTATTTTACATGATATTAAACCTTATATCACTGAATTACATTATAATAAGACGCCTAAAATTCTa                                                                                                                    >  2:169788/1‑141 (MQ=255)                  gACGAAGATACAAAAAATGTTGAAGCATTTTTAAAAGATAATAAAAACTTAATCGAACATCGATTTGTTGTTGACAGTATTTTACATGATATTAAACCTTATATCACTGAATTACATTATAATAAGACGCCTAAAATTCTa                                                                                                                    >  2:169787/1‑141 (MQ=255)                    cGAAGATACAAAAAATGTTGAAGCATTTTTAAAAGATAATAAAAACTTAATCGAACACCGATTTGTTGTTGACAGTATTTTACATGATATTAAACCTTATATCACTGAATTACATTATGATAAGACGCCTAAAATTCTaaa                                                                                                                  <  1:189316/141‑1 (MQ=255)                       aGATACAAAAAATGTTGAAGCATTTTTAAAAGATAATAAAAACTTAATCGAACATCGATTTGTTGTTGACAGTATTTTACATGATATTAAACCTTATATCACTGAATTACATTATGATAAGACGCCTAAAATTCTaaaaaa                                                                                                               <  2:64576/141‑1 (MQ=255)                        gATACAAAAAATGTTGAAGCATTTTTAAAAGATAATAAAAACTTAATCGAACATCGATTTGTTGTTGACAGTATTTTACATGATATTAAACCTTATATCACTGAATTACATTATGATAAGACGCCTAAAATTCTAAAAAAt                                                                                                              >  2:406060/1‑141 (MQ=255)                                   tGTTGAAGCATTTTTAAAAGATAATAAAAACTTAATCGAACATCGATTTGTTGTTGACAGTATTTTACATGATATTAAACCTTATATCACTGAATTACATTATGATAAGACGCCTAAAATTCTAAAAAATGATCATTtata                                                                                                   <  1:485517/141‑1 (MQ=255)                                    gTTGAAGCATTTTTAAAAGATAATAAAAACTTAATCGAACATCGATTTTTTGTTTACAGTATTTTACATGATATTAAACCTTATATCACTGAATTACATTATGATAATACGCCTAAAATTCTAAAAAATGATCATTTatat                                                                                                  >  2:357611/1‑141 (MQ=255)                                      tGAAGCATTTTTAAAAGATAATAAAAACTTAATCGAACATCGATTTGTTGTAGAAAGTATTTTACATGATATTAAACCTTATATCACTGAATTACATTATGATAAGACGACTAAAATTCTAAAAAATGATCATTTATATCa                                                                                                >  2:8405/1‑141 (MQ=255)                                          gCATTTTTAAAAGATAATAAAAACTTAATCGAACATCGATTTGTTGTTGACAGTATTTTACATGATATTAAACCTTATATCACTGAATTACATTATGATAAGACGCCTAAAATTCTAAAAAATGGTCATTTATATCACTTg                                                                                            >  1:406587/1‑141 (MQ=255)                                               tttAAAAGATAATAAAAACTTAATCGAACATCGATTTGTTGTTGACAGTATTTTACATGATATTAAACCTTATATCACTGAATTACATTATGATAAGACGCCTAAAATTCTAAAAAATGATCATTTATATCACTTGTacac                                                                                       <  1:138992/141‑1 (MQ=255)                                                     aGATAATAAAAACTTAATCGAACATCGATTTGTTGTTGACAGTATTTTACATGATATTAAACCTTATATCACTGAATTACATTATGATAAGACGCCTAAAATTCTAAAAAATGATCATTTATATAACTTGTAAAATGAAAt                                                                                 >  2:425794/1‑141 (MQ=255)                                                           taaAAACTTAATCGAACATCGATTTGTTGTTGACAGTATTTTACATGATATTAAACCTTATATCACTGAATTACATTATGATAAGACGCCTAAAATTCTAAAAAATGATCATTTATATCACTTGTACACTGAAATAAAGgc                                                                           <  1:406060/141‑1 (MQ=255)                                                           taaAAACTTAATCGAACATCGATTTGTAGTTGACAGTATTTTACATGATATTAAACCTTATATCACTGAATTACATTATGATAAGACGCCTAAAATTCTAAAAAATGATCATTTATATCACTTGTACACTGAAATAAAGgc                                                                           >  2:398860/1‑141 (MQ=255)                                                                   tAATCGAACATCGATTTGTTGTTGACAGTATTTTACATGATATTAAACCTTATATCACTGAATTACATTATGATAAGACGCCTAAAATTCTAAAAAATGATCATTTATATCACTTGTACACTGAAATAAAGGCGCCACTGa                                                                   <  1:218111/141‑1 (MQ=255)                                                                    aaTCGAACATCGATTTGTTGTTGACAGTATTTTACATGATATTAAACCTTATATCACTGAATTACATTATGATAAGACGCCTAAAATTCTAAAAAATGATCATTTATATCACTTGTACACTGAAATAAAGGCGCCACTGaa                                                                  >  1:93128/1‑141 (MQ=255)                                                                              cGATTTGTTGTTGACAGCCTGGGACATGATATTAAACCTTATATCACTGAATTACATTATGATAAGGCGCCTAAAATTCTAAAAAATGATCATTTATATCACTTGTACACTGAAATAAAGGCGCCACTGAAGGATGATTCg                                                        <  1:4012/141‑1 (MQ=255)                                                                                aTTTGTTGTTGACAGTATTTTACATGATATTAAACCTTATATCACTGAATTACATTATGATAAGACGCCTAAAATTCTAAAAAATGATCATTTATATCACTTGTACACTGAAATAAAGGCGCCACTGAAGGATGATTCGta                                                      <  1:464960/141‑1 (MQ=255)                                                                                 tttGTTGTTGACAGTATTTTACATGATATTAAACCTTATATCACTGAATTACATTATGATAAGACGCCTAAAATTCTAAAAAATGATCATTTATATCACTTGTACACTGAAATAAAGGCGCCACTGAAGGATGATTCGtat                                                     >  1:452853/1‑141 (MQ=255)                                                                                   tgttgttgACAGTATTTTACATGATATTAAACCTTATATCACTGAATTACATTATGATAAGACGCCTAAAATTCTAAAAAATGATCATTTATATCACTTGTACACTGAAATAAAGGCGCCACTGAAGGATGATTCGTatat                                                   <  2:34368/141‑1 (MQ=255)                                                                                    gttgttgACAGTATTTTACATGATATTAAACCTTATATCACTGAATTACATTATGATAAGACGCCTAAAATTCTAAAAAATGATCATTTATATCACTTGTACACTGAAATAAAGGCGCCACTGAAGGATGATTCGTATAtt                                                  <  1:55414/141‑1 (MQ=255)                                                                                        ttgACAGTATTTTACATGATATTAAACCTTATATCACTGAAATACATTATGATAAGACGCCTAAACTTCTAAAAAATGATCATTTATATCACTTGTACACTGAAATAAAGGCGCCACTGAAGGATTAATCGTATATTAGtt                                              <  2:314602/141‑1 (MQ=255)                                                                                            cAGTATTTTACATGATATTAAACCTTATATCACTGAATTACATTATGATAAGACGCCTAAAATAATAAAAAATGATCATTTATATCACTTGTACACTGAAATAAAGGCGCCACTGAAGGATGATTCGTATATTAGTTTAAt                                          <  2:150368/141‑1 (MQ=255)                                                                                             aGTATTTTACATGATATTAAACCTTATATCACTGAATTACATTATGATAAGACGCCTAAAATTCTAAAAAATGATCATTTATATCACTTGTACACTGAAATAAAGGCGCCACTGAAGGATGATTCGTATATTAGTTTAAtt                                         <  1:469913/141‑1 (MQ=255)                                                                                               tATTTTACATGATATTAAACCTTATATCACTGAATTACATTATGATAAGACGCCTAAAATTCTAAAAAATGATCATTTATATCACTTGTACACTGAAATAAAGGCGCCACTGAAGGATGATTCGTATATTAGTTTAATTGa                                       <  1:252557/141‑1 (MQ=255)                                                                                               tATTTTACATGATATTAAACCTTATATCACTGAATTACATTATAATAAGACGCCTAAAATTCTAAAAAATGATCATTTATATCACTTGTACACTGAAATAAAGGCGCCACTGAAGGATGATTCGTATATTAGTTTAATTGa                                       <  1:169787/141‑1 (MQ=255)                                                                                               tATTTTACATGATATTAAACCTTATATCACTGAATTACATTATAATAAGACGCCTAAAATTCTAAAAAATGATCATTTATATCACTTGTACACTGAAATAAAGGCGCCACTGAAGGATGATTCGTATATTAGTTTAATTGa                                       <  1:169788/141‑1 (MQ=255)                                                                                                 ttttACATGATATTAAACCTTATATCACTGAATTACATTATGATAAGACGCCTAAAATTCTAAAAAATGATCATTTATATCACTTGTACACTGAAATAAAGGCGCCACTGAAGGATGATTCGTATATTAGTTTAATTGATa                                     >  2:337726/1‑141 (MQ=255)                                                                                                  tttACATGATATTAAACCTTATATCACTGAATTACATTATGATAAGACGCCTAAAATTCTAAAAAATGATCATTTATATCACTTGTACACTGAAATAAAGGCGCCACTGAAGGATGATTCGTATATTAGTTTAATTGATaa                                    <  1:265928/141‑1 (MQ=255)                                                                                                   ttACATGATATTAAACCTTATATCACTGAATTACATTATGATAAGACGCCTAAAATTCTAAAAAATGATCATTTATATCACTTGTACACTGAAATAAAGGCGCCACTGAAGGATGATTCGTATATTAGTTTAATTGATAAt                                   >  2:308123/1‑141 (MQ=255)                                                                                                        tGATATTAAACCTTATATCACTGAATTACATTATGATAAGACGCCTAAAATTCTAAAAAATGATCATTTATATCACTTGTACACTGAAATAAAGGCGCCACTGAAGGATGATTCGTATATTAGTTTAATTGATAATTTACa                              >  2:294880/1‑141 (MQ=255)                                                                                                        tGATATTAAACCTTATATCACTGAATTACATTATGATAAGACGCCTAAAATTCTAAAAAATGATCATTTATATCACTTGTACACTGAAATAAAGGCGCCACTGAAGGATGATTCGTATATTAGTTTAATTGATAATTTACa                              >  1:101033/1‑141 (MQ=255)                                                                                                                aaCCTTATATCACTGAATTACATTATGATAAGACGCCTAAAATTCTAAAAAATGATCATTTATATCACTTGTACACTGAAATAAAGGCGCCACTGAAGGATGATTCGTATATTAGTTTAATTGATAATTTACATCCAacac                      >  2:84123/1‑141 (MQ=255)                                                                                                                  ccTTATATCACTGAATTACATTATGATAAGACGCCTAAAATTCTAAAAAATGATCATTTATATCACTTGTACACTGAAATAAAGGCGCCACTGAAGGATGATTCGTATATTAGTTTAATTGATAATTTACATCCAACACCt                    <  2:433364/141‑1 (MQ=255)                                                                                                                   cTTATATCACTGAATTACATTATGATAAGACGCCTAAAATTCTAAAAAATGATCATTTATATCACTTGTACACTGAAATAAAGGCGCCACTGAAGGATGATTCGTATATTAGTTTAATTGATAATTTACATCCAACACCTg                   <  1:84123/141‑1 (MQ=255)                                                                                                                    ttATATCACTGAATTACATTATGATAAGACGCCTAAAATTCTAAAAAATGATCATTTATATCACTTGTACACTGAAATAAAGGCGCCACTGAAGGATGATTCGTATATTAGTTTAATTGATAATTTACATCCAACACCTGc                  <  1:356367/141‑1 (MQ=255)                                                                                                                      atatCACTGAATTACATTATGATAAGACGCCTAAAATTCTAAAAAATGATCATTTATATCACTTGTACACTGAAATAAAGGCGCCACTGAAGGATGATTCGTATATTAGTTTAATTGATAATTTACATCCAACACCTGCtt                <  2:26994/141‑1 (MQ=255)                                                                                                                          cACTGAATTACATTATGATAAGACGCCTAAAATTCTAAAAAATGATCATTTATATCACTTGTACACTGAAATAAAGGCGCCACTGAAGGATGATTCGTATATTAGTTTAATTGATAATTTACATCCAACACCTGCTTTAgg            >  2:287263/1‑141 (MQ=255)                                                                                                                              gAATTACATTATGATAAGACGCCTAAAATTCTAAAAAATGATCATTTATATCACTTGTACACTGACATAAAGGCGCCACTGAAGGATGATTCGTATATTAGTTTAATTGATAATTTACATCCAACACCTGCTTTAGGTGGc        <  1:98375/141‑1 (MQ=255)                                                                                                                              gAATTACATTATGATAAGACGCCTAAAATTCTAAAAAATGATCATTTATATCACTTGTACACTGAAATAAAGTCGCCACTGAAGGATGATTCGTATATTAGTTTAATTGATAATTTACATCCAACACCTGCTTTAGGTGGc        <  1:425794/141‑1 (MQ=255)                                                                                                                              gAATTACATTATGATAAGACGCCTAAAATTCTAAAAAATGATCATTTATATCACTTGTACACTGAAATAAAGGCGCCACTGAAGGATGATTCGTATATTAGTTTAATTGATAATTTACATCCAACACCTGCTTTAGGTGGc        <  1:394337/141‑1 (MQ=255)                                                                                                                              gAATTACATTATGATAAGACGCCTAAAATTCTAAAAAATGATCATTTATATCACTTGTACACTGAAATAAAGGCGCCACTGAAGGATGATTCGTATATTAGTTTAATTGATAATTTACATCCAACACCTGCTTTAGGTGGc        <  1:277145/141‑1 (MQ=255)                                                                                                                                    cATTATGATAAGACGCCTAAAATTCTAAAAAATGATCATTTATATCACTTGTACACTGAAATAAAGGCGCCACTGAAGGATGATTCGTATATTAGTTTAATTGATAATTTACATCCAACACCTGCTTTAGGTGGCTATCCa  <  1:457440/141‑1 (MQ=255)                                                                                                                                          |                                                                                                                                       AACGTTCACAAGATGAGGACGAAGATACAAAAAATGTTGAAGCATTTTTAAAAGATAATAAAAACTTAATCGAACATCGATTTGTTGTTGACAGTATTTTACATGATATTAAACCTTATATCACTGAATTACATTATGATAAGACGCCTAAAATTCTAAAAAATGATCATTTATATCACTTGTACACTGAAATAAAGGCGCCACTGAAGGATGATTCGTATATTAGTTTAATTGATAATTTACATCCAACACCTGCTTTAGGTGGCTATCCA  >  CP000730/1047818‑1048089

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