Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	2,123,156	T→C	16.3%	N402D (AAT→GAT)	wcaD ←	putative colanic acid polymerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	2,123,156	0	T	C	16.3%	55.4 / 3.0	26	N402D (AAT→GAT)	wcaD	putative colanic acid polymerase
Reads supporting (aligned to +/- strand):  ref base T (12/9);  new base C (2/2);  total (14/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACAACAATCCATTCGAAGCTGATATCTTCCACCTGCGCCAGATGCGCCAGCGAGGCATGTGTTTTGACTATCCCTTCGAGGTTACGAAACGCGACAGTGATTATGCTAAGCAACATGTTCTTATTCTTACCTCGTAATATTTAACGCTTTTCGCAAAATAAACGGACAGACGATTAAAAATGCATATTCCGGGCTAAATATCGAACCGGTAAAAAACAGCGATACCGGAGTAAAAAGATATAATTGCACGCGAAAATTGCGGTTATCACCAAAAGCGTTGA  >  CP009273/2123016‑2123296                                                                                                                                             |                                                                                                                                             acaacaATCCATTCGAAGCTGATATCTTCCACCTGCGCCAGATGCGCCAGCGCGGCATGTGTTTTGACTATCCCTTCGAGGTTACGAAACGCCACAGTGATTATGCTAAGCAACATGTTCTTCTTCTTCCCTCGCAATATTTAACGCtt                                                                                                                                      >  1:313377/1‑149 (MQ=255)         ccATTCGAAGCTGATATCTTCCACCTGCGCCAGATGCGCCAGCGAGGCATGTGTTTTGACTATCCCTTCGAGGTTACGAAACGCGACAGTGATTATGCTAAGCAACATGTTCTTATTCTTACCTCGTAATATTTAACGCTTTTCGCaaa                                                                                                                              <  1:152585/149‑1 (MQ=255)           aTTCGAAGCTGATATCTTCCACCTGCGCCAGATGCGCCAGCGAGGCATGTGTTTTGACTATCCCTTCGAGGTTACGAAACGCGACAGTGATTATGCTAAGCAACATGTTCTTATTCTTACCTCGTAATATTTAACGCTTTTCg                                                                                                                                  <  2:2081/143‑1 (MQ=255)                    tGATATCTTCCACCTGCGCCAGATGCGCCAGCGAGGCATGTGTTTTGACTATCCCTTCGAGGTTACGAAACGCGACAGTGATTATGCTAAGCAACATGTTCTTATTCTTACCTCGTACTATCTAACGCTTTTCGCAAAATAAACGGACa                                                                                                                   >  1:247430/1‑149 (MQ=255)                    tGATATCTTCCACCTGCGCCAGATGCGCCAGCGAGGCATGTGTTTTGACTATCCCTTCGAGGTTACGAAACGCGACAGTGATTATGCTAAGCAACATGTTCTTATTCTTACCTCGTACTATCTAACGCTTTTCGCAAAATAAACGGACa                                                                                                                   >  1:247439/1‑149 (MQ=255)                                    cgcCAGATGCGCCAGCGCGGCATGTGTTTTGCCTATCCCTTCCAGGTTCCGAAACGCCACAGTGATTATGCTACGCAACATGTTCTTATTCTTACCCCGTAATATTTACCGCTTCTCGCAAAATAAACGGACAGACCATTaaa                                                                                                         <  2:168962/143‑1 (MQ=255)                                        aGATGCGCCAGCGAGGCATGTGTTTTGACTATCCCTTCGAGGTTACGAAACGCGACAGTGATTATGCTAAGCAATGTATTGCAGTACTTGCCTCGTAATAGATAACGCTTTTCGCAAAATAAACGGACAGACGATTAAAAATg                                                                                                     <  2:318533/143‑1 (MQ=255)                                        aGATGCGCCAGCGAGGCATGTGTTTTGACTATCCCTTCGAGGTTACGAAACGCGACAGTGATTATGCTAAGCAACATGTTCTTATTCTTACCTCGTAATATTTAACGCTTTTCGCAAAATAAACGGACAGACGATTAAAAATg                                                                                                     <  2:199066/143‑1 (MQ=255)                                         gATGGGCCAGCGAGGCATGTGTTTTGACTATCCCTTCGAGGTTACGAAACGCGACAGTGATTATGCTAAGCAACATGTTCTTATTCTTACCTCGTAATATTTAACGCTTTTCGCAAAATAAACGGACAGACGATTAAAAATGc                                                                                                    <  2:282734/143‑1 (MQ=255)                                                cAGCGAGGCATGTGTTTTGACTATCCCTTCGAGGTTACGAAACGCGACAGTGATTATGCTAAGCAACATGTTCTTATTCTTACCTCGTAATATTTAACGCTTTTCGCAAAATAAACGGACAGACGATTAAAAATGCATATTcc                                                                                             >  2:232078/1‑143 (MQ=255)                                                            tgtTTTGACTATCCCTTCGAGGTTACGAAACGCGACAGTGATTATGCTAAGCAACATGTTCTTATTCTTACCTCGTAATATTTAACGCTTTTCGCAAAATAAACGGACAGACGATTAAAAATGCATATTCCGGGCTAAATATCGAACCg                                                                           >  1:365809/1‑149 (MQ=255)                                                            tgtTTTGACTATCCCTTCGAGGGTACGAAACGCGACAGTGATTATGCTAAGCACCATGTTCTTACGCTTACCTCGTAATATTTAACGCTTTTCGCAAAATAAACGGACAGACGATTAAAAATGCATATTCCGGGCTAAATATc                                                                                 >  2:149622/1‑143 (MQ=255)                                                              ttttGACTATCCCTTCGAGGTTACGAAACGCGACAGTGATTATGCTAAGCAACATGTTCTTATTCTTACCTCGTAATATTTAACGCTTTTCGCAAAATAAACGGACAGACGATTAAAAATGCCTATTCCGGGCTAAATAACGa                                                                               >  2:11087/1‑143 (MQ=255)                                                               tttGACTATCCCTTCGAGGTTACGAAACGCGACAGTGATTATGCTAAGCAACATGTTCTTATTCTTACCTCGTAATATTTAACGCTTTTCGCAAAATAAACGGACAGACGATTAAAAATAAATATTCAGGGTTAAATATCGaaggcgta                                                                        >  1:58048/1‑143 (MQ=255)                                                                    cTATCCCTTCGAGGTTACGAAACGCGACAGTGATTATGCTAAGCAACATGTTCTTATTCTTACCTCGTAATATTTAACGCTTTTCGCAAAATAAACGGACAGACGATTAAAAATGCATATTCCGGGCTAAATATCGAACCGGTaaaaaa                                                                   <  1:232292/149‑1 (MQ=255)                                                                          cgtCGAGGTTACGAAACGCGACAGTGATTATGCTAAGCAACATGTTCTTATTCTTACCTCGTAATATTTAACGCTTTTCGCAAAATAAACGGACAGACGATTAAAAATGCATATTCCGGGCTAAATATCGAACCGGTaaaaaa                                                                   <  2:194453/141‑1 (MQ=255)                                                                                                aGTGATTATGCTAAGCAACATGTTCTTATTCTTACCTCGTAATATTTAACGCTTTTCGCAAAATAAACGGACAGACGATTAAAAATGCATATTCCGGGCTAAATATCGAACCGGTAAAAAACAGCGCAACCGGGGAAAAAAGa                                             >  2:181713/1‑143 (MQ=255)                                                                                                 gTGATTATGCTAAGCAACATGTTCTTATTCTTACCTCGTAATATTTAACGCTTTTCGCAAAATAAACGGACAGACGATTAAAAATGCATATTCCGGGCTAAATATCGAACCGGTAAAAAACAGCGATACCGGAGTAAAAAGATATAAtt                                      <  1:343452/149‑1 (MQ=255)                                                                                                  tGATTATGCTAAGCAACATGTTCTTATTCTTACCTCGTAATATTTAACGCTTTTCGCAAAATAAACGGACAGACGATTAAAAATGCATATTCCGGGCTAAATATCGAACCGGTAAAAAACAGCGATACCGGAGTAAAAAGata                                           >  2:105307/1‑143 (MQ=255)                                                                                                        tGCTAAGCAACATGTTCTTATTCTTACCTCGTAATATTTAACGCTTTTCGCAAAATAAACGGACAGACGATTAAAAATGCATATTCCGGGCTAAATATCGAACCGGTAAAAAACAGCGATACCGGAGTAAAAAGATATAATTGCAcgcg                               >  1:354970/1‑149 (MQ=255)                                                                                                           tAAGCAACATGTTCTTATTCTTACCTCGTAATATTTAACGCTTTTCGCAAAATAAACGGACAGACGATTAAAAATGCATATTCCGGGCTAAATATCGAACCGGTAAAAAACAGCGATACCGGAGTAAAAAGATATAATTGCaa                                  >  2:158335/1‑142 (MQ=255)                                                                                                                 aCATGTTCTTATTCTTACCTCGTAATATTTAACGCTTTTCGCAAAATAAACGGACAGACGATTAAAAATGCATATTCCGGGCTAAATATCGAACCGGTAAAAAACAGCGATACCGGAGTAAAAAGATATAATTGCACGCGAAAATTGCg                      <  1:20667/149‑1 (MQ=255)                                                                                                                       tcttattcttaCCTCGTAATATCTAACGCTTTTCGCAAAATAAACGGACAGACGATTAAAAATGCATATTCCGGGCTAAATATCGAACCGGTAAAAAACAGCGATACCGGAGTAAAAAGATATAATTGCACGCGAAAATTGCg                      <  2:220400/143‑1 (MQ=255)                                                                                                                       tcttattcttaCCTCGTAATATCTAACGCTTTTCGCAAAATAAACGGACAGACGATTAAAAATGCATATTCCGGGCTAAATATCGAACCGGTAAAAAACAGCGATACCGGAGTAAAAAGATATAATTGCACGCGAAAATTGCg                      <  2:220409/143‑1 (MQ=255)                                                                                                                                  ccTCGTAATATTTAACGCTTTTCGCAAAATAAACGGACAGACGATTAAAAATGCATATTCCGGGCTAAATATCGAACCGGTAAAAAACAGCGATACCGGAGTAAAAAGATATAATTGCACGCGCAAATTGCGGTTATCACCAAAAGCGt     >  1:130914/1‑149 (MQ=255)                                                                                                                                           atTTAACGCTTTTCGCAAAATAAACGGACAGACGATTAAAAATGCATATTCCGGGCTAAATATCGAACCGGTAAAAAACAGCGATACCGGAGTAAAAAGATATAATTGCACGCGAAAAATGCGGTTATCACCAAAAGCGTAGa  >  2:263164/1‑143 (MQ=255)                                                                                                                                             |                                                                                                                                             ACAACAATCCATTCGAAGCTGATATCTTCCACCTGCGCCAGATGCGCCAGCGAGGCATGTGTTTTGACTATCCCTTCGAGGTTACGAAACGCGACAGTGATTATGCTAAGCAACATGTTCTTATTCTTACCTCGTAATATTTAACGCTTTTCGCAAAATAAACGGACAGACGATTAAAAATGCATATTCCGGGCTAAATATCGAACCGGTAAAAAACAGCGATACCGGAGTAAAAAGATATAATTGCACGCGAAAATTGCGGTTATCACCAAAAGCGTTGA  >  CP009273/2123016‑2123296

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