Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,031,2480TA51.6% 25.7 / 38.4 34T312T (ACT→ACA) pcaYaromatic acid chemoreceptor
Reads supporting (aligned to +/- strand):  ref base T (4/11);  new base A (16/0);  total (23/11)
Fisher's exact test for biased strand distribution p-value = 1.61e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.42e-04
Rejected as consensus: Frequency below/above cutoff threshold.

GACACTCTGGGGCAGATCAACAACTGTGCCCGGCAGCTGGACGCCGCCGCCACGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAACCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCACTGCCG‑TGACCGAAATGACCACGGCAGTGGAGGAAGTTGCGCGCAACGCGATCACCACCTCGCAAACCACCAGTGAATCCAACCAGCTGGCCGCGCAAAGCCGCCGGCAAGTCAGCGAAAACATCGACGGC  >  NC_002947/3031120‑3031377
                                                                                                                                |                                                                                                                                  
gACACTCTGGGGCAGATCAACAACTGTGCCCGGCAGCTGGACGCCGCCGCCACGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAACCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCACTGCCG‑TGa                                                                                                                            <  1:481526/136‑1 (MQ=255)
    ctctGGGGCAGATCAACAACTGTGCCCGGCAGCTGGACGCCGCCGCCACGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAACCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCACAGCCG‑GGGCCgc                                                                                                                        >  2:1069339/1‑135 (MQ=255)
          tgCAGATCAACAACTGTGCCCGGCAGCTGGACGCCGCCGCCACGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAACCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCCCGGCCC‑GGGCCGAAAttgg                                                                                                                  >  1:85912/2‑133 (MQ=255)
          ggCAGATCAACAACTGTGCCCGGCAGCTGGACGCCGCCGCCACGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAACCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCACAGCCG‑TGACCCaaaagcc                                                                                                                  >  2:327807/1‑132 (MQ=255)
            cAGATCAACATCTGTGCCGGGCAGCGGGACGCCGCCGCCACGGCGTTGAACGCCGTCACCGAGTAAAGCGCCAACAACCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCACTGCCG‑TGACCGAAATGACCa                                                                                                                <  1:550352/136‑1 (MQ=255)
            cAGATCAACAACTGTGCCCGGCATCTGGACGCCGCCGCCACGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAACCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCACTGCCG‑TGACCGAAATGACCa                                                                                                                <  1:509537/136‑1 (MQ=255)
              gATCAACAACTGTGCCCGGCAGCTGGACGCCGCCGCCACGGCGTTGAACGCCGTCACCGAGGAACGCGCCAACAACCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCACTGCCG‑TGACCGAAATGACCACg                                                                                                              <  2:217279/136‑1 (MQ=255)
                  aacaacTGTGCCCGGCAGCTGGACGCCGCCGCCACGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAACCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGGCACTTCCG‑TGACCCAAAAGACCACgggcg                                                                                                          >  2:313685/1‑133 (MQ=255)
                  aacaacTGTGCCCGGCAGCTGGACGCCGCCGCCACGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAACCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCCCAGACG‑GGGCCGCAAAGACCCCCGCAg                                                                                                          >  1:19281/1‑136 (MQ=255)
                  aacaacTGTGCCCGGCAGCTGGACGCCGCCGCCACGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAACCTGCGTCAGCACGGGCACGAGCTGGAGCAGGCCGCCCCAGGCGATGACC‑AAAAGACCACCGCac                                                                                                          >  2:975561/1‑135 (MQ=255)
                            cccGGCAGCTGGACGCCGCCGCCCCGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAACCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCCCGGCCG‑GGGCCCAAAAGACCCCCGCCGTGGGGGGAGt                                                                                                >  2:1002473/1‑136 (MQ=255)
                            cccGGCAGCTGGACGCCGCCGCCACGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAACCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCACAGCCG‑TGGCCGAAATGACCCCGGCAGTGGAGGAAGt                                                                                                >  2:305415/1‑136 (MQ=255)
                            cccGGCAGCTGGAAGCCGCCGCCACGGCGTTGAACGCCGTAACCGAGGAAAGCGCAAAAAACCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCCCAGGCG‑TTACCGCAATGACCACCGCAGTGGCGGAAgg                                                                                                >  2:17107/1‑135 (MQ=255)
                             ccGGCAGCTGGACGCCGCCGCCACGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAACCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCCCAGCCC‑TGTCCCAAATGACCACCGGAGTGGGGGGAGtt                                                                                               >  1:311946/1‑136 (MQ=255)
                             ccGGCAGCTGGACGCCGCCGCCACGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAACCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCACTGCCG‑TGACCGAAATGACCACGGCAGTGGAGGAAGtt                                                                                               >  2:360835/1‑136 (MQ=255)
                                  aGCTGGACGCCGCCGCCACGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAACCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCCCAGCCG‑GGACCGAAAAGACCACCGCCGTGGAGGAAGTTgcgcg                                                                                          >  2:861185/1‑136 (MQ=255)
                                  aGCTGGACGCCGCCGCCACGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAACCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCACAGCCG‑GGGCCGCAAAGACCCCCGCAGTGGAGGAAGTTtcgcg                                                                                          >  2:42781/1‑136 (MQ=255)
                                  aGCTGGACGCCGCCGCCACGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAACCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCACAGCCG‑GGCCCGCAATGACCCCGGCAGGGGGGGCAGTTtcgcg                                                                                          >  2:780756/1‑136 (MQ=255)
                                    cTGGACGCCGCCGCCACGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAACCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCACTGCCG‑TGACCGAAATGACCACGGCAGTGGAGGAAGTTGCGCGCa                                                                                        <  1:1036920/136‑1 (MQ=255)
                                                 ccACGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAACCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCACAGCCG‑GGACCGAAAAGACCACGGCAGGGGAGGAAGTTGCGCGCAACGCGAAaaccac                                                                           >  2:723607/1‑136 (MQ=255)
                                                    cGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAACCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCACTGCCG‑TGACCGAAATGACCACGGCAGTGGAGGAAGTTGCGCGCAACGCGATCACCACCTc                                                                        <  1:1263885/136‑1 (MQ=255)
                                                       ggTTGAACGCGGTCACCGAGGAAAGCGCCAACAACCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCACTGCCG‑TGACCGAAATGACCACGGCAGTGGAGGAAGTTGCGCGCAACGCGATCACCACCTCGCa                                                                     <  1:1097889/135‑1 (MQ=255)
                                                                ccGTCACCGAGGAAAGCGCCAACAACCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCACTGCCG‑TGACCGAAATGACCACGGCAGTGGAGGAAGTTGCGCGCAACGCGATCACCACCTCGCAAACCACCAg                                                            <  1:1037813/136‑1 (MQ=255)
                                                                    cACCGAGGAAAGCGCCAACAACCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCACTGCCG‑TGACCGAAATGACCACGGCAGTGGAGGAAGTTGCGCGCAACGCGATCACCACCTCGCAAACCACCAGTGaa                                                        >  2:514134/1‑136 (MQ=255)
                                                                    cACCGAGGAAAGCGCCAACAACCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCACGACCC‑GGGCCGCCATGACCCCGGCCGTGGGGGTAGTTTGGGGCACCGCGCTCCCCCCCGCGGCAAACACCCGGGaa                                                        >  2:1164724/1‑136 (MQ=255)
                                                                    cACCGAGGAAAGCGCCAACAACCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCACATCCG‑TGACCGAAATTACCACGGCAGTGGAGGAAGTTGCGCGCAACGCGATCACCACCAAGCAAACCCCCAGTGaa                                                        >  2:863214/1‑136 (MQ=255)
                                                                    cACCGAGGAAAGCGCCAACAACCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCACAGCCG‑TGACCGAAAATCCCACGGCACTGGAGGAAGTTTCGCGCCACGCCCTCCCCCCCACGCAAACCACCAGTGaa                                                        >  1:605090/1‑136 (MQ=255)
                                                                                cgcCAACAACCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCACAGCCG‑TGACCGAAATGACCACGGCAGTGGAGGAAGTTGCGCGCAACGCGCTCACCACCTCGCAAACCACCAGTGAATCCAACCAACTg                                            >  2:348607/1‑136 (MQ=255)
                                                                                     acaacCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCACTGCCG‑TGACCGAAATGACCACGGCAGTGGAGGAAGTTGCGCGCAACGCGATCACCACCTCGCAAACCACCAGTGAATCCAACCAGCTGGCcgc                                       <  1:1324801/136‑1 (MQ=255)
                                                                                          cTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCACTGCCG‑TGACCGAAATGACCACGGCAGTGGAGGAAGTTGCGCGCAACGCGATCACCACCTCGCAAACCACCAGTGAATCCAACCAGCTGGCCGCGCaaa                                  <  2:1333837/136‑1 (MQ=255)
                                                                                                                tGGAGCAGGCCGCCACAGCCG‑GGGCCGAAATGACCACGGCAGTGGAGGAAGTTTCGCGCAACGCGATCACCCCCTCGCAAAACACCAGTGAATACCACCACCTGGCCGCGCAAAGCCGCCCGCAAGTCAGCGaaaa            >  2:1040898/1‑136 (MQ=255)
                                                                                                                     cAGGCCGCCCCAGCCG‑GGGCCGAAAATTCCACCGCAGTGGTGGAAGTTGGGCCCCACGCCCGGACCACCGCGCCAACCCCCAGTGAATACCACCCCCTGGCCGCGCAAAGCCGCCGGCCAGTGAGCGGAAACAtag       >  2:18809/1‑134 (MQ=255)
                                                                                                                      aGGCCGCCACTGCCG‑TGACCGAAATGACCACGGCAGTGGAGGAAGTTGCGCGCAACGCGATCACCACCTCGCAAACCACCAGTGAATCCAACCAGCTGGCCGCGCAAAGCCGCCGGCAAGTCAGCGAAAACATCGa      <  1:514134/136‑1 (MQ=255)
                                                                                                                          cgccACTGCCG‑GGACCGAAATGACCACGGCAGTGGAGGAAGTTGCGCGCAACGCGCTCACCACCTCGCAAACCACCAGGGAATCCAACCAGCTGGCCCCGCCACGCCGCCGGCCAGTCAGCGAAAACATCGACggg  >  2:722452/1‑135 (MQ=255)
                                                                                                                                |                                                                                                                                  
GACACTCTGGGGCAGATCAACAACTGTGCCCGGCAGCTGGACGCCGCCGCCACGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAACCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCACTGCCG‑TGACCGAAATGACCACGGCAGTGGAGGAAGTTGCGCGCAACGCGATCACCACCTCGCAAACCACCAGTGAATCCAACCAGCTGGCCGCGCAAAGCCGCCGGCAAGTCAGCGAAAACATCGACGGC  >  NC_002947/3031120‑3031377

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