Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	2,676,343	G→C	8.7%	Y62* (TAC→TAG)	ilvB2 ←	pyruvate oxidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	2,676,343	0	G	C	8.7%	132.3 / 6.5	47	Y62* (TAC→TAG)	ilvB2	pyruvate oxidase
Reads supporting (aligned to +/- strand):  ref base G (16/26);  new base C (2/2);  total (19/28)
Fisher's exact test for biased strand distribution p-value = 6.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.44e-01

GCTGTACTATTCGTTTGTCCAGATAATATTAATTGCGGTACATTATCCATTTTGGCATCATACATACCATTTAATAAATGAATTAAACCAGGGCCACCGATACTTAATGCCACACCGATTTTACCAGTTAATTTTGTGTAACCAGCAGCCGCTAAGCTTGCTACTTCTTCATGACGTACATGATAAAATTTAAATTGATCTCTCACTGTACGTAAACTATCGACTACTGCGTCGATTGAGTCTCCTGGAATACCATACAAGTGATCTA  >  CP000730/2676206‑2676473                                                                                                                                          |                                                                                                                                   gCTGTACTATTCGTTTGTCCAGATAATATTAATTGCGGTACATTATCCATTTTGGCATCATACATACCATTTAATAAATGAATTAAACCAGGGCCACCGATACTTAATGCCACACCGATTTTACCAGTTAATTTTGTGTaa                                                                                                                                 >  1:75396/1‑141 (MQ=255)    gTACTATTCGTTTGTCCAGATAATATTAATTGCGGTACATTATCCATTTTGGCATCATACATACCATTTAATAAATGAATTAAACCAGGGCCACCGATACTTAATGCCCCACCGATTTTACCAGTTAATTTTTTTTAAcaa                                                                                                                              >  2:8239/1‑139 (MQ=255)    gTACTATTCGTTTGTCCAGATAATATTAATTGCGGTACATTATCCATTTTGGCATCATACATACCATTTAATAAATGAATTAAACCAGGGCCACCGATACTTAATGCCACACCGATTTTACCAGTTAATTTTGTGTAACca                                                                                                                              >  1:113738/1‑141 (MQ=255)    gTACTATTCGTTTGTCCAGATAATATTAATTGCGGTACATTATCCATTTTGGCATCATACATACCATTTAATAAATGAATTAAACCAGGGCCACCGATACTTAATGCCACACCGATTTTACCAGTTAATTTTGTGTAACca                                                                                                                              >  2:302892/1‑141 (MQ=255)        tATTCGTTTGTCCAGATAATATTAATTGCGGTACATTATCCATTTTGGCATCATACATACCATTTAATAAATGAATTAAACCAGGGCCACCGATACTTAATGCCACACCGATTTTACCAGTTAATTTTGTGTAACcagcag                                                                                                                          <  1:341169/141‑1 (MQ=255)         aTTCGTTTGTCCAGATAATATTAATTGCGGTACATTATCCATTTTGGCATCATACATACCATTTAATAAATGAATTAAACCAGGGCCACCGATACTTAATGCCACACCGATTTTACCAGTTAATTTTGTGTAACCagcagc                                                                                                                         >  1:239815/1‑141 (MQ=255)           tCGTTTGTCCAGATAATATTAATTGCGGTACATTATCCATTTTGGCATCATACATACCATTTAATAAATGAATTAAACCAGGGCCACCGATACTTAATGCCACACCGATTTTACCAGTTAATTTTGTGTAACCAGCAGCCg                                                                                                                       <  2:215430/141‑1 (MQ=255)           tCGTTTGTCCAGATAATATTAATTGCGGTACATTATCCATTTTGGCATCATACATACCATTTAATAAATGAATTAAACCAGGGCCACCGATACTTAATGCCACACCGATTTTACCAGTTAATTTTGTGTAACCAGCAGCCg                                                                                                                       <  1:103084/141‑1 (MQ=255)           tCGTCTGTCCAGATAATATTAATTGCGGTACATTATCCATTTTGGCATCATACATATCATTTAATAAATGAATTCAACCTGTGCCACCGATACTTAATGCCACACCGATTTTACCAGTTAATTTTGTGTAACCAGCAGCCg                                                                                                                       <  2:59923/141‑1 (MQ=255)            cGTTTGTCCAGATAATATTAATTGCGGTACATTATCCATTTTGGCATCATACATACCATTTAATAAATGAATTAAACCAGGGCCACCGATACTTAATGCCACACCGATTTTACCAGTTAATTTTGTGTAACCAGCAGCCGc                                                                                                                      <  2:178396/141‑1 (MQ=255)             gTTTGTCCAGATAATATTAATTGCGGTACATTATCCATTTTGGCATCATACATACCATTTAATAAATGAATTAAACCAGGGCCACCGATACTTAATGCCACACCGATTTTACCAGTTAGTTTTGTCTAACCAGCAGCCGcc                                                                                                                     >  2:141799/1‑140 (MQ=255)             gTTTGTCCAGATAATATTAATTGCGGTACATTATCCATTTTGGCATCATACATACCATTTAATAAATGAATTAAACCAGGGCCACCGATACTTAATGCCACACCGATTTTACCAGTTAATTTTGTCTAACCAGCAGCCGCt                                                                                                                     >  2:141808/1‑141 (MQ=255)               ttGTCCAGATAATATTAATTGCGGTACATTATCCATTTTGGCATCATACATACCATTTAATAAATGAATTAAACCAGGGCCACCGATACTTAATGCCACACCGATTTTACCAGTTAATTTTGTGTAACCAGCAGCCGCTaa                                                                                                                   >  1:361720/1‑141 (MQ=255)                    cAGATAATATTAATTACGGTACATTGTCCATTTTGGCATCATACATACCATTTAATAAATGAATTAAACCAGGGCCACCGATACTTAATGCCACACCGATTTTACCAGTTAATTTTGTGTAACCAGCAGCCGCTAAGCTTg                                                                                                              >  1:186098/1‑141 (MQ=255)                     aGAAAATATTAATTGCCGTACATTACCCATTTTGGCATCATACATACCATTTAATAAATGAATTAAACCAGGGCCACCGATACTTAATGCCACACCGATTTTACCAGTTAATTTTGTGTAACCAGCAGCCGCTAAGCTTGc                                                                                                             <  1:20185/141‑1 (MQ=255)                      gATAATATTAATTGCGGTACATTATCCATTTTGGCATCATACATACCATTTAATAAATGAATTAAACCAGGGCCACCGATACTTAATGCCACACCGATTTTACCAGTTAATTTTGTGTAACCAGCAGCCGCTAAGCTTGCt                                                                                                            <  2:59025/141‑1 (MQ=255)                      aaaaataTTAATTGCGGTACATTATTCATTTTGGCATCATACATACCATTTAATAAATGAATTAAACCAGTGCCACCGATACTTAATGCCACACCGATTTTACCAGTTAATTTTGTGTAACCAGCAGCCGCTAAGCTTCCt                                                                                                            <  2:362456/138‑1 (MQ=255)                               taTTGCGGTACCTTATCCATTTTGGCATCATACATACCATTTAATAAATGAATTAAACCAGGGCCACCGATACTTAATGCCACCCGGGTTTTACCAGTTAATTTTGTGTAACCAGCAGCCGCTAAGCTTGCTACTTCTTCa                                                                                                   <  2:33709/140‑1 (MQ=255)                                       tACATTATCCATTTTGGCATCATACATACCATTTAATAAATGAATTAAACCAGGGCCACCGATACTTAATGCCACACCGATTTTACCAGTTAATTTTGTGTAACCAGCAGCCGCTAAGCTTGCTACTTCTTCATGACGTAc                                                                                           <  2:357243/141‑1 (MQ=255)                                       tACATTATCCATTTTGGCATCATACATACCATTTAATAAATGAATTAAACCAGGGCCACCGATACTTAATGCCACACCGATTTTACCAGTTAATTTTGTGTAACCAGCAGCCGCTAAGCTTGCTACTTCTTCATGACGTAc                                                                                           <  1:337507/141‑1 (MQ=255)                                                   tttGGCAGCATACATACCATTTAATAAATGAATTAAACCAGGGCCACCGATTCTTAATGCCACCCCGATTTTACCAGTTCATTTTGTGTAACCAGCAGCCGCTAAGCTTGCTACTTCTTCATGACCTACACGATAACAttt                                                                               <  2:165060/141‑1 (MQ=255)                                                   tggtgCATCAGACATACCATTTAATAAATGAATTAAACCAGGGCCACCGATACTTAATGCCACACCGATTTTACCAGTTAATTTTGTGTAACCAGCAGCCGCTAAGCTTGCTACTTCTTCATGACGTACATGATAAAAttt                                                                               <  2:182583/137‑1 (MQ=255)                                                                       ttAATAAATGAATTAAACCAGGGCCACCGATACTTAATGCCACACCGATTTTACCAGTTCATTTTGTGTAACCAGCAGCCGCTAAGCTTGCTACTTCTTCATGACGTACATGATAAAATTTAAATTGATCTCTCACTGTaa                                                           >  2:26344/1‑140 (MQ=255)                                                                        gactaaCTGAATTACACGAGGGCCCCCGATACTTCATGCCAAACCGATCTTACCAGTTAATTTTGTGTAACCAGCACTCGCTAAGCTTGCTACTTCTTCATGACGTACATGATAAAATTTAAATTGATCTCTCACTGTACg                                                          <  2:101781/138‑1 (MQ=255)                                                                          ataaATGAATTAAACCAGGGCCACCGATACTTAATGCCACACCGATTTTACCAGTTAATTTTGTGTAACCAGCAGCCGCTAAGCTTGCTACTTCTTCATGACGTACATGATAAAATTTAAATTGATCTCTCACTGTACGTa                                                        <  2:244331/141‑1 (MQ=255)                                                                           taaATGAATTAAACCAGGGCCACCGATACTTAATGCCACACCGATTTTACCAGTTAATTTTGTCTAACCAGCAGCCGCTAAGCTTGCTACTTCTTCATGACGTACATGATAAAATTTAAATTGATCTCTCACTGTACGTaa                                                       <  1:141808/141‑1 (MQ=255)                                                                           taaATGAATTAAACCAGGGCCACCGATACTTAATGCCACACCGATTTTACCAGGTAAGTTTGTCTAACCAGCAGCCGCTAAGCTTGCTACTTCTTCATGACGTACATGATAAAATTTAAATTGATCTCTCACTGTACGTaa                                                       <  1:141799/141‑1 (MQ=255)                                                                            aaaTGAATTAAACCAGGGCCACCGATACTTAATGCCACACCGATGTTACCAGTTAATTTTGTGTAACCATCAGCCGATACGATCGTCACGTGCCCACTCCGTGCGTCATAAACTTTAAATTGATCTCTCACCTTAAGTaaa                                                      >  2:262885/1‑141 (MQ=255)                                                                              aTGAATTAAACCAGGGCCACCGATACTTAATGCCACACCGATTTTACCAGTTAATTTTGTGTAACCAGCAGCCGCTAAGCTTGCTACTTCTTCATGACGTACATGATAAAATTTAAATTGATCTCTCACTGTACGTAAACt                                                    <  1:181424/141‑1 (MQ=255)                                                                                    tAAACCAGGGCCACCGATACTTAATGCCACACCGATTTTACCAGTTAATTTTGTGTAACCAGCAGCCGCTAAGCTTGCTACTTCTTCATGACGCACATGATAAAATTTAAATTGATCTCTCACTGTACGTATACTATCGac                                              <  2:239815/141‑1 (MQ=255)                                                                                     aaaCCAGGGCCACCGATACTTAATGCCACACCGATTTTACCAGTTAATTTTTTGTAACCAGCAGCCGCTAAACTTGCTACTTCTTCATTACGTAAATGATAAAATTTAAATTGATCTCTCACTGTACGTAAACTATCGact                                             >  2:103102/1‑141 (MQ=255)                                                                                          aGGGCCACCGATACTTAATGCCACACCGATTTTACCAGTTAATTTTGTGTAACCAGCAGCCGCTAAGCTTGCTACTTCTTCATGACGTACATGATAAAATTTAAATTGATCTCTCACTGTCCGTAAACTATCGACTACTGc                                        <  2:361720/141‑1 (MQ=255)                                                                                           gggCCACCGATACTTAATGCCACACCGATTTTACCAGTTAATTTTGTGTAACCAGCAGCCGCTAAGCTTGCTACTTCTTCATTACGTACATGATAAAATTTAAATTGATCTCTCACTGTACGTAAACTATCGACTACTGCg                                       >  2:16027/1‑141 (MQ=255)                                                                                                aCCGATACTTAATGCCACACCGATTTTACCAGTTAATTTTGTGTAACCAGCAGCCGCTAAGCTTGCTACTTCTTCATGACGTACATGATAAAAATTAAATTGATCTCTCACTGTACGTAAACTATCGACTACTGCGTCGAt                                  >  1:339344/1‑141 (MQ=255)                                                                                                   gATACTTAATGCCACACCGATTTTACCAGTTAATTTTGTGTAACCATCAGCCGCTAAGCTTGCTACTTCTTCATGACGTACATGATAAAATTTAAATTGATCTCTCACTGTACGTAAACTATCGACTACTGCGTCGATTGa                               <  2:291013/141‑1 (MQ=255)                                                                                                        ttAATGCCACACCGATTTTACCAGTTAATTTTGTGTAACCAGCAGCCGCTAAGCTTGCTACTTCTTCATGACGTACATGATAAAATTTAAATTGATCTCTCACCGAACGCAAACTATCGACCACCGCGTCGATTTAAtctc                          >  2:310052/1‑141 (MQ=255)                                                                                                          aaTGCCACACCGATTTTACCAGTTAATTTTGTGTAACCAGCAGCCGCTAAGCTTGCTACTTCTTCATGACGTACATGATAAAATTTAAATTGATCTCTCACTGTACGTAAACTATCGACTACTGCGTCGATTGAGTCTCCt                        <  1:304871/141‑1 (MQ=255)                                                                                                           attccACCCCGTTTTTACCAGTTAATTTTGTGTAACCAGCAGCCGCTAAGCTTGCTACTTCTTCATGACGTACATGATAAAATTTAAATTGATCTCTCACTGTACGTAAACTATCGACTACTGCGTCGATTGAGTCTCCTg                       <  2:161569/138‑1 (MQ=255)                                                                                                               cacaCCGATTTTACCAGTTAATTTTGTGTAACCACCAGCAGCTAATCTTGGTACGTCTTCATGACGTACATGATAAAATTTAAATTGATCTCTCACTGTACGTAACCTCTCGCCTCCTGCGTCGATTGAGTCTCCTgcact                   >  2:280565/1‑137 (MQ=255)                                                                                                                   ccGATTTTACCAGTTAATTTTGTGTAACCAGCAGCCGCTAAGCTTGCTACTTCTTCATGACGTACATGATAAAATTTAAATTGATCTCTCCCCGTACGTAAACTATCGACTACTGCGTCGCTTGAGTCTCCTGGACTAccc               >  1:209405/1‑140 (MQ=255)                                                                                                                   ccGATTTTACCAGTTAATTTTGTGTAACCAGCAGCCGCTAAGCTTGCTACTTCTTCATGACGTACATGATAAAATTTAAATTGATCTCTCACTGTACGTAAACTATCGACTACTGCGTCGATTGAGTCTTCTGGAATACCa               >  1:241193/1‑141 (MQ=255)                                                                                                                    cGATTTGACCAGTTAATTTGGTGTAACCAGCAGCCGCTAAGCTTGCTACTTCTTCATGACGTACATGATAAAATTTAAATTGATCTCTCACTGTACGTAAACTATCTACTACTGCGTCGATTGAGTCTCCTGGAATACCAt              <  2:183638/141‑1 (MQ=255)                                                                                                                      tttgtACCAGTTAATTTTGTGTAACCAGCAGCCGCTAAGCTTGCTACTTCCTCATGACGTACATGATAAAATTTAAATTGATCTCTCACTGTACGTAAACTATCGACTACTGCGTCGATTGAGTCTCCTGGAATACCATAc            <  1:107403/137‑1 (MQ=255)                                                                                                                       ttttACCAGTTAATTTTGTGTAACCAGCAGCCGCTAAGCTTGCTACTTCTTCATGACGTACATGATAAAATTTAAATTGATCTCTCACTGTACGTAAACTATCGACTACTGTGGGGTTTGAGTCTCCTGGAATACCATACa           <  1:180177/141‑1 (MQ=255)                                                                                                                         ttACCAGTTAATTTTGTGTAACCAGCAGCCGCTAAGCTTGCTACTTCTTCATGACGTACATGATAAAATTTAAATTGATCTCTCACTGTACGTAAACTATCGACTACTGCGTCGATTGAGTCTCCTGGAATACCATACAAg         >  1:41265/1‑141 (MQ=255)                                                                                                                           aCCAGTTAATTTTGTGTAACCAGCAGCCGCTAAGCTTGCTACTTCTTCATGACGTACATGATAAAATTTAAATTGATCTCTCACTGTACGTAAACTATCGACTACTGCGTCGATTGAGTCTCCTGGAATACCATACAAGTg       <  1:59346/141‑1 (MQ=255)                                                                                                                                ttAATTTTGTGTAACCAGCAGCCGCTAAGCTTGCTACTTCTTCATGACGTACATGATAAAATTTAAATTGATCTCTCACTGTACGTAAACTATCGACTACTGCGTCGATTGAGTCTCCTGGAATACCATACAAGTGATCta  <  1:268591/141‑1 (MQ=255)                                                                                                                                          |                                                                                                                                   GCTGTACTATTCGTTTGTCCAGATAATATTAATTGCGGTACATTATCCATTTTGGCATCATACATACCATTTAATAAATGAATTAAACCAGGGCCACCGATACTTAATGCCACACCGATTTTACCAGTTAATTTTGTGTAACCAGCAGCCGCTAAGCTTGCTACTTCTTCATGACGTACATGATAAAATTTAAATTGATCTCTCACTGTACGTAAACTATCGACTACTGCGTCGATTGAGTCTCCTGGAATACCATACAAGTGATCTA  >  CP000730/2676206‑2676473

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