Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1223,995,4180TG52.9% 3.8 / 11.4 17intergenic (+90/‑44)dmoA_3/allEDimethyl‑sulfide monooxygenase/(S)‑ureidoglycine aminohydrolase
Reads supporting (aligned to +/- strand):  ref base T (0/8);  new base G (0/9);  total (0/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.66e-02
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GGCACGTTGGACGTCAGGTCAGCATTCAGTGATCAAATACCGGGTATTCCAACATTTAATTTGTGGATGCCACGCACTCGACCCTATTCTGTGGCCACATCCCCTGCGCCGCACCGGCCGGTGCGGTTCAAGGCATGGGGAGAGAACAACAAGTCGAGATCGTCCATGTCGAACCATTCATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCT  >  ACB122/3995297‑3995549
                                                                                                                          |                                                                                                                                   
ggCACGTTGGACGTCAGGTCAGCATTCAGTGATCAAAAACCGGGTATTCCAACATTTAATTTGTGGATGCCACGCACTCGACCCTATTCTGTGGCCACATCCCCTGCGCCGCACCGGCCGGTGCGGTTCAAGGCAt                                                                                                                       <  8:114667/136‑1 (MQ=255)
          aggtcaggtacgCATTATGTGATCAATACCCGGTATTTCCAACATTTAATTTGGGTATGCCCCGCACTCCACCATATTCTGGGGCCACAGCCCCTGCGCCGCACCGGCGGGGGCGGTTCAAGGCATGGGGTGAGaa                                                                                                             <  6:260100/134‑1 (MQ=255)
          aCGTCAGGTCAGCATTCAGTGATCAAATACCGGGTATTCCAACATTTAATTTGTGGATGCCACGCACTCGACCCTATTCTGTGGCCACATCCCCTGCGCCGCACCGGCCGGTGCGGTTCAAGGCATGGGGAGAGaa                                                                                                             <  5:130278/136‑1 (MQ=255)
           cGTCAGGTCAGCATTCAGTGATCAAATACCGGGTATTCAAACATTTAATTTGTGGATGCCACGCACTCGACCCTATTCTGTGGCCACATCCCCGGCGCCGCACCGGCCGGTGCGGTTCAAGGCATGGGGAGAGaac                                                                                                            <  4:50879/136‑1 (MQ=255)
                                 cAAATACCGGGTATTCCAACATTTAATTTGTGGATGCCACGCACTCGACCCTATTCTGTGGCCACATCCCCTGCGCCGCACCGGCCGGGGCGGTTCAAGGCATGGGGAGAGAACAACAAGTCGAGATCGTCCATGt                                                                                      <  3:39733/136‑1 (MQ=255)
                                 cAAATACCGGGTATTCCAAAATTTAAGTTGTGGATGCCACGCACTCGACCCTATTCTGGGGCCACATCCCCGGCGCCGCACCGGCGGGTGCGGTTCAAGGCATGGGGAGAGAACAACAAGTCGAGATCGTCCATGt                                                                                      <  5:13767/136‑1 (MQ=255)
                                 aaaaacccGGGAATACAAAAATTTAATTGGTGGATCCCACGCACGCGCCCCTATTCTGGGCCCACACCCCCTGCGCCGCCCCGGCCGGGGCGGTTCAAGGCATGGGGAGAGAACAACAAGTCGAGATCGTCCATGt                                                                                      <  6:333002/130‑1 (MQ=255)
                                  aaaTCCGGGGTATTCCAACATTTATTTTGTGGATGCCACGCACTCGACCCTATTCTGTGGCCACATCCCCTGCGCCGCACCGGCCGGTGCGGTTCAAGGCATGGGGAGAGAACAACAAGTCGAGATCGTCCATGt                                                                                      <  4:88374/135‑1 (MQ=255)
                                                                                gACCGTTTTCTGTGCCCACTGCCCCGGCGCCGCCCCGCCGGGGGCGGTTGAAGGCATGGGGAGAGAACAACAAGTCGAGATCGTCCATGTCGAACCATTCATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAAc                                        <  4:286776/136‑1 (MQ=255)
                                                                                  cccTATTCTGTGGCCACAGCCCCGGCGCCGCACGGGCCGGGGCGGTTCAAGGCATGGGGAGAGAACAACAAGTCGAGATCGTCCATGTCGAACCATTCATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCg                                      <  7:76566/136‑1 (MQ=255)
                                                                                  cccTATGGTGTGGCCACATCCCCTGCGCCGCACGGGCGGGTGCGGTTCAAGGCATGGGGAGAGAACAACAAGTCGAGATCGTCCATGTCGAACCATTCATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCg                                      <  4:269234/136‑1 (MQ=255)
                                                                                             ggCCACATCCCCTGCGCCGCACCGGCCGGTGCGGTTCAAGGCATGGGGAGAGAACAACAAGTCGAGATCGTCCATGTCGAACCATTCATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGctgac                            <  2:232394/135‑1 (MQ=255)
                                                                                                 ccagccccGGCGCCCCACCGGCCGGGGCGGTTCAAGGCATGGGGAGAGAACAACAAGTCGAGATCGTCCATGTCGAACCATTCATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCtgactgac                        <  2:136892/131‑1 (MQ=255)
                                                                                                    tCCCCTGCGCCGCACCGGCCGGGGCGGTTCAAGGCATGGGGAGAGAACAACAAGTCGAGATCGTCCATGTCGAACCATTCATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGt                     <  2:50713/135‑1 (MQ=255)
                                                                                                              ccccccggcgggGGCGGTTCAAGGCATGGGGAGAGAACAACAAGTCGAGATCGTCCATGTCGAACCATTCATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCa           <  1:328930/131‑1 (MQ=255)
                                                                                                                      cgggTGCGGTTCAAGGCATGGGGAGAGAACAACAAGTCGAGATCGTCCATGTCGAACCATTCATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCt  <  6:104880/134‑1 (MQ=255)
                                                                                                                      ccggGGCGGTTCAAGGCATGGGGAGAGAACAACAAGTCGAGATCGTCCATGTCGAACCATTCATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCt  <  8:37577/136‑1 (MQ=255)
                                                                                                                          |                                                                                                                                   
GGCACGTTGGACGTCAGGTCAGCATTCAGTGATCAAATACCGGGTATTCCAACATTTAATTTGTGGATGCCACGCACTCGACCCTATTCTGTGGCCACATCCCCTGCGCCGCACCGGCCGGTGCGGTTCAAGGCATGGGGAGAGAACAACAAGTCGAGATCGTCCATGTCGAACCATTCATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCT  >  ACB122/3995297‑3995549

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: