Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 820,942 G→A 10.5% intergenic (+19/+170) glxK → / ← pepT glycerate kinase/M20 family peptidase T

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730820,9420GA10.5% 97.9 / 2.5 38intergenic (+19/+170)glxK/pepTglycerate kinase/M20 family peptidase T
Reads supporting (aligned to +/- strand):  ref base G (18/16);  new base A (2/2);  total (20/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.38e-01

TAATACATTTATCGATATGCCAGAAACTTATACTGACTTTAAAATGGGGTTACAAATTAGGCATTATACGGTTCAGTCTTTAAAACTGTTGAAAACACATTTTAATGTTGAGGTTTAGTAAAGAAGGACTAAATTGGTGATGCTGTCATGATGGTTAATAACATTTATGATGGTTAGCAAAACGAATTAGAAGATCGAAAGTATACGTAAAAAATATGAAAAATCACGCTATCATTGCACTGAATGTTAGCGTGATTTTTATATATTAATTAAG  >  CP000730/820806‑821079
                                                                                                                                        |                                                                                                                                         
ttaaaCATTTATCGATATGCCAGAAACTTATACTGACTTTAAAATGGGGTTACAAATTAGGCATTATACGGTTCAGTCTTTAAAACTGTTGAAAACACATTTTAATGTTGAGGTTTAGTAAAGAAGGACTAAATTGGTGAt                                                                                                                                       >  1:66377/5‑141 (MQ=255)
   tACATTTATCGATATGCCAGAAACTTATACTGACTTTAAAATGGGGTTACAAATTAGGCATTATACGGTTCAGTCTTTAAAACTGTTGAAAACACATTTTAATGTTGAGGTTTAGTAAAGAAGGACTAAATTGGTGATGCt                                                                                                                                    >  2:269109/1‑141 (MQ=255)
         tATCGATATGCCAGAAACTTATACTGACTTTAAAATGGGGTTACAAATTAGGCATTATACGGTTCAGTCTTTAAAACTGTTGAAAACACATTTTAATGTTGAGGTTTAGTAAAGAAGGACTAAATTGGTGATGCTGTCatg                                                                                                                              >  1:145777/1‑141 (MQ=255)
         tATCGATATGCCAGAAACTTATACTGACTTTAAAATGGGATTACAAATTAGGCATTATACGGTTCAGTCTTTAAAACTGTTGAAAACACATTTTAATGTTGAGGTTTAGTAAAGAAGGACTAAATTGATGATGCTGTCatg                                                                                                                              >  1:245177/1‑141 (MQ=255)
          aTCGATATGCCAGAAACTTATACTGACTTTAAAATGGGGTTACAAATTAGGCATTATACGGTTCAGTCTTTAAAACTGTTGAAAACACATTTTAATGTTGAGGTTTAGTAAAGAAGGACTAAATTGGTGATGCTGTCatga                                                                                                                             <  2:260271/141‑1 (MQ=255)
             gATATGCCAGAAACTTATACTGACTTTAAAATGGGGTTACAAATTAGGCATTATACGGTTCAGTCTTTAAAACTGTTGAAAACACATTTTAATGTTGAGGTTTAGTAAAGAAGGACTAAATTGGTGATGCTGTCATGATgg                                                                                                                          >  1:13762/1‑141 (MQ=255)
               tatGCCAGAAACTTATACTGACTTTAAAATGGGGTTACAAATTAGGCATTATACGGTTCAGTCTTTAAAACTGTTGAAAACACATTTTAATGTTGAGGTTTAGTAAAGAAGGACTAAATTGGTGATGCTGTCATGATGGtt                                                                                                                        >  2:215555/1‑141 (MQ=255)
                  gCCAGAAACTTATACTGACTTTAAAATGGGATTACAAATTAGGCATTATACGGTTCAGTCTTTAAAACTGTTGAAAACACATTATAATGTTGAGGTTTAGTAAAGAAGGACTAAATTGATGATGCTGTCATGATGGTtaat                                                                                                                     <  2:318251/141‑1 (MQ=255)
                  gCCAGAAACTTATACTGACTTTAAAATGGGATTACAAATTAGGCATTATACGGTTCAGTCTTTAAAACTGTTGAAAACACATTATAATGTTGAGGTTTAGTAAAGAAGGACTAAATTGATGATGCTGTCATGATGGTtaat                                                                                                                     <  2:318300/141‑1 (MQ=255)
                      gAAACTTATACTGACTTTAAAATGGGGTTACAAATTAGGCATTATACGGTTCAGTCTTTAAAACTGTTGAAAACACATTTTAATGTTGAGGTTTAGTAAAGAAGGACTAAATTGGTGATGCTGTCATGATGGTTAATAACa                                                                                                                 >  1:381648/1‑141 (MQ=255)
                      gAAACTTATACTGACGTTAAAATGGGGTTACAAATTAGGCATTATACGGTTCAGTCTTTGAAACTGTTGAAAACACATTTTAATGTTGAGGTTTAGTAAAGAAGGACTAAATTGGTGATGCTGTCATGATGGTTAATAACa                                                                                                                 >  2:385862/1‑141 (MQ=255)
                                cTGACTTTAAAATGGGGTTACAAATTAGGCATTATACGGTTCAGTCTTTAAAACTGTTGAAAACACATTTTAATGTTGAGGTTTAGTAAAGAAGGACTAAATTGGTGATGCTGTCATGATGGTTAATAACATTTATGATgg                                                                                                       >  2:127453/1‑141 (MQ=255)
                                 tGACTTTAAAATGGGGTTACAAATTAGGCATTATACGGTTCAGTCTTTAAAACTGTTGAAAACACATTTTAATGTTGAGGTTTAGTAAAGAAGGACTAAATTGGTGATGCTGTCATGATGGTTAATAACATTTATGATGGt                                                                                                      <  1:233632/141‑1 (MQ=255)
                                    cTTTAAAATGGGGTTACAAATTAGGCATTATACGGTTCAGTCTTTAAAACTGTTGAAAACACATTTTAATGTTGAGGTTTAGTAAAGAAGGACTAAATTGGTGATGCTGTCATGATGGTTAATAACATTTATGATGGTTAg                                                                                                   <  1:38983/141‑1 (MQ=255)
                                       tAAAAGGGCATTACAAATTAGGCATTATACGGTTCAGTCTTTAAAACTGTTGAAAACACATTTTAATGTTGAGGTTTAGTAAAGAAGGACTAAATTGGTGATGCTGTCATGATGGTTAATAACATTTATGATGGTTAGCaa                                                                                                <  2:364982/141‑1 (MQ=255)
                                        aaaaTGGGGTTACAAATTAGGCATTATACGGTTCAGTCTTTAAAACTGTTGAAAACACATTTTGGTGTTGAGGTTTAGTAAAGAAGGACTAAATTGGTGATGCTGTCATGATGGTTAATAACATTTATGATGGTTAGCaaa                                                                                               >  2:263688/1‑141 (MQ=255)
                                              gggTTACAAATTAGGCATTATACGGTTCAGTCTTTAAAACTGTTGAAAACACATTTTCTTGTTGAGGTTTAGTACGGAAGGAAAAAGTATGTGATGTTGTCATGATGGTTAATAACATTTATGATGGTTAGCAAAACGAAt                                                                                         <  1:433446/141‑1 (MQ=255)
                                              gggTTACAAATTAGGCATTATACGGTTCAGTCTTTAAAACTGTTGAAAACACATTTTAATGTTGAGGTTTAGTAAAGAAGGACTAAATTGGTGATGCTGTCATGATGGTTAATAACATTTATGATGGTTAGCAAAACGAAt                                                                                         >  1:83040/1‑141 (MQ=255)
                                                          aGGCATTATACGGTTCAGTCTTTAAAACTGTTGAAAACACATTTTAATGTTGAGGTTTAGTAAAGAAGGACTAAATTGGTGATGCTGTCATGATGGTTAATAACATTTATGATGGTTAGCAAAACGAATTAGAAGATCGaa                                                                             >  2:78981/1‑141 (MQ=255)
                                                                    cGGGTCAGTCTTTAAAACTGTTGAAAACACATTTTAATGTTGAGGTTTAGTAAAGAAGGACTAAATTGGTGAAGCTGTCATGATGGTTAACAACATTTATGATGGTTAGCAAAACGAATTAGAAGATCGAAAGTATACGTa                                                                   >  2:390627/1‑141 (MQ=255)
                                                                       ttCAGTCTTTAAAACTGTTGAAAACACATTTTAATGTTGAGGTTTAGTAAAGAAGGACTAAATTGGTGATGCTGTCATGATGGTTAATAACATTTATGATGGTTAGCAAAACGAATTAGAAGATCGAAAGTATACGTaaaa                                                                <  1:15675/141‑1 (MQ=255)
                                                                       ttCAGTCTTTAAAACTGTTGAAAACACATTTTAATGTTGAGGTTTAGTAAAGAAGGACTAAATTGGTGATGCTGTCATGATGGTTAATAACATTTATGATGGTTAGCAAAACGAATTAGAAGATCGAAAGTATACGTaaaa                                                                >  2:280143/1‑141 (MQ=255)
                                                                       ttCAGTCTTTAAAACTGTTGAAAACACATTTTAATGTTGAGGTTTAGTAAAGAAGGACTAAATTGGTGATGCTGTCATGATGGTTAATAACATTTATGATGGTTAGAAAAACGAATTAGAAGATCGAAAGTATACGTaaaa                                                                >  1:188354/1‑141 (MQ=255)
                                                                                  aaaCTGTTGAAAACACATTTTAATGTTGAGGTTTAGTAAAGAAGGACTAAATTGATGATGCTGTCATGATGGTTAATAACATTTATGATGGTTAGCAAAACGAATTAGAAGATCGAAAGTATACGTAAAAAATATGaaaaa                                                     >  1:260140/1‑141 (MQ=255)
                                                                                          gAAAACACATTTTAATGTTGAGGTTTAGTAAAGAAGGACTAAATTGGTGATGCTGTCATGATGGTTAATAACATTTATGATGGTTAGCAAAACGAATTAGAAGATCGAAAGTATACGTAAAAAATATGAAAAATCACGCTa                                             <  2:311512/141‑1 (MQ=255)
                                                                                            aaaCACATTTTAATGTTGAGGTTTAGTAAAGAAGGACTAAATTGGTGATGCTGTCATGATGGTTAATAACATTTATGATGGTTAGCAAAACGAATTAGAAGATCGAAAGTATACGTAAAAAATATGAAAAATCACGCTATc                                           >  1:232354/1‑141 (MQ=255)
                                                                                                     ttAATGTTGAGGTTTAGTAAAGAAGGACTAAATTGGTGATGCTGTCATGATGGTTAATAACATTTATGATGGTTAGCAAAACGAATTAGAAGATCGAAAGTATACGTAAAAAATATGAAAAATCACGCTATCATTGCACTg                                  <  2:50320/141‑1 (MQ=255)
                                                                                                     ttAATGTTGAGGTTTAGTAAAGAAGGACTAAATTGGTGATGCTGTCATGACGGTTAATAACATTTATGTTGGTTAGCAAAACGAATTAGAAGATCGAAAGTATACGTAAAAAATATGAAAAATCACGCTATCATTGCACTg                                  <  2:391143/141‑1 (MQ=255)
                                                                                                         tGTTGAGGTTTAGTAAAGAAGGACTAGATTGGTGATGCTGTCATGATGGTTAATAACATTTATGATGGTTAGCAAAACGAATTAGAAGATCGAAAGTATACGTAAAAAATATGAAAAATCACGCTATCATTGCACTGAATg                              <  2:32151/141‑1 (MQ=255)
                                                                                                          gTTGAGGTTTAGTAAAGAAGGACTAAATTGGTGATGCTGTCAGGATGGTTAATAACATTTATGATGGTTAGCAAAACGAATTAGAAGATCGAAAGTATAAAAAAAAAAAAAAAAAAATCACGCTATCATTGCACTGAATGt                             <  2:257954/141‑1 (MQ=255)
                                                                                                           ttGAGGTTTAGTAAAGAAGGACTAAATTGGTGATGCTGTCATGATGGTTAATAACATTTATGATGGTTAGAAAAAAGAATTAGAAGATAGAAAGTATACGTAAAAAAAATAAAAAATGACGCTCTCATTAAACTGAATGtt                            >  1:363312/1‑141 (MQ=255)
                                                                                                              aGGTTTAGTAAAGAAGGACTAAATTGGTGATGCTGTCATGATGGTTAATAACATTTATGATGGTTAGCAAAACGAATTAGAAGATCGAAAGTATACGTAAAAAATATGAAAAATCACGCTATCATTGCACTGAATGTTAGc                         <  2:279481/141‑1 (MQ=255)
                                                                                                                gTTTAGTAAAGAAGGACTAAATTGGTGATGCTGTCATGATGGTTAATAACATTTATGATGGTTAGCAAAACGAATTAGAAGATCGAAAGTATACGTAAAAAATATGAAAAATCACGCTATCATTGCACTGAATGTTAGCGt                       <  1:375744/141‑1 (MQ=255)
                                                                                                                    aGTAAAGAAGGACTAAATTGGTGATGCTGTCATGATGGTTAATAACATTTATGATGGTTAGCAAAACGAATTAGAAGATCGAAAGTATACGTAAAAAATATGAAAAATCACGCTATCATTGCACTGAATGTTAGCGTGAtt                   >  2:175714/1‑141 (MQ=255)
                                                                                                                             ggACTAAATTGGTGATGCTGTCATGATGGTTAATAACATTTATGATGGTTAGCAAAACGAATTAGAAGATCGAAAGTATACGTAAAAAATATGAAAAATCACGCTATCATTGCACTGAATGTTAGCGTGATTTTTatatat          >  2:406858/1‑141 (MQ=255)
                                                                                                                                 tAAATTGGTGATGCTGTCATGATGGTTAATAACATTTATGATGGTTAGCAAAACGAATTAGAAGATCGAAAGTATACGTAAAAAATATGAAAAATCACGCTATCATTGCACTGAATGTTAGCGTGATTTTTATATattaat      <  2:195241/141‑1 (MQ=255)
                                                                                                                                     ttGGTGATGCTGTCATGATGGTTAATAACATTTATGATGGTTAGCAAAACGAATTAGAAGATCGAAAGTATACGTAAAAAATATGAAAAATCACGCTATCATTGCACTGAATGTTAGCGTGATTTTTATATATTGATTAAg  <  2:110355/141‑1 (MQ=255)
                                                                                                                                     ttGGTGATGCTGTCATGATGGTTAATAACATTTATGATGGTTAGCAAAACGAATTAGAAGATCGAAAGTATACGTAAAAAATATGAAAAATCACGCTATCATTGCACTGAATGTTAGCGTGATTTTTATATATTAATTAAg  <  2:451723/141‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                         
TAATACATTTATCGATATGCCAGAAACTTATACTGACTTTAAAATGGGGTTACAAATTAGGCATTATACGGTTCAGTCTTTAAAACTGTTGAAAACACATTTTAATGTTGAGGTTTAGTAAAGAAGGACTAAATTGGTGATGCTGTCATGATGGTTAATAACATTTATGATGGTTAGCAAAACGAATTAGAAGATCGAAAGTATACGTAAAAAATATGAAAAATCACGCTATCATTGCACTGAATGTTAGCGTGATTTTTATATATTAATTAAG  >  CP000730/820806‑821079

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