Predicted mutation
evidence seq id position mutation annotation gene description
RA CP006881 2,198,659 A→T P108P (CCT→CCA)  U712_11015 ← Chemotaxis protein methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068812,198,6590AT93.8% 88.0 / ‑1.8 32P108P (CCT→CCA) U712_11015Chemotaxis protein methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base T (18/12);  total (19/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-01

CTTCCTGTAAAGACCTCTCCTGATAGACGCCTTTTTTCGCTTTCTCCAATGCTTTTTCATCAATATCAGTCGCCAGGATCTGAAAGCCGGGAAGACTTTTTTGCTGATCCAAAAGCATGGCAAGCGTGTAAGGCTCCTCACCCGTTGAGCAGGCGGCACTCCATATCTTTAAAGGCTTACCGTTCTTGATTAACGGCAGAATAGCAGTCTCAAGAACTTCCCATCTTTTATAATTACGATAAAATTCTGAAACATTGATCGTCATT  >  CP006881/2198529‑2198794
                                                                                                                                  |                                                                                                                                       
cTTCCTGTAAAGACCGCTCCTGATAAACGCCTTTTTTCGCTTTTTCCAATGCTTTTTCATCAATATCCGTCGCTAAAATCTGATAGCCCGGAAGACCTTTTTGCTGGTCCAAAAGCATGGCAAGCGTGTATGGctcctc                                                                                                                                 <  7:372215/139‑1 (MQ=255)
 ttCCTGTAAAGACCGCTCCTGATAAACGCCTTTTTTCGCTTTTTCCAATGCTTTTTCATCAATATCCGTCGCTAAAATCTGATAGCCCGGAAGACCTTTTTGCTGGTCCAAAAGCATGGCAAGCGTGTATGGctcctcc                                                                                                                                >  5:181863/1‑138 (MQ=255)
   ccTGTAAAGACCGCTCCTGATAAACGCCTTTTTTCGCTTTTTCCAATGCTTTTTCATCAATATCCGTCGCTAAAATCTGATAGCCCGGAAGACCTTTTTGCTGGTCCAAAAGCATGGCAAGCGTGTATGGCTCCTcccc                                                                                                                              >  5:389967/1‑139 (MQ=255)
           gACCGCTCCTGATAAACGCCTTTTTTCGCTTTTTCCAATGCTTTTTCATCAATATCCGTCGCTAAAATCTGATAGCCCGGAAGACCTTTTTGCTGGTCCAAAAGCATGGCAAGCGTGTATGGCTCCTCCCCCGTTGAGc                                                                                                                      >  3:112376/1‑139 (MQ=255)
            aCCGCTCCTGATAAACGCCTTTTTTCGCTTTTTCCAATGCTTTTTCATCAATATCCGTCGCTAAAATCTGATAGCCCGGAAGACCTTTTTGCTGGTCCAAAAGCATGGCAAGCGTGTATGACTCCTCCCCCGTTgcg                                                                                                                       >  2:184476/1‑135 (MQ=255)
              cgctcCTGATAAACGCCTTTTTTCGCTTTTTCCAATGCTTTTTCATCAATATCCGTCGCTAAAATCTGATAGCCCGGAAGACCTTTTTGCTGGTCCAAAAGCATGGCAAGCGTGTATGGCTCCTCCCCCGTTGAGCAgg                                                                                                                   >  7:303128/3‑139 (MQ=255)
                           cGCCTTTTTTCGCTTTTTCCAATGCTTTTTCATCAATATCCGTCGCTAAAATCTGATAGCCCGGAAGACCTTTTTGCTGGTCCAAAAGCATGGCAAGCGTGTATGGCTCCTCCCCCGTTGAGCAGGCGGCACTCCaaat                                                                                                      >  2:283264/1‑139 (MQ=255)
                           cGCCTTTTTTCGCTTTTTCCAATGCTTTTTCATCAATATCCGTCGCTAAAATCTGATAGCCCGGAAGACCTTTTTGCTGGTCCAAAAGCATGGCAAGCGTGTATGGCTCCTCCCCCGTTGAGCAGGCGGCACTCCaaat                                                                                                      >  6:355671/1‑139 (MQ=255)
                             ccTTTTTTCGCTTTTTCCAATGCTTTTTCATCAATATCCGTCGCTAAAATCTGATAGCCCGGAAGACCTTTTTGCTGGTCCAAAATCAAGGCAAGCGTGTATGGCTCCTCCCCCGTTGCGCAGGCGGCCCTCCaaatat                                                                                                    >  2:441274/1‑137 (MQ=255)
                             ccTTTTTTCGCTTTTTCCAATGCTTTTTCATCAATATCCGTCGCTAAAATCTGATAGCCCGGAAGACCTTTTTGCTGGTCCAAAAGCATGGCAAGCGTGTATGGCTCCTCCCCCGTTGAGCAGGCGGCACTCCaaattt                                                                                                    >  6:223418/1‑137 (MQ=255)
                                 ttttCGCTTTTTCCAATGCTTTTTCATCAATATCCGTCGCTAAAATCTGATAGCCCGGAAGACCTTTTTGCTGGTCCAAAAGCATGGCAAGCGTGTATGGCTCCTCCCCCGTTGAGCAGGCGGCACTCCAAATTTTTa                                                                                                 <  6:150741/138‑1 (MQ=255)
                                 ttttCGCTTTTTCCAATGCTTTTTCATCAATATCCGTCGCTAAAATCTGATAGCCCGGAAGACCTTTTTGCTGGTCCAAAAGCATGGCAAGCGTGTATGGCTCCTCCCCCGTTGAGCAGGCGGCACTCCAAATTTTTa                                                                                                 <  6:389967/138‑1 (MQ=255)
                                  tttCGCTTTTTCCAATGCTTTTTCATCAATATCCGTCGCTAAAATCTGATAGCCCGGAAGACCTTTTTGCTGGTCCAAAAGCATGGCAAGCGTGTATGGCTCCTCCCCCGTTGAGCAGGCGGCACTCCAAATTTTTaa                                                                                                <  1:441274/138‑1 (MQ=255)
                                         ttttcCAATGCTTTTTCATCAATATCCGTCGCTAATATCAGATAGCCCGGAAGACCTTTTTGCGGGTCCAATAGCATGGCAAGCGTGTATGGCTCCTCCCCCGTTGAGCAGGCGGCACTCCAAATTTTTAAAGGCCTag                                                                                        <  8:314438/136‑2 (MQ=255)
                                          tttcCAATGCTTTTTCATCAATATCCGTCGCTAAAATCTGATAGCCCGGAAGACCTTTTTGCTGGTCCAAAAGCATGGCAAGCGTGTATGGCTCCTCCCCCGTTGAGCAGGCGGCACTCCAAATTTTTAAAGGCCTag                                                                                        <  8:428353/136‑2 (MQ=255)
                                           ttcCAATGCTTTTTCATCAATATCCGTCGCTAAAATCTGATAGCCCGGAAGACCTTTTTGCTGGTCCAAAAGCATGGCAAGCGTGTATGGctc                                                                                                                                    >  2:428877/2‑93 (MQ=255)
                                           ttcCAATGCTTTTTCATCAATATCCGTCGCTAAAATCTGATAGCCCGGAAGACCTTTTTGCTGGTCCAAAAGCATGGCAAGCGTGTATGGctc                                                                                                                                    <  1:428877/92‑1 (MQ=255)
                                                      tttCATCAATATCCGTCGCTAAAATCTGATAGCCCGGAAGACCTTTTTGCTGGTCCAAAAGCAGGGCAAGCGTGTATGGCTCCTCCCCCGTTGAGCAGGCGGCACTCCAAATTTTTAAAGGCCTAGAGGTCTTGAtt                                                                             <  4:112376/137‑1 (MQ=255)
                                                             aaTATCCGTCGCTAAAATCTGATAGCCCGGAAGACCTTTTTGCTGGTCCAAAAGCATGGCAAGCGTGTATGGCTCCTCCCCCGTTGAGCAGGCGGCACGCCAAATTTTTTAAGGCCTAGAGGTCTTGATTAACGGCAg                                                                     >  4:1523/1‑138 (MQ=255)
                                                             aaTATCAGTCGCCAGGATCTGAAAGCCGGGAAGACTTTTTTGCTGATCCAAAAGCATGGCAAGCGTGTAAGGCTCCTCACCCGTTGAGCAGGCGGCACTCCATATCTTTAAAGGCTTACCGTTCTTGATTAACGGCAg                                                                     >  6:427556/1‑138 (MQ=255)
                                                                           aaaTCTGATAGCCCGGAAGACCTTTTTGCTGGTCCAAAAGCATGGCAAGCGTGTATGGCTCCTCCCCCGTTGAGCAGGCGGCACTCCAAATTTTTAAAGGCCTAGAGGTCTTGATTAACGGCAGAATTGCAGTCTCAAg                                                      <  6:282011/137‑1 (MQ=255)
                                                                                                 tttttGCTGGTCCAAAAGCATGGCAAGCGTGTATGGCTCCTCCCCCGTTGAGCAGGCGGCACTCCAAATTTTTAAAGGCCTAGAGGTCTTGATTAACGGCAGAATTGCAGTCTCAAGAACTTCCCATCTTTTATAATTa                                >  4:206269/1‑139 (MQ=255)
                                                                                                     tGCTGGTCCAAAAGCATGGCAAGCGTGTATGGCTCCTCCCCCGTTGAGCAGGCGGCTCTCCAAATTTTTAAAGGCCTAGAGGTCTTGATTAACGGCAGAATTGCAGTCTCAATAACTTCCCATCTTTTATTATTACGa                             >  3:11860/1‑138 (MQ=255)
                                                                                                          gtCCAAAAGCATGGCAAGCGTGTATGGCTCCTCCCCCGTTGAGCAGGCGGCACTCCAAATTTTTAAAGGCCTAGAGGTCTTGATTAACGGCAGAATTGCAGTCTCAAGAACTTCCCATCTTTTATAATTACGATAAAAt                       >  4:85657/2‑139 (MQ=255)
                                                                                                                     tGGCAAGCGTGTATGGCTCCTCCCCCGTTGAGCAGGCGGCACTCCAAATTTTTAAAGGCCTAGAGGTCTTGATTAACGGCAGAATTGCAGTCTCAAGAACTTCCCATCTTTTATAATTACGATAAAATTCTGAAACAtt            >  3:66275/1‑139 (MQ=255)
                                                                                                                     tGGCAAGCGTGTATGGCTCCTCCCCCGTTGAGCAGGCGGCACTCCAAATTTTTAAAGGCCTAGAGGTCTTGATTAACGGCAGAATTGCAGTCTCAAGAACTTCCCATCTTTTATAATTACGATAAAATTCTGAAACAtt            <  1:51100/139‑1 (MQ=255)
                                                                                                                     tGGCAAGCGTGTATGGCTCCTCCCCCGTTGAGCAGGCGGCACTCCAAATTTTTAAAGGCCTAGAGGTCTTGATTAACGGCAGAATTGCAGTCTCAAGAACTTCCCATCTTTTATAATTACGATAAAATTCTGAAACAtt            >  1:437995/1‑139 (MQ=255)
                                                                                                                     tGGCAAGCGTGTAAGGCTCCTCACCCGTTGAGCAGGCGGCACTCCATATCTTTAAAGGCTTACCGTTCTTGATTAACGGCAGAATAGCAGTCTCAAGAACTTCCCATCTTTTATAATTACGATAAAATTCTGAAACAtt            <  5:427556/139‑1 (MQ=255)
                                                                                                                            cGTGTATGGCTCCTCCCCCGTTGAGCAGGCGGCACTCCAAATTTTTAAAGGCCTAGAGGTCTTGATTAACGGCAGAATTGCAGTCTCAAGAACTTCCCATCTTTTATAATTACGATAAAATTCTGAAACATTGATCGt      >  4:177009/1‑138 (MQ=255)
                                                                                                                             gtgtATGGCTCCTCCCCCGTTGAGCAGGCGGCACTCCAAATTTTTAAAGGCCTAGAGGTCTTGATTAACGGCAGAATTGCAGTCTCAAGAACTTCCCATCTTTTATAATTACGATAAAATTCTGAAACATTGATCGTCa    <  6:71533/139‑1 (MQ=255)
                                                                                                                              tgtATGGCTCCTCCCCCGTTGAGCAGGCGGCACTCCAAATTTTTAAAGGCCTAGAGGTCTTGATTAACGGCAGAATTGCAGTCTCAAGAACTTCCCATCTTTTATAATTACGATATAATTCTGAAACATTGATCGTCAt   >  6:188295/1‑139 (MQ=255)
                                                                                                                               gtATGGCTCCTCCCCCGTTGAGCAGGCGGCACTCCAAATTTTTAAAGGCCTAGAGGTCTTGATTAACGGCAGAATTGCAGTCTCAAGAACTTCCCATCTTTTATAATTACGATAAAATTCTGAAACATTGATCGTCAtt  <  4:66275/139‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                       
CTTCCTGTAAAGACCTCTCCTGATAGACGCCTTTTTTCGCTTTCTCCAATGCTTTTTCATCAATATCAGTCGCCAGGATCTGAAAGCCGGGAAGACTTTTTTGCTGATCCAAAAGCATGGCAAGCGTGTAAGGCTCCTCACCCGTTGAGCAGGCGGCACTCCATATCTTTAAAGGCTTACCGTTCTTGATTAACGGCAGAATAGCAGTCTCAAGAACTTCCCATCTTTTATAATTACGATAAAATTCTGAAACATTGATCGTCATT  >  CP006881/2198529‑2198794

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