Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	336,153	A→G	9.1%	G97G (GGA→GGG)	essC →	virulence protein EssC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	336,153	0	A	G	9.1%	110.6 / 4.0	44	G97G (GGA→GGG)	essC	virulence protein EssC
Reads supporting (aligned to +/- strand):  ref base A (21/19);  new base G (2/2);  total (23/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.20e-01

TTGGCAAGCGAATCATACTTCTATTAATAAGGTGCTTGTTAGAAAAGGTGACCTTGATGACATTACATTACAGCTTTATACAGAAGCTGATTATGCATCATTTGCGTATCCTTCAATTCAAGATACGATGACAATTGGACCAAATGCGTATGATGATATGGTTATTCAAAGCTTGATGAATGCCATCATTATTAAAGATTTTCAATCAATACAAGAATCACAATACGTACGCATTGTGCACGATAAAAATACAGATGTGTATATTAACTATGAA  >  CP000730/336015‑336288                                                                                                                                           |                                                                                                                                        ttGGCAAGCGAATCATACTTCTATTAATAAGGTGCTTGTTAGAAAAGGTGACCTTGATGACATTACATTACAGCTTTATACAGAAGCTGATTATGCATCATTTGCGTATCCTTCAATTCAAGATACGATGACAATTGGAcc                                                                                                                                       >  1:114495/1‑141 (MQ=255)    gCAAGCGAATCATACTTCTATTAATAAGGTGCTTGTTAGAAAAGGTGACCTTGATGACATTACATTACAGCTTTATACAGAAGCTGATTATGCATCATTTGCGTATCCTTCAATTCAAGATACGATGACAATTGAACCaaa                                                                                                                                    >  2:191898/1‑141 (MQ=255)     cAAGCGAATCATACTTCTATTAATAAGGTGCTTGTTAGAAAAGGTGACCTTGATGACATTACATTACAGCTTTATACAGAAGCTGATTATGCATCATTTGCGTATCCTTCAATTCAAGATACGATGACAATTGGGCCAAAt                                                                                                                                   <  1:345879/141‑1 (MQ=255)     cAAGCGAATCATACTTCTATTAATAAGGTGCCGGTTAGAAAAGGTGACCTTGATGACATTACATTACAGCTTTATACAGAAGCGGATTATGCATCATTTGCGTATCCTTCAATTCAAGATACGATGACAATTGGACCAAAt                                                                                                                                   >  2:379569/1‑141 (MQ=255)      aaGCGAATCATACTTCTATTAATAAGGTGCTTGTTAGAAAAGGTGACCTTGATGACATTACATTACAGCTTTATACAGAAGCTGATTATGCATCATTTGCGTATCCTTCACTTCAAGATACGATGACAATTGGACCAAATg                                                                                                                                  <  1:409923/141‑1 (MQ=255)          gAATCATACTTCTATTAATAAGGTGCTTGTTAGAAAAGGTGACCTTGATGACATTACATTACAGCTTTATACAGAAGCTGATTATGCATCATTTGCGTATCCTTCAATTCAAGATACGATGACAATTGGACCAAATGCGTa                                                                                                                              >  1:59203/1‑141 (MQ=255)           aaTCATACTTCTATTAATAAGGTGCTTGTTAGAAAAGGTGACCTTGATGACATTACATTACAGCTTTATACAGAAGCTGATTATGCATCATTTGCGTATCCTTCAATTCAAGATACGATGACAATTGGACCAAATGCGTat                                                                                                                             <  1:167161/141‑1 (MQ=255)            aTCATACTTCTATTAATAAGGTGCTTGTTAGAAAAGGTGACCTTGATGACATTACATTACAGCTTTATACAGAAGCTGATTATGCATCATTTGCGTATCCTTCAATTCAAGATACGATGACAATTGGACCAAATGCGTatg                                                                                                                            <  2:182070/141‑1 (MQ=255)            aTCACACTTCTATTAATAAGGTGCTTGTTAGAAAAGGTGACCTTGATGACATTACATTACAGCTTTATACAGAAGCTGATTATGCATCATTTGCGTATCCTTCAATTCAAGATACGATGACAATTGGACCAAATGCGTatg                                                                                                                            <  1:411810/141‑1 (MQ=255)                  cTTCTATTAATAAGGTGCTTGTTAGAAAAGGTGACCTTGATGACATTACATTACAGCTTTATACAGAAGCTGATTATGCATCATTTGCGTATCCTTCAATTCAAGATACGATGACAATTGGACCAAATGCGTATGATGata                                                                                                                      <  1:216783/141‑1 (MQ=255)                     cTATTAATAAGGTGCTTGTTAGAAAAGGTGACCTTGATGACATTACATTACAGCTTTATACAGAAGCTGATTATGCATCATTTGCGTATCCTTCAATTCAAGATACGATGACAATTGGACCAAATGCGTATGATGATATgg                                                                                                                   <  2:320371/141‑1 (MQ=255)                      tATTAATAAGGTGCTTGTTAGAAAAGGTGACCTTGATGACATTACATTACAGCTTTATACAGAAGCTGATTATGCATCATTTGCGTATCCTTCAATTCAAGATACGATGACAATTGGACCAAATGCGTATGATGATATGGt                                                                                                                  >  1:50914/1‑141 (MQ=255)                         taataaGGTGCTTGTTAGAAAAGGTGACCTTGATGACATTACATTACAGCTTTATACAGAAGCTGATTATGCATCATTTGCGTATCCTTCAATTCAAGATACGATGACAATTGGACCAAATGCGTATGATGATATGGTTAt                                                                                                               >  2:259176/1‑141 (MQ=255)                         taataaGGTGCTTGTTAGAAAAGGTGACCTTGATGACATTACATTACAGCTTTATACAGAAGCTGATTATGCATCATTTGCGTATCCTTCAATTCAAGATACGATGACAATTGGACCAAATGCGTATGATGATATGGTTAt                                                                                                               <  2:408881/141‑1 (MQ=255)                         taataaGGTGCTTGTTAGAAAAGGTGACCTTGATGACATTACATTACAGCTTTATACAGAAGCTGATTATGCATCATTTGCGTATCCTTCAATTCAAGATACGATGACAATTGGACCAAATGCGTATGATGATATGGTTAt                                                                                                               <  2:60904/141‑1 (MQ=255)                           ataaGGTGCTTGTTAGAAAAGGTGACCTTGATGACATTACATTACAGCTTTATACCGAATCTGATTATGCATCATTTGCGTATCCTTCAATTCAAGATACGATGACAATTGGACCAAATGCGTATGATGATATGGTTATTc                                                                                                             <  1:96976/141‑1 (MQ=255)                             aaGGTGCTTGTTAGAAAAGGTGACCTTGATGACATTACATTACAGCTTTATACAGAAGCTGATTATGCATCATTTGCGTATCCTTCAATTCAAGATACGATGACAATTGGACCAAATGCGTATGAGGATATGGTTATTCaa                                                                                                           >  2:86807/1‑141 (MQ=255)                                gTGCTTGTTAGAAAAGGTGACCTTGATGACATTACATTACAGCTTTATACAGAAGCTGATTATGCATCATTTGCGTATCCTTCAATTCAAGATACGATGACAATTGGACCAAATGCGTATGATGATATGGTTATTCAAAGc                                                                                                        >  1:328223/1‑141 (MQ=255)                                gTGCTTGTTAGAAAAGGTGACCTTGATGACATTACATTACAGCTTTATACAGAAGCTGATTATGCATCATTTGCGTATCCTTCAATTCAAGATACGATGACAATTGGACCAAATGCGTATGATGATATGGTTATTCAAAGc                                                                                                        <  1:257996/141‑1 (MQ=255)                                 tGCTTGTTAGAAAAGGTGACCTTGATGACATTACATTACAGCTTTATACAGAAGCTGATTATGCATCATTTGCGTATCCTTCAATTCAAGATACGATGACAATTGGACCAAATGCGTATGATGATATGGTTATTCAAAGCt                                                                                                       <  1:273317/141‑1 (MQ=255)                                             aaGGTGACCTTGATGACATTACATTACAGCTTTATACAGAAGCTGATTATGCATCATTTGCGTATCCTTCAATTCAAGATACGATGACAATTGGACCAAATGCGTATGATGATATGGTTATTCAAAGCTTGATGAATGCca                                                                                           >  1:357402/1‑141 (MQ=255)                                              aGGTGACCTTGATGACATTACATTACAGCTTTATACAGAAGCTGATTATGCATCATTTGCGTATCCTTCAATTCAAGATACGATGACAATTGGACCAAATGCGTATGATGATATGGTTATTCAAAGCTTGATGAATGCcat                                                                                          >  1:44751/1‑141 (MQ=255)                                               ggagACCTTGCTGTCGTTACATGACAGCATGATACAGAAGCTGGTTAAGCATCTTTTGCGGGTCCGTCAATTCAAGATACGATGACAATTGGACCAAATGCGTATGATTATTTGTGTATTCAAAGCTTTTTAAATGCcatc                                                                                         >  2:78093/4‑141 (MQ=255)                                               ggTGACCTTGATGACATTACATTACAGCTTTATACAGAAGCTGATTATGCATCATTTGCGTATCCTTCAATTCAAGATACGATGACAATTGGACCAAATGCGTATGATGATATGGTTATTCAAAGCTTGATGAATGCcatc                                                                                         >  2:397809/1‑141 (MQ=255)                                                     cTTGATGACATTACATTACAGCTTTATACAGAAGCTGATTATGCATCATTTGCGTATCCTTCAATTCAAGATACGATGACAATTGGACCAAATGCGTATGATGATATGGTTATTCAAAGCTTGATGAATGCCATCattatt                                                                                   <  2:445876/141‑1 (MQ=255)                                                        gatgaCATTACATTACTGCTTTATCCAGAGGCTGATTATGCATCATTTGCGTATCCTTCAATTCAAGATACGTTGACAATTGGATCAAATGCGTATGATGATATGGTTATTCAAAGCTTGATGAATGCCATAATTATTaaa                                                                                <  2:29912/141‑1 (MQ=255)                                                                        aGCTTTATACAGAAGCTGATTATGCATCATTTGCGTATCCTTCAATTCAAGATACGATGACAATTGGACCAAATGCGTATGATGATATGGTTATTCAAAGCTTGATGAATGCCATCATTATTAAAGATTTTCAATCAATAc                                                                >  1:40230/1‑141 (MQ=255)                                                                        aGCTTTATACAGAAGCTGATTATGCATCATTTGCGTATCCTTCAATTCAAGATACGATGACAATTGGACCAAATGCGTATGATGATATGGTTATTCAAAGCTTGATGAATGCCATCATTATTAAAGATTTTCAATCAATAc                                                                >  1:186861/1‑141 (MQ=255)                                                                               gtataaaGCTGATTATGCATCATTTGCGTATCCTTCAATTCAAGATACGATGACAATTGGGCCAAATGCGTATGATGTTATGGTTATTCAAAGCTTGATGAATGCCATCATTATTAAAGATTTTCAATCAATACAAGAATc                                                         >  2:101851/6‑141 (MQ=255)                                                                                 cAGAAGCTGATTATGCATCATTTGCGTATCCTTCAATTCAAGATACGATGACAATTGGACCAACTGCGTATGATGATATGGCTATTCAAAGCTTGATCCATGCCATCATTATTAAAGATTTTCAATCAATACAAGAATcac                                                       >  2:207021/1‑141 (MQ=255)                                                                                         gATTATGCATCATTTGCGTATCCTTCAATTCAAGATACGATGACAATTGGACCAAATGCGTATGATGATATGGTTATTCAAAGCTTGATGAATGCCATCATTATTAAAGATTTTCAATCAATACAAGAATCACAAtacgta                                               >  2:322001/1‑141 (MQ=255)                                                                                           ttATGCATCATTTGCGTATCCTTCAATTCAAGATACGATGACAATTGGACCAAATGCGTATGATGATATGGTTATTCAAAGCTTGATGATTGCCATCATTATTAAAGATTTTCAATCAATACAAGAATCACAAtacgtacg                                             <  1:349268/141‑1 (MQ=255)                                                                                           ttATGCATCATTTGCGTATCCTTCAATTCAAGATACGATGACAATTGGACCAAATGCGTATGATGATATGGTTATTCAAAGCTTGATGAATGCCATCATTATTAAAGATTTTCAATCAATACAAGAATCACAAtacgtacg                                             >  1:243252/1‑141 (MQ=255)                                                                                           ttATGCATCATTTGCGTATCCTTCAATTCAAGATACGATGACAATTGGACCAAATGCGTATGATGATATGGTTATTCAAAGCTTGATGAATGCCATCATTATTAAAGATTTTCAATCAATACAAGAATCACAAtacgtacg                                             <  2:113339/141‑1 (MQ=255)                                                                                                        gCGTATCCTTCAATTCAAGATACGATGACAATTGGACCAAATGCGTATGATGATATGGTTATTCAAAGCTTGATGAATGCCATCATTATTAAAGATTTTCAATCAATACAAGAATCACAATACGTACGCATTGTGCACGAt                                <  1:38255/141‑1 (MQ=255)                                                                                                           tATCCTTCAATTCAAGATACGATGACAATTGGGCCAAATGCGTATGATGATATGGTTATTCAAAGCTTGATGAATGCCATCATTATTAAAGATTTTCAATCAATACAAGAATCACAATACGTACGCATTGTGCACGATaaa                             >  1:203542/1‑141 (MQ=255)                                                                                                           tATCCTTCAATTCAAGATACGATGACAATTGGACCAAATGCGTATGATGATATGGTTATTCAAAGCTTGATGAATGCCATCATTATTAAAGATTTTCAATCAATACAAGAATCACAATACGTACGCATTGTGCACGATaaa                             <  1:326808/141‑1 (MQ=255)                                                                                                                         aGATACGATGACAATTGGACCAAATGCGTATGATGATATGGTTATTCAAAGCTTGATGAATGCCATCATTATTAAAGATTTTCAATCAATACAAGAATCACAATACGTACGCATTGTGCACGATAAAAATACAGATGTGta               >  2:259758/1‑141 (MQ=255)                                                                                                                         aGATACGATGACAATTGGACCAAATGCGTATGATGATATGGTTATTCAAAGCTTGATGAATGCCATCATTATTAAAGATTTTCAATCAATACAAGAATCACAATACGTACGCATTGTGCACGATAAAAATACAGATGTGta               >  2:220884/1‑141 (MQ=255)                                                                                                                          gATACGATGACAATTGGACCAAATGCGTATGATGATATGGTTATTCAAAGCTTGATGAATGCCATCATTATTAAAGATTTTCAATCAATACAAGAATCACAATACGTACGCATTGTGCACGATAAAAATACAGATGTGtat              >  1:209016/1‑141 (MQ=255)                                                                                                                              cGATGACAATTGGACCAAATGCGTATGATGATATGGTTATTCAAAGCTTGATGAATGCCATCATTATTAAAGATTTTCAATCAATACAAGAATCACAATACGTACGCATTGTGCACGATAAAAATACAGATGTGTATATTa          >  1:205116/1‑141 (MQ=255)                                                                                                                                    cAATTGGGCCAAATGCGTATGATGATATGGTTATTCAAAGCTTGATGAATGCCATCATTATTAAAGATTTTCAATCAATACAAGAATCACAATACGTACGCATTGTGCACGATAAAAATACAGATGTGTATATTAACTATg    <  1:431061/141‑1 (MQ=255)                                                                                                                                    cAATTGGACCAAATGCGTATGATGATATGGTTATTCAAAGCTTGATGAATGCCATCATTATTAAAGATTTTCAATCAATACAAGAATCACAATACGTACGCATTGTGCACGATAAAAATACAGATGTGTATATTAACTATg    <  1:383605/141‑1 (MQ=255)                                                                                                                                      aTTGGACCAAATGCGTATGATGATATGGTTATTCAAAGCTTGATGAATGCCATCATTATTAAAGATTTTCAATAAATACAAGAATCACAATACGTACGCATTGTGCACGATAAAAATACAGATGTGTATATTAACTATGaa  <  1:178141/141‑1 (MQ=255)                                                                                                                                           |                                                                                                                                        TTGGCAAGCGAATCATACTTCTATTAATAAGGTGCTTGTTAGAAAAGGTGACCTTGATGACATTACATTACAGCTTTATACAGAAGCTGATTATGCATCATTTGCGTATCCTTCAATTCAAGATACGATGACAATTGGACCAAATGCGTATGATGATATGGTTATTCAAAGCTTGATGAATGCCATCATTATTAAAGATTTTCAATCAATACAAGAATCACAATACGTACGCATTGTGCACGATAAAAATACAGATGTGTATATTAACTATGAA  >  CP000730/336015‑336288

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