Predicted mutation
evidence seq id position mutation annotation gene description
RA CP006881 2,198,604 2 bp→AA coding (378‑379/771 nt) U712_11015 ← Chemotaxis protein methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068812,198,6040GA98.0% 145.7 / ‑4.5 49L127L (CTG→TTA) U712_11015Chemotaxis protein methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (21/27);  total (22/27)
Fisher's exact test for biased strand distribution p-value = 4.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01
*CP0068812,198,6050GA97.9% 142.3 / ‑4.4 49I126I (ATC→ATTU712_11015Chemotaxis protein methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (21/26);  total (22/27)
Fisher's exact test for biased strand distribution p-value = 4.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.40e-01

CATAACTTCTATTTGCATTTTGAGTAAAATAACGATCCTTCACAGATGAAGGTACTTCCTGTAAAGACCTCTCCTGATAGACGCCTTTTTTCGCTTTCTCCAATGCTTTTTCATCAATATCAGTCGCCAGGATCTGAAAGCCGGGAAGACTTTTTTGCTGATCCAAAAGCATGGCAAGCGTGTAAGGCTCCTCACCCGTTGAGCAGGCGGCACTCCATATCTTTAAAGGCTTACCGTTCTTGATTAACGGCAG  >  CP006881/2198475‑2198727
                                                                                                                                 ||                                                                                                                          
cATAACTTCTATTTGCGTTTTGCGTAAAATAACGATCTTTCACAGATAGCGGCACTTCCTGTAAAGACCGCTCCTGATAAACGCCTTTTTTCGCTTTTTCCAATGCTTTTTCATCAATATCCGTCGCTAAAATCTGa                                                                                                                      <  3:331592/137‑1 (MQ=255)
 aTAACTTCTATTTGCGTTTTGCGTAAAATAACGATCTTTCACAGATAGCGGCACTTCCTGTAAAGACCGCTCCTGATAAACGCCTTTTTTCGCTTTTTCCAATGCTTTTTCATCAATATCCGTCGCTAAAATCTGATAg                                                                                                                   <  7:254042/139‑1 (MQ=255)
 aTAACTTCTATTTGCGTTATGCGTAATATAACGATCTTTCACAGATAGCGGCACTTCCTGTAAATACCGCTCCTGATAATCGCCTTTTTTCGCTTTTTCCAATGCTTTTTCATCAATATCCGTCGCTAAAATCTGATAg                                                                                                                   <  2:365270/139‑1 (MQ=255)
        cTATTTGCGTTTTGCGTAAAATAACGATCTTTCACAGATAGCGGCACTTCCTGTAAAGACCGCTCCTGATAAACGCCTTTTTTCGCTTTTTCCAATGCTTTTTCATCAATATCCGTCGCTAAAATCTGata                                                                                                                    <  4:315547/131‑3 (MQ=255)
        cTATTTGCGTTTTGCGTAAAATAACGATCTTTCACAGATAGCGGCACTTCCTGTAAAGACCGCTCCTGATAAACGCCTTTTTTCGCTTTTTCCAATGCTTTTTCATCAATATCCGTCGCTAAAATCTGata                                                                                                                    >  3:315547/1‑129 (MQ=255)
           tttgcgttttgcgtAAAATAACGATCTTTCACAGATAGCGGCACTTCCTGTAAAGACCGCTCCTGATAAACGCCTTTTTTCGCTTTTTCCAATGCTTTTTCATCAATATCCGTCGCTAAAATCTGATAGCCCGGaa                                                                                                            <  1:347301/136‑1 (MQ=255)
           tttgcgttttgcgtAAAATAACGATCTTTCACAGATAGCGGCACTTCCTGTAAAGACCGCTCCTGATAAACGCCTTTTTTCGCTTTTTCCAATGCTTTTTCATCAATATCCGTCGCTAAAATCTGATAGCCCGGaa                                                                                                            <  5:380305/136‑1 (MQ=255)
           tttgcgttttgcgtAAAATAACGATCTTTCACAGATAGCGGCACTTCCTGTAAAGACCGCTCCTGATAAACGCCTTTTTTCGCTTTTTCCAATGCTTTTTCATCAATATCCGTCGCTAAAATCTGATAGCCCGGaa                                                                                                            <  5:68061/136‑1 (MQ=255)
           tttgcgttttgcgtAAAATAACGATCTTTCACAGATAGCGGCACTTCCTGTAAAGACCGCTCCTGATAAACGCCTTTTTTCGCTTTTTCCAATGCTTTTTCATCAATATCCGTCGCTAAAATCTGATAGCCCGGaa                                                                                                            <  6:19195/136‑1 (MQ=255)
           tttgcgttttgcgtAAAATAACGATCTTTCACAGATAGCGGCACTTCCTGTAAAGACCGCTCCTGATAAACGCCTTTTTTCGCTTTTTCCAATGCTTTTTCATCAATATCCGTCGCTAAAATCTGATAGCCCGGaa                                                                                                            <  3:296291/136‑1 (MQ=255)
           tttgcgttttgcgtAAAATAACGATCTTTCACAGATAGCGGCACTTCCTGTAAAGACCGCTCCTGATAAACGCCTTTTTTCGCTTTTTCCAATGCTTTTTCATCAATATCCGTCGCTAAAATCTGATAGCCCGGAAGAc                                                                                                         >  2:51100/1‑139 (MQ=255)
              gcgttttGCGTAAAATAACGATCTTTCACAGATAGCGGCACTTCCTGTAAAGACCGCTCCTGATAAACGCCTTTTTTCGCTTTTTCCAATGCTTTTTCATCAATATCCGTCGCTAAAATCTGATAGCCCGGAAGAcct                                                                                                       <  4:98203/135‑3 (MQ=255)
                gttttGCGTAAAATAACGATCTTTCACAGATAGCGGCACTTCCTGTAAAGACCGCTCCTGATAAACGCCTTTTTTCGCTTTTTCCAATGCTTTTTCATCAATATCCGTCGCTAAAATCTGATAGCCCGGAAGACctttt                                                                                                    >  7:140764/2‑139 (MQ=255)
                    ttgcgaaaaTAACGATCATTCACAGATAGCGGCACTTCCTGTAAAGACCGCTCCTGATAAACGCCTTTTTTCGCTTTATCCAATGCTTTTTCATCAATATCCGTCGCTAAAATCTGATAGCCCGGAAGACCTTTTTGc                                                                                                 <  4:276149/133‑1 (MQ=255)
                     gcgTAAAATAACGATCTTTCACAGATAGCGGCACTTCCTGTAAAGACCGCTCCTGATAAACGCCTTTTTTCGCTTTTTCCAATGCTTTTTCATCAATATCCGTCGCTAAAATCTGATAGCCCGGAAGACCTTTTTGCt                                                                                                <  7:4063/136‑1 (MQ=255)
                     gcgTAAAATAACGATCTTTCACAGATAGCGGCACTTCCTGTAAAGACCGCTCCTGATAAACGCCTTTTTTCGCTTTTTCCAATGCTTTTTCATCAATATCCGTCGCTAAAATCTGATAGCCCGGAAGACCTTTTTGCt                                                                                                <  8:140764/136‑1 (MQ=255)
                                cGATCTTTCACAGATAGCGGCACTTCCTGTAAAGACCGCTCCTGATAAACGCCTTTTTTCGCTTTTTCCAATGCTTTTTCATCAATATCCGTCGCTAAAATCTGATAGCCCGGAAGACCTTTTTGCTGGTCCAAAAGCa                                                                                    >  7:123313/1‑139 (MQ=255)
                                cGATCTTTCACAGATAGCGGCACTTCCTGTAAAGACCGCTCCTGATAAACGCCTTTTTTCGCTTTTTCCAATGCTTTTTCATCAATATCCGTCGCTAAAATCTGATAGCCCGGAAGACCTTTTTGCTGGTCCAAAAGCa                                                                                    >  2:249941/1‑139 (MQ=255)
                                  aTCTTTCACAGATAGCGGCACTTCCTGTAAAGACCGCTCCTGATAAACGCCTTTTTTCGCTTTTTCCAATGCTTTTTCATCAATATCCGTCGCTAAAATCTGATAGCCCGGAAGACCTTTTTGCTGGTCCAAAAGCATg                                                                                  >  5:150741/1‑139 (MQ=255)
                                     tttCACAGATAGCGGCACTTCCTGTAAAGACCGCTCCTGATAAACGCCTTTTTTCGCTTTTTCCAATGCTTTTTCATCAATATCCGTCGCTAAAATCTGATAGCCCGGATGACCTTTTTGCTGGTCCAAAAGCATGGc                                                                                >  6:413373/2‑138 (MQ=255)
                                     tttCACAGATAGCGGCACTTCCTGTAAAGACCGCTCCTGATAAACGCCTTTTTTCGCTTTTTCCAATGCTTTTTCATCAATATCCGTCGCTAAAATCTGATAGCCCGGAAGACCTTTTTGCTGGTCCAAAAGCATGGc                                                                                >  4:438552/2‑138 (MQ=255)
                                      ttCACAGATAGCGGCACTTCCTGTAAAGACCGCTCCTGATAAACGCCTTTTTTCGCTTTTTCCAATGCTTTTTCATCAATATCCGTCGCTAAAATCTGATAGCCCGGAAGACCTTTTTGCTGGTCCAAAAGCATGGc                                                                                >  4:128305/1‑137 (MQ=255)
                                       tCACAGATAGCGGCACTTCCTGTAAAGACCGCTCCTGATAAACGCCTTTTTTCGCTTTTTCCAATGCTTTTTCATCAATATCCGTCGCTAAAATCTGATAGCCCGGAAGACCTTTTTGCTGGTCCAAAAGCATgg                                                                                 <  2:321569/135‑1 (MQ=255)
                                       tCACAGATAGCGGCACTTCCTGTAAAGACCGCTCCTGATAAACGCCTTTTTTCGCTTTTTCCAATGCTTTTTCATCAATATCCGTCGCTAAAATCTGATAGCCCGGAAGACCTTTTTGCTGGTCCAAAAGCATGGCaa                                                                              <  7:169309/138‑1 (MQ=255)
                                       tCACAGATAGCGGCACTTCCTGTAAAGACCGCTCCTGATAAACGCCTTTTTTCGCTTTTTCCAATGCTTTTTCATCAATATCCGTCGCTAAAATCTGATAGCCCGGAAGACCTTTTTGCTGGTCCAAAAGCATGGCaa                                                                              <  4:60538/138‑1 (MQ=255)
                                       tCACAGATAGCGGCACTTCCTGTAAAGACCGCTCCTGATAAACGCCTTTTTTCGCTTTTTCCAATGCTTTTTCATCAATATCCGTCGCTAAAATCTGATAGCCCGGAAGACCTTTTTGCTGGTCCAAAAGCATGGCaa                                                                              <  4:65381/138‑1 (MQ=255)
                                       tCACAGATAGCGGCACTTCCTGTAAAGACCGCTCCTGATAAACGCCTTTTTTCGCTTTTTCCAATGCTTTTTCATCAATATCCGTCGCTAAAATCTGATAGCCAGGAAGACCTTTTTGCTGGTCCAAAAGCATGGCaa                                                                              <  3:72123/138‑1 (MQ=255)
                                                      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)
                                                                                 cGCCTTTTTTCGCTTTTTCCAATGCTTTTTCATCAATATCCGTCGCTAAAATCTGATAGCCCGGAAGACCTTTTTGCTGGTCCAAAAGCATGGCa                                                                               <  8:212700/95‑1 (MQ=255)
                                                                                 cGCCTTTTTTCGCTTTTTCCAATGCTTTTTCATCAATATCCGTCGCTAAAATCTGATAGCCCGGAAGACCTTTTTGCTGGTCCAAAAGCATGGCa                                                                               >  7:212700/1‑95 (MQ=255)
                                                                                 cGCCTTTTTTCGCTTTTTCCAATGCTTTTTCATCAATATCCGTCGCTAAAATCTGATAGCCCGGAAGACCTTTTTGCTGGTCCAAAAGCATGGCAAGCGTGTATGGCTCCTCCCCCGTTGAGCAGGCGGCACTCCaaat                                   >  6:355671/1‑139 (MQ=255)
                                                                                 cGCCTTTTTTCGCTTTTTCCAATGCTTTTTCATCAATATCCGTCGCTAAAATCTGATAGCCCGGAAGACCTTTTTGCTGGTCCAAAAGCATGGCAAGCGTGTATGGCTCCTCCCCCGTTGAGCAGGCGGCACTCCaaat                                   >  2:283264/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)
                                                                                                                                 ||                                                                                                                          
CATAACTTCTATTTGCATTTTGAGTAAAATAACGATCCTTCACAGATGAAGGTACTTCCTGTAAAGACCTCTCCTGATAGACGCCTTTTTTCGCTTTCTCCAATGCTTTTTCATCAATATCAGTCGCCAGGATCTGAAAGCCGGGAAGACTTTTTTGCTGATCCAAAAGCATGGCAAGCGTGTAAGGCTCCTCACCCGTTGAGCAGGCGGCACTCCATATCTTTAAAGGCTTACCGTTCTTGATTAACGGCAG  >  CP006881/2198475‑2198727

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