Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB143 1,630,540 A→C 20.6% L419R (CTG→CGG)  czcS ← Sensor protein CzcS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1431,630,5400AC20.6% 50.3 / 2.6 40L419R (CTG→CGG) czcSSensor protein CzcS
Reads supporting (aligned to +/- strand):  ref base A (16/14);  new base C (7/2);  total (24/16)
Fisher's exact test for biased strand distribution p-value = 2.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCTCCAGCCGCGAACGCCCCGAACCT‑ACCCGGTAAAACCGCTCGAACAAGCGCGACTGGTGCTCTGCAGGCACACCCGGACCCTGGTTCTCCACGGATAGAATGACTTCGCTGCCCGAGCGTTCGATGCGCAGAGCCA  >  ACB143/1630410‑1630672
                                                                                                                                  |                                                                                                                                     
cGATCACTCGGCAGGAAACTCCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCTCCAGCcgc                                                                                                                                  <  5:122892/136‑1 (MQ=255)
   tCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCTCCAGCCGCGaa                                                                                                                               >  1:32798/1‑136 (MQ=255)
     aCTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCTCCAGCCGCGAACg                                                                                                                             <  6:128692/136‑1 (MQ=255)
     aCTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCAACCCGCCCGCCATGCAACTGCATGATCGACATCACGATCGCCAACCCCAACCCGGTCCCGCCCTCCAGCCGCGAACg                                                                                                                             <  7:255358/136‑1 (MQ=255)
         ggCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCTCCAGCCGCGAACGcccc                                                                                                                         <  1:282349/136‑1 (MQ=255)
          gCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCTCCAGCCGCGAACGCCCCg                                                                                                                        <  2:301472/136‑1 (MQ=255)
             ggAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCTCCAGCCGCGAACGCCCCGAAc                                                                                                                     <  3:105591/136‑1 (MQ=255)
              gAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCACCAGCCGCGAACCCCCCCAAc                                                                                                                     >  2:177133/1‑135 (MQ=255)
                           aaaCCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCCACCCCCACCCGGGGCCGCCCGCCCGCCGCGCACGCCCCGACCCG‑ACCCGGGaaacc                                                                                                      >  2:299334/1‑134 (MQ=255)
                           aaaCCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCTCCGGCCGCGAACGCCCCGAACCT‑ACCCGGTAAAAc                                                                                                      >  6:69777/1‑136 (MQ=255)
                                gggTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCTCCCGCCGCGCACGCCCCGAACCT‑TCCCGGTAAAACCGCTc                                                                                                 >  6:288536/1‑136 (MQ=255)
                                     ggCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCTCCAGCCGCGAACGCCCCGAACCT‑ACCCCGTAAAACCGCTCGAAc                                                                                             >  3:134233/1‑135 (MQ=255)
                                              gggCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGCTCGCCAACCCCACCCCGGTCCCGCCCGCCCGCCGCGCACGCCCCGAACCT‑ACCCGGTAACACCCCTCGACCAAGCGCGCCt                                                                                   >  5:77492/1‑136 (MQ=255)
                                              gggCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATAGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCTCCCGCCGCGACCGCCCCGCACCG‑ACCCGGTTAATCCGCGAGCACCAGcgcgg                                                                                     >  6:89555/1‑133 (MQ=255)
                                               ggCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCTCCAGCCGCGAACGCCCCGAACCT‑ACCCGGTAAAACCGCTCGAACAAGCGCGACTg                                                                                  <  8:254900/136‑1 (MQ=255)
                                                    ctgctgACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCCACCCCGTCCCGCCCGCCCCCCGCCGACCCCCCCAACCG‑ACCCCGGAAAACCCCGCGCACAAGCGCGCCTGGTGct                                                                             >  2:155152/1‑136 (MQ=255)
                                                    ctgctgACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCTCCAGCCGCGAACGCCCCGAACCT‑ACCCGGTAAAACCGCTCGAACAAGCGCGACTGGTGct                                                                             >  2:23953/1‑136 (MQ=255)
                                                    ctgctgACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCTCCACCCCCGAACGCCCCCAACCT‑ACCCGGTAAAACCGCTAGAAAAAGCGCGACTGGTGc                                                                              >  3:107146/1‑135 (MQ=255)
                                                    ctgctgACCTAGACCCGCCCGCCATGCAACTGCATGATCGAATTCACGATCGCCAACCCCCACCCGGTACCGCCCTCCAGCCGCGGACGCCCCGAACCT‑CACCCGTAAAACCGCTCGAACCAGCGCGACTGGGGct                                                                             >  2:40938/1‑136 (MQ=255)
                                                                aCCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGACCCGCCCTCCAGCCGCGAACGCCCCGAACCT‑ACCCGGTAAAACCGCTCGAACAAGCGCGACTGGTGCTCTGCAGGCacac                                                                 <  8:20982/136‑1 (MQ=255)
                                                                  ccgcccgcCATGCAACTGCATGCTCGACTTCACGATCGCCAACCCCAACCCCGGCCCGCCCTCCCCCCGCGCACCCCCCCACCCTAACCCGGAAAAACC‑CTCGCACAAGCGCCCAGTGTGCTCTGCAGGCcccccc                                                               >  6:30290/1‑136 (MQ=255)
                                                                  ccgcccgcCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCCCTCCCGCCCTCCAGCCGCGAACGCCCCGAACCA‑ACCCGGTAAAACCACTCCAACAAGCCCGACTGGTGCTCTGCAGGCACAccc                                                               >  6:113139/1‑136 (MQ=255)
                                                                               aaCTGCATGATCTACTTCACCATCGCCAACCCCACCCGGGTCCCCCCCTCCAGCCGCTACCGCGCCGAACCT‑ACCGGGTAAAACCGCTCCAACAAGCGCGACTGGTGCTCTGCAGGCACACCCGGACCCTGGTtct                                                  <  1:177133/136‑1 (MQ=255)
                                                                                       gATCGACTTCACGATCGCCAACCCCAACCCCGTCCCCCCCTACACCCGCCAACGCCCCGCACCA‑ACCCGGTAAAAACCCTCGCACAAGCGCGCCTGGTGCTAGGCACGCACACCCGGACCCTGGTTGTCCACGGa                                           >  4:8907/1‑135 (MQ=255)
                                                                                                cccgcccGCCACCCCCCACCGGGCCCCGCCCGCCCGCCGCGCACGCCCCGACCTT‑CCCCGGTAAAACGGCTCGAACACGCGCGACTGGTGCTCTGCAGGCACACCCGGACCCTGGTTCTCCACGGATAGAATGa                                   <  7:130709/128‑1 (MQ=255)
                                                                                                cACGATCGCCAACCCCAACCCGGTCCCGCCCTCCAGCCGCGAACGCCCCGAACCT‑ACCCGGAAAAACCGCTCGAACACGCGCGACTGGTGCTCTGCAGGCACACCCGGACCCTGGTTCTCCACGGATAGAATGa                                   <  4:107146/134‑1 (MQ=255)
                                                                                                 aCGATCGCCAACCCCAACCCGGTCCCGCCCTCCAGCCGCGAACGCCCCGAACCT‑ACCCGGTAAAACCGCTCGAACAAGCGCGACTGGTGCTCTGCAGGCACCCCCGGACCCTGGTTCTCCACGGATAGAATGACt                                 >  6:256299/1‑135 (MQ=255)
                                                                                                   gATCGCCAACCCCAACCCGGTCCCGCCCTCCATCCGCGAACGCCCCGAACCT‑ACCCGGTAAAACCGCTCGAACAAGCGCGACTGGTGCTCTGCAg                                                                       >  5:132308/1‑95 (MQ=255)
                                                                                                   gATCGCCAACCCCAACCCGGTCCCGCCCTCCAGCCGCGAACGCCCCGAACCT‑ACCCGGTAAAACCGCTCGAACAAGCGCGACTGGTGCTCTGCAg                                                                       <  6:132308/95‑1 (MQ=255)
                                                                                                        ccaaccccaaccCGGTCCCGCCCTCCAGCCGCGAACGCCCCGAACCT‑ACCCCGTAAAACAGCTCGAACAAGCGCGACTGGTGCTCTGCAGGCACACCCGGACCCTGGTTGTCCACGGATAGAATGACTTAGCTGcc                         >  4:115585/1‑136 (MQ=255)
                                                                                                                    cGGTCCCGCCCTCCAGCCGCGCACGCCCCGAACCT‑ACCCGGTAAAACCGCTCGAACAAGCGCGACTGGTGCTCTGCAGGCACACCCGGACCCTGGTTCTCCACGGATAGAATGACTTCGCTGCCCGAGCGTTCGAt             >  5:258782/1‑136 (MQ=255)
                                                                                                                    acgTCCCGCCCTCCAGCCGCGAAAGAAAAAAAACA‑ACCCCGTAAAACCGCTCGAACAAGCGCGACTGGTGCTCTGCAGGCACACCCGGCCCCTGGTTCTCCCCGGATAGAATGACTTCGCTGCCCGAGCGTTCGAt             >  4:159605/3‑136 (MQ=255)
                                                                                                                      gTCCCGCCCTCCAGCCGCGAACGCCCCGAACCT‑ACCCGGTAAAACCGCTCGAACACGCGCGACTGGTGCTCTGCAGGCACACCCGGACCCTGGTTCTCCACGGATAGAATGACTTCGCTGCCCGAGCGTTCGATgc           <  1:72460/136‑1 (MQ=255)
                                                                                                                      gTCCCGCCCTCCAGCCGCGAACGCCCCGAACCT‑ACCCGGTAAAACCGCTCGAACAAGCGCGACTGGTGCTCTGCAGGCACACCCGGACCCTGGTTCTCCACGGATAGAATGACTTCGCTGCCCGAGCGTTCGATgc           <  4:312134/136‑1 (MQ=255)
                                                                                                                          cgcccgccCCCCCCGAACGCCCCCAACCT‑ACCCGGTTAAACCGCTCGAACAACCGCCACTGGTGCTCTGCAGGCACACCCCGACCCGGGTTAGGCACCGAAAGAATGACTTCGCAGCCCGCGCCTTCGAAGCCCa        >  4:164306/1‑135 (MQ=255)
                                                                                                                          cGCCCTCCAGCCGCGAACGCCCCGAACCT‑ACCCGGTAAAACCGCTCGAACAAGCGCGACTGGTGct                                                                             <  3:116745/66‑1 (MQ=255)
                                                                                                                          cGCCCTCCAGCCGCGAACGCCCCGAACCT‑ACCCGGTAAAACCGCTCGAACAAGCGCGACTGGTGct                                                                             >  4:116745/1‑66 (MQ=255)
                                                                                                                            ccgtCCCGCCGCGACCGCCCCGAACCT‑CCCCGGTAAAACGGCTCGAACAAGCGCGACTGGTGCTCTGCAGGCACACCCGGACCCTGGTTCTCCACGGATAGAATGACTTCGCTGCCCGAGCGTTCGATGCGCagag     <  4:185062/133‑1 (MQ=255)
                                                                                                                             ccTCCAGCCGCGAACGCCCCGAACCT‑ACCCGGTAAAACCGCTCGAACAAGCGCGACTGGTGCTCTGCAGGCACACCCGGACCCTGGTTCTCCACGGATAGAATGACTTCGCTGCCCGAGCGTTCGATGCGCAGAGc    >  3:139340/1‑136 (MQ=255)
                                                                                                                               tCCAGCCGCGAACGCCCCGAACCT‑ACCCGGTAAAACCGCTCGAACAAGCGCGACTGGTGCTCTGCAGGCACACCCGGACCCTGGTTCTCCACGGATAGAATGACTTCGCTGCCCGAGCGTTCGATGCGCAGAGCCa  >  4:234511/1‑136 (MQ=255)
                                                                                                                                  |                                                                                                                                     
CGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCTCCAGCCGCGAACGCCCCGAACCT‑ACCCGGTAAAACCGCTCGAACAAGCGCGACTGGTGCTCTGCAGGCACACCCGGACCCTGGTTCTCCACGGATAGAATGACTTCGCTGCCCGAGCGTTCGATGCGCAGAGCCA  >  ACB143/1630410‑1630672

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