Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 2,343,788 A→C 27.4% L253R (CTT→CGT)  U712_11850 ← Farnesyl diphosphate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068812,343,7880AC27.4% 58.6 / 3.4 31L253R (CTT→CGT) U712_11850Farnesyl diphosphate synthase
Reads supporting (aligned to +/- strand):  ref base A (20/2);  new base C (2/7);  total (22/9)
Fisher's exact test for biased strand distribution p-value = 4.22e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTTAGTGATCTCTTGCCGCAATTAAATCACAAAGCTCATATAAAAGGTCTTTTTGAAGAGAGAGTCCGCCGATCAATCGCTTCGCCTCTTTTATATGAACATCCAATTTATGTTTGGCCCCTTCCAATGAAAGAAGCGACGGGTATGTCGATTTGTCATTTGTGGTATCCGAGCCGACACGTTTGCCGATTTTCTCTTCACTGCCTTCTAAATCTAAAATATCGTCTCTGATTTGAAATCCGATTCCAATGTGGCTGCTGAAGGTACGCA  >  CP006881/2343653‑2343923
                                                                                                                                       |                                                                                                                                       
ttttAGTGATCTCTTGCCGCAATTAAATCACAAAGCTCATATAAAAGGTCTTTTTGAAGAGAGAGTCCGCCGATCAATCGCTTCGCCTCTTTTATATGAACATCCAATTTATTTTTGGACCCTTCCAATGAAAGACGc                                                                                                                                       >  5:31518/1‑138 (MQ=255)
          ctctTGCCGCAATTAAATCACAAAGCTCATATAAAAGGTCTTTTTGAAGAGAGAGTCCGCCGATCAATCGCTTCGCCTCTTTTATATGAACATCCAATTTATTTTTGGCCCCTTCCAATGAAAGAAGCGACGGGTATGt                                                                                                                            >  1:378645/1‑139 (MQ=255)
                cccaaaTTAAAACAAAACGATAATATATATGGTTTTTTTTACGAGCGAGGCCCCCTATAAATCGCTGCGCCTATTTTATATCAACATCAAATTTTTGTTTGCCCCATACCAATCAAACACGCGACGGGTATGTCGAttt                                                                                                                      <  2:242159/135‑1 (MQ=255)
                     aTTAAATCACAAAGCTCATATAAAAGGTCTTTTTGAAGAGAGAGTCCGCCGATCAATCGCTTCGCCTCTTTTATATGAACATCCAATTTATGTTTGGCCCCTTCCAATGAAAGAAGCGACGGGTATGTCGATTTGTCAt                                                                                                                 >  2:279750/1‑139 (MQ=255)
                     aTTAAATCACAAAGCTCATATAAAAGGTCTTTTTGAAGAGAGAGTCCGCCGATCAATCGCTTCGCCTCTTTTATATGAACATCCAATTTATGTTTGGCCCCTACCAATGAAAGAAGCGACGGGTATGTCGATTTGTCAt                                                                                                                 >  3:170638/1‑139 (MQ=255)
                              caAAGCTCATATAAAAGGTCTTTTTGAAGAGAGAGTCCGCCGATCAATCGCTTCGCCGCTTTTATATGAACATCCAATTTATGTTTGGCCCCTTCCAATGAAAGAAGCGACGGGTATGTCGATTTGTCATTTGTGGTAt                                                                                                        >  5:366164/1‑139 (MQ=255)
                                      atatAAAAGGTCTTTTTGAAGAGAGAGTCCGCCGATCAATCGCTTCGCCTCTTTTATATGAACATCCAATTTATGTTTGGCCCCTTCCAATGAAAGAAGCGACGGGTATGTCGATTTGTCATTTGTGGTATCCGAGcc                                                                                                 >  3:343493/1‑138 (MQ=255)
                                             aGGTCTTTTTGAAGAGAGAGTCCGCCGATCAATCGCTTCGCCTCTTTTATATGAACATCCAATTTATGTTTGGCCCCTTCCAATGAAAGAAGCGACGGGTATGTCGATTTGTCATTTGTGGTATCCGAGCCGCCACGtt                                                                                         >  7:322813/1‑139 (MQ=255)
                                                 cTTTTTGAAGAGAGAGTCCGCCGATCAATCGCTTCGCCTCTTTTATATGAACATCCAATTTATGTTTGGCCCCTTCCAATGAAAGAAGCGACGGGTATGTCGATTTGTCATTTGTGGTATCCGAGCCGACACGTTTGcc                                                                                     >  1:354529/1‑139 (MQ=255)
                                                        ccgcgcgagCCCGCCGACCAAGCGCTCCGCCTATTTTATATGAACATCCAATTTTTGTTGGCCCCCAACAAATCAACCACGCGACGGGTATGTCGATTTGTCATTTGTGGTATCCGAGCCGACACGTTTGCCGAtttt                                                                               <  1:80095/132‑1 (MQ=255)
                                                             agagTCCGCCGATCAATCGCTTCGCCGCTGTTATATGAACATCCAATTTATTTTTGGCCCCTTCCAATGAAAGCCGCGACGGGGATGTCGATTTGTCATTTGTGGTATACGAGCCGCCCCGTTTGCCGCTTTTCTCtt                                                                          >  1:334928/1‑138 (MQ=255)
                                                              gagCCCGCCGAGCAACCGCTACGCCTCTTTAATATGAACATCCAATTTATGTTGGGCCCCTCCCAAGAAAACACGCGACGGGTATGTCGATTTGTCATTTGTGGTATCCGAGCCGACACGTTTGCCGATTTTCTCTTCa                                                                        <  6:163905/139‑1 (MQ=255)
                                                                         acatcggCTGCGCCTTTTTTATATTAACATCCTTTTTATTTTGGGCCCCTTAAATTTAAATACGCGACGGGTATGTCGATTTGTCATTTGTGGTATCCGAGCCGACACGTTTGCCGATTTTCTCTTCACTGCCTtctaa                                                             <  6:366164/133‑1 (MQ=255)
                                                                          cAATCGCTTCGCCTGTTTTATATGAACATCCAATTTATGTTTGGCCCCTTCCAAGGAAAGAAGCGACGGGTATGTCGATTTGTCATTTGTGGTTTCCGAGCCGACACGTTTGCCGATTTTCTCTTCACTGCCTtctaaa                                                            >  4:162280/1‑139 (MQ=255)
                                                                                cTTCGCCTCTTTTATATGAACATCCAATTTATGTTTGGCCCCTTCCAATGAAAGAAGCGACGGGTATGTCGATTTGTCATTTGTGGTATCCGAGCCGACACGTTTGCCGATTTTCTCTTCACTGCCTtctaaatctaa                                                       >  4:184576/1‑138 (MQ=255)
                                                                                  tCGCCTCTTTTATATGAACATCCAATTTATGTTTGGCCCCTTCCAATGAAAGAAGCGACGGGTATGTCGATTTGTCATTTGTGGTATCCGAGCCGACACGTTTGCCGATTTTCTCTTCACTGCCTTCTAAATCTaaaa                                                     >  6:149126/1‑138 (MQ=255)
                                                                                           ttATATGAACATCCAATTTATGTTTGGCCCCTTGCAATGAAAGAAGAGACGGGTATGTCGATTTGGAATTTGTGGTATACGAGCCGACACGATTTCCGATTTTCTCTTCACTGCCTTCTAAATCTAAAATATCGActct                                           >  2:180475/1‑139 (MQ=255)
                                                                                                   aCATCCAATTTATGTTTGGCCCCTTCCAATGAAAGAAGCGACGGGTATGTCGATTTGTCATTTGTGGTATCCGAGCCGACACGTTTGCCGATTTTCTCTTCACTGCCTTCTAAATCTAAAATATCGTCTCTGATTTGa                                    >  6:105759/1‑138 (MQ=255)
                                                                                                   aCATCCAATTTATGTTTGGCCCCTTCCAATGAAAGAAGCGACGGGTATGTCGATTTGTCATTTGTGGTATCCGAGCCGACACGTTTGCCGATTTTCTCTTCACGGCCTTCTAAATCTAAAATATAGTCTCTGATTTGa                                    >  1:194062/1‑138 (MQ=255)
                                                                                                   aCATCCAATTTATGTTTGGCCCCTGCCAATGAAAGAAGCGACGGGTATGTCGATTTGTCATTTGTGGTATCCGAGCCGCGACGTTTGCCGATTTTGACGTCACTGCCGTCGAAATCTAAAATTAAGGCTCTTATTTGa                                    >  1:192679/1‑138 (MQ=255)
                                                                                                    cATCAAATTAATGTTTGCCCCCATACAATGAACGACGCGACGGGTATGTCGATTTGTCATTTGTGGTATCCGAGCCGACACGTTTGCCGATTTTCTCTTCACTGCCTTCTAAATCTAAAATATCGTCTCTGATTTGaaa                                  <  5:335696/139‑1 (MQ=255)
                                                                                                               tGTTTGGCCCCTTCCAATGAAAGAAGCGACGGGTATGTCGATTTGTCATTTGTGGTATCCGAGCCGACACGTTTGCCGATTTTCTCTTCACTGCCTTCGAAATCTAAAATATCCTCTCTGATTTGATATCCGAttgaaa                       >  8:237466/1‑135 (MQ=255)
                                                                                                                   tGGCCCCTTCCAATGAAAGAAGCGACGGGTATGTCGATTTGTCATTTGTGGTATCCGAGCCGACACGTTTGCCGATTTTCTCTTCACTGCCTTCTAAATCTAAAATATCGTCTCTGATTTGAAATCCGATTCCAAtgtg                   >  5:215383/1‑139 (MQ=255)
                                                                                                                    ggCCCCTTCCAATGAAAGAAGCGACGGGTATGTCGATTTGTCATTTGTGGTATCCGAGCCGACACGTTTGCCGATTTTCTCTTCACTGCCTTCTAAATCTAAAATATCGTCTCTGATTTGAAATCCGATTCCAATGTgg                  >  3:111968/1‑139 (MQ=255)
                                                                                                                    ggCCCATTCCAATGAAAGAAGAGACGGGTATGTCGATTTGGCATTTGTGGTATCCGAGCCGACACGTTTGCAGATTTTCTCTTCTATGCCTTCTAAATCTAAAATATAGTCTCTTATTTGATATCATATTACAAttagg                  >  5:137388/1‑135 (MQ=255)
                                                                                                                     gCCCCTTCCAATGAAAGAAGCGACGGGTATGTCGATTTGTCATTTGTGGTATACGAGCCGACACGTTTGCCGATTTTCTCTTCACTGCCTTCTAAATCTAAAATATCGTCTCTGATTTGAAATCCGATTCCAATGTGgc                 >  5:129778/1‑139 (MQ=255)
                                                                                                                        ccTTCCAATCATACAAGCGACGGGTATGTCGATTTGTCATTTGTGGTATCCGAGCCGACACGTTTGCCGATTTTCTCTTCACTGCCTTCTAAATCTAAAATATCGTCTCTGATTTGAAATCCGATTCCAATGTGgctgc              <  4:170638/139‑1 (MQ=255)
                                                                                                                        cataccAAAGAAACAAGCGACGGGTATGTCGATTTGTCATTTGTGGTATCCGAGCCGACACGTTTGCCGATTTTCTCTTCACTGCCTTCTAAATCTAAAATATCGTCTCTGATTTGAAATCCGATTCCAATGTGgctgc              <  4:343493/135‑1 (MQ=255)
                                                                                                                          ttCCAATGAAAGAAGCGACGGGTATGTCGATTTGTCATTTGTGGTATCCGAGCCGACACGTTTGCCGATTTTCTCTTCACTGCCTTCTAAATCTAAAATATCGTCTCTGATTTGAAATCCGATTCCAATGTGgctgct             >  1:320595/1‑138 (MQ=255)
                                                                                                                            ccAATGAACCACGCGACGGGTATGTCGATTTGTCATTTGTGGTATCCGAGCCGACACGTTTGCCGATTTTCTCTTCACTGCCTTCTAAATCTAAAATATCGTCTCTGATTTGAAATCCGATTCCAATGTGGCTGCTGaa          <  5:315731/139‑1 (MQ=255)
                                                                                                                                  aaaGACGCGACGGGTATGTCGATTTGTCATTTGTGGTATCCGAGCCGACACGTTTGCCGATTTTCTCTTCACTGCCTTCTAAATCTAAAATATCGTCTCTGATTTGAAATCCGATTCCAATGTGGCTGCTGAAGGTACg    <  5:105759/139‑1 (MQ=255)
                                                                                                                                   aacaagCGACGGGTATGTCGATTTGTCATTTGTGGTATCCGAGCCGACACGTTTGCCGATTTTCTCTTCACTGCCTTCTAAATCTAAAATATCGTCTCTGATTTGAAATCCGATTCCAATGTGGCTGCTGAAGGTACGc   <  6:69802/136‑1 (MQ=255)
                                                                                                                                    agaagCGACGGGTATGTCGATTTGTCATTTGTGGTATCCGAGCCGACCCGTTTGCCGATTTTCTCTTCACTGCCTTCTAAATCTAAAATATCGTCTATGATTTGAAATCCGATTCCAATGTGGCTGCTGAAGGTACg    >  6:15421/1‑137 (MQ=255)
                                                                                                                                    agaagCGACGGGTATGTCGATTTGTCATTTGTGGTATCCGAGCCGACACGTTTGCCGATTTTCTCTTCACTGCCTTCTAAATCTAAAATATCGTCTCTGATTTGAAATCCGATTCCAATGTGGCTGCTGAAGGTACGCa  <  3:184576/139‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                       
TTTTAGTGATCTCTTGCCGCAATTAAATCACAAAGCTCATATAAAAGGTCTTTTTGAAGAGAGAGTCCGCCGATCAATCGCTTCGCCTCTTTTATATGAACATCCAATTTATGTTTGGCCCCTTCCAATGAAAGAAGCGACGGGTATGTCGATTTGTCATTTGTGGTATCCGAGCCGACACGTTTGCCGATTTTCTCTTCACTGCCTTCTAAATCTAAAATATCGTCTCTGATTTGAAATCCGATTCCAATGTGGCTGCTGAAGGTACGCA  >  CP006881/2343653‑2343923

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