Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0077791,253,4040AC43.2% 69.7 / 12.2 44L254F (TTA→TTC) dhaRDNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (20/5);  new base C (0/19);  total (20/24)
Fisher's exact test for biased strand distribution p-value = 3.34e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AGCAGGGTAATTTGCAATTTATTAATGCCCAGGCGGCGCGGGTCTTGCGCCTTGACGCGACGGCAAGTCAGGGACGGGCAATCACTGAACTCTTAACGTTACCCGCCGTATTGCAACAAGCAATAAAACAGGCACATCCGCTCAAACACGTAGAAGCAACCTTTGAAAGCCAGCACCAGTTTATTGATGCGGTGATAA  >  NC_007779/1253304‑1253501
                                                                                                    |                                                                                                 
aGCAGGGTAATTTGCAATTTATTAATGCCCAGGCGGCGCGGGTCTTGCGCCTTGACGCGACGGCAAGTCAGGGACGGGCAATCACTGAACTCTTAACGTTa                                                                                                   <  3:606982/101‑1 (MQ=255)
aGCAGGGTAATTTGCAATTTATTAATGCCCAGGCGGCGCGGGTCTTGCGCCTTGACGCGACGGCAAGTCAGGGACGGGCAATCACTGAACTCTTAACGTTa                                                                                                   <  3:1248721/101‑1 (MQ=255)
aGCAGGGTAATTTGCAATTTATTAATGCCCAGGCGGCGCGGGTCTTGCGCCTTGACGCGACGGCAAGTCAGGGACGGGCAATCACTGAACTCTTAACGTTa                                                                                                   >  3:1241253/1‑101 (MQ=255)
aGCAGGGTAATTTGCAATTTATTAATGCCCAGGCGGCGCGGGTCTTGCGCCTTGACGCGACGGCAAGTCAGGGACGGGCAATCACTGAACTCTTAACGTTa                                                                                                   <  3:862997/101‑1 (MQ=255)
   aGGGTAATTTGCAATTTATTAATGCCCAGGCGGCGCGGGTCTTGCGCCTTGACGCGACGGCAAGTCAGGGACGGGCAATCACTGAACTCTTAACGTTAccc                                                                                                <  3:238883/101‑1 (MQ=255)
   aGGGTAATTTGCAATTTATTAATGCCCAGGCGGCGCGGGTCTTGCGCCTTGACGCGACGGCAAGTCAGGGACGGGCAATCACTGAACTCTTAACGTTAccc                                                                                                <  3:1202163/101‑1 (MQ=255)
           ttGCAATTTATTAATGCCCAGGCGGCGCGGGTCTTGCGCCTTGACGCGACGGCAAGTCAGGGACGGGCAATCACTGAACTCTTAACGTTACCCGCCGTAtt                                                                                        >  3:1359096/1‑101 (MQ=255)
             gCAATTTATTAATGCCCAGGCGGCGCGGGTCTTGCGCCTTGACGCGACGGCAAGTCAGGGACGGGCAATCACTGAACTCTTAACGTTACCCGCCGTATTGc                                                                                      >  1:195605/1‑101 (MQ=255)
               aaTTTATTAATGCCCAGGCGGCGCGGGTCTTGCGCCTTGACGCGACGGCAAGTCAGGGACGGGCAATCACTGAACTCTTAACGTTACCCGCCGTATTGcaa                                                                                    >  4:119041/1‑101 (MQ=255)
                  ttattaATGCCCAGGCGGCGCGGGTCTTGCGCCTTGACGCGACGGCAAGTCAGGGACGGGCAATCACTGAACTCTTAACGTTACCCGCCGTATTGcaacaa                                                                                 >  1:1489152/1‑101 (MQ=255)
                      taATGCCCGGGGGGCGCGGGTCTTGCGCCTTGACGCGACGGCAAGTCAGGGACGGGCAATCACTGACCTCTTAACGTTACCCGCCGTATTGCAACAAGCaa                                                                             <  3:1500258/101‑1 (MQ=255)
                      taATGCCCGGGCGGCGGGGGTCTTGCGCCTTGACGCGACGGCAAGTCAGGGACGGGCAATCACTGAACTCTTAACGTTCCCCGCCGTATTGCAACAAGCaa                                                                             <  3:838375/101‑1 (MQ=255)
                      taATACCCAGGCGGCGCGGGTCTTGCCCCTTGACGCGCCGAAAAGTAAGGGACGGGCAATCACTGACCTCTAATCGTTCCCCGCCGTATTGCAACAAGCaa                                                                             <  1:1143446/101‑1 (MQ=255)
                       aaTGCCCAGGCGGCGCGGGTCTTGCGCCTTGACGCGACGGCAAGTCAGGGACGGGCAATCACTGAACTCTTAACGTTACCCGCCGTATTGCAACAAGCAAt                                                                            >  3:1355663/1‑101 (MQ=255)
                           cccAGGCGGCGCGGGTCTTGCGCCTTGACGCGACGGCAAGTCAGGGACGGGCAATCACTGAACTCTTAACGTTACCCGCCGTATTGCAACAAGCAATaaaa                                                                        >  3:557871/1‑101 (MQ=255)
                           cccAGGCGGCGCGGGTCTTGCGCCTTGACGCGACGGCAAGTCAGGGACGGGCAATCACTGAACTCTTAACGTTACCCGCCGTATTGCAACAAGCAATaaaa                                                                        >  3:309197/1‑101 (MQ=255)
                            ccAGGCGGCGCGGGTCTTGCGCCTTGACGCGACGGCAAGTCAGGGACGGGCAATCACTGAACTCTTAACGTTACCCGCCGTATTGCAACAAGCAATAAAAc                                                                       >  1:685803/1‑101 (MQ=255)
                              gggcggggcggTTCTTGGGCTTTGACGCGACGGGAAGGCGGGGGCGGGAAACCACTGAACTTTAAACGTTACCCGCCGTATTGCAACAAGCAATAAAACAg                                                                     <  3:753480/100‑1 (MQ=255)
                                  ggggggggTCTGGCGCCTTGACGCGACGGCAGGTCAGGGACGGGCAATCACTAAACTCTTAACGTTCCCCGCCGTATTGCAACAAGCAATAAAACAGGcac                                                                 <  3:983020/96‑1 (MQ=255)
                                  ggcgCGGGTCTTGCGCCTTGACGCGACGGAAAGTCAGGGACGGGCAATCACTGAACTCTTAACGTTCCCCGCCGTATTGCAACAAGCAATAAAACAGGcac                                                                 <  3:931124/101‑1 (MQ=255)
                                     gcgGGTCTTGCGCCTTGACGCGACGGCAAGTCAGGGACGGGCAATCCCTGAACTCTTAACGTTACCCGCCGTATTGCAACAAGCAATAAAACAGGCACATc                                                              >  1:680937/1‑101 (MQ=255)
                                       gggTCTTGCGCCTTGCCGCGCCGGCAGGTCGGGGACGGGCAACCACTAAACTTTTAACGTTCCCCGCCGTATTGCAACAAGCAATAAAACAGGCACATCCg                                                            <  3:1400759/101‑1 (MQ=255)
                                       gggTCTTGCGCCTTGACGCGACGGCAAGTCAGGGACGGGCAATCACTGAACTCTTAACGTTACCCGCCGTATTGCAACAAGCAATAAAACAGGCACATCCg                                                            >  3:1449654/1‑101 (MQ=255)
                                       gggTCTTGCGCCTTGACGCGACGGCAAGTCAGGGACGGGCAATCACTAAACTCTTAACGTTCCCCGCCGTATTGCAACAAGCAATAAAACAGGCACATCCg                                                            <  3:1256613/101‑1 (MQ=255)
                                        ggTCTTGCGCCTTGACGCGACGGCAAGTCAGGGACGGGCAATCACTGAACTCTTAACGTTACCCGCCGTATTGCAACAAGCAATAAAACAGGCACATCCGc                                                           >  1:1008445/1‑101 (MQ=255)
                                           cTTGCGCCTTGACGCGACGGCAAGTCAGGGACGGGCAATCACTGAACTCTTAACGTTACCCGCCGTATTGCAACAAGCAATAAAACAGGCACATCCGCTCa                                                        >  1:47840/1‑101 (MQ=255)
                                               cgcCTTGACGCGACGGCAAGTCAGGGACGGGCAATCACTGAACTCTTAACGTTCCCCGCCGTATTGCAACAAGCAATAAAACAGGCACATCCGCTCAAaca                                                    <  1:75188/101‑1 (MQ=255)
                                               cgcCTTGACGCGACGGAAAGTCAGGGACGGGCAATCACTGAACTCTTAACGTTCCCCGCCGTATTGCAACAAGCAATAAAACAGGCACATCCGCTCAAaca                                                    <  3:391069/101‑1 (MQ=255)
                                                      aCGCGACGGCAAGTCAGGGACGGGCAACCACTGAACTCTTAACGTTACCCGCCGTATTGCAACAAGCAATAAAACAGGCACATCCGCTCAAACACGTAGaa                                             <  3:704519/101‑1 (MQ=255)
                                                      aCGCGACGGAAAGTCGGGGACGGGCAATCACTGACCTTTTACCGTTCCCCGCCGTATTGCAACAAGCAATAAAACAGGCACATCCGCTCAAACACGTAGaa                                             <  1:1111110/101‑1 (MQ=255)
                                                       cgcgACGGCAAGTCAGGGACGGGCAATCACTGAACTCTTAACGTTACCCGCCGTATTGCAACAAGCAATAAAACAGGCACATCCGCTCAAACACGTAGAAg                                            >  1:102691/1‑101 (MQ=255)
                                                       cgcgACGGCAAGTCAGGGACGGGCAATCACTGAACTCTTAACGTTACCCGCCGTATTGCAACAAGCAATAAAACAGGCACATCCGCTCAAACACGTAGAAg                                            >  3:231079/1‑101 (MQ=255)
                                                         cgACGGAAAGCCGGGGACGGGCAATCACTGAACTCTTAACGTTCCCCGCCGTATTGCAACAAGCAATAAAACAGGCACATCCGCTCAAACACGTAGAAGCa                                          <  1:1309574/101‑1 (MQ=255)
                                                             ggaaaGTCAGGGACGGGCAATCACTAACCTCTTAACTTTCCCCGCCGTATTGCAACAAGCAATAAAACAGGCACATCCGCTCAAACACGTAGAAGCAACCt                                      <  1:265536/98‑1 (MQ=255)
                                                              gCAAGGCGGGGACGGGCAATCACGGAACTTTTAACGTTCCCCGCCGTATTGCAACAAGCAATAAAACAGGCACATCCGCTCAAACACGTAGAAGCAACCtt                                     <  1:1199235/101‑1 (MQ=255)
                                                               cAAGTCAGGGACGGGCAATCACTGAACTCTTAACGTTACCCGCCGTATTGCAACAAGCAATAAAACAGGCACATCCGCTCAAACACGTAGAAGCAACCttt                                    >  2:141660/1‑101 (MQ=255)
                                                               aaaGTCAGGGACGGGCAATCACTGAACTCTTAGCGTTCCCCGCCGTATTGCAACAAGCAATAAAACAGGCACATCCGCTCAAACACGTAGAAGCAACCttt                                    <  1:662623/100‑1 (MQ=255)
                                                                aagtacgggACGGGCCCTCCCTGAACTCTTAACGTTACCCGCCGTATTGCAACAAGCAATAAAACAGGCACATCCGCTCAAACACGTAGAAGCAACCTTTg                                   >  4:45299/7‑101 (MQ=255)
                                                                aaGTCAGGGACGGGCAATCACTGAACTCTTAACGTTACCCGCCGTATTGCAACAAGCAATAAAACAGGCACATCCGCTCAAACACGTAGAAGCAACCTTTg                                   >  4:55897/1‑101 (MQ=255)
                                                                           gggCAATCACTGAACTCTTAACGTTACCCGCCGTATTGCAACAAGCAATAAAACAGGCACATCCGCTCAAACACGTAGAAGCAACCTTTGAAAGCCAGCAc                        >  3:1421778/1‑101 (MQ=255)
                                                                                  cACTGAACTCTTAACGTTCCCCGCCGTATTGCAACAAGCAATAAAACAGGCACATCCGCTCAAACACGTAGAAGCAACCTTTGAAAGCCAGCACCAGTTTa                 <  1:650462/101‑1 (MQ=255)
                                                                                    cTGAACTCTTAACGTTCCCCGCCGTATTGCAACAAGCAATAAAACAGGCACATCCGCTCAAACACGTAGAAGCAACCTTTGAAAGCCAGCACCAGTTTAtt               <  1:1229247/101‑1 (MQ=255)
                                                                                       aaCTCTTAACGTTACCCGCCGTATTGCAACAAGCAATAAAACAGGCACATCCGCTCAAACACGTAGAAGCAACCTTTGAAAGCCAGCACCAGTTTATTGAt            >  1:503712/1‑101 (MQ=255)
                                                                                          tctaagcGTTCCCCGCCGTATTGCAACAAGCAATAAAACAGGCACATCCGCTCAAACACGTAGAAGCAACCTTTGAAAGCCAGCACCAGTTTATTGATGCg         <  3:965917/95‑1 (MQ=255)
                                                                                          tctTAACGTTCCCCGCCGTATTGCAACAAGCAATAAAACAGGCACATCCGCTCAAACACGTAGAAGCAACCTTTGAAAGCCAGCACCAGTTTATTGATGCg         <  3:748244/101‑1 (MQ=255)
                                                                                             taccGTTCCCCGCCGTATTGCAACAAGCAATAAAACAGGCACATCCGCTCAAACACGTAGAAGCAACCTTTGAAAGCCAGCACCAGTTTATTGATGCGGTg      <  3:1424930/98‑1 (MQ=255)
                                                                                             tAACGTTACCCGCCGTATTGCAACAAGCAATAAAACAGGCACATCCGCTCAAACACGTAGAAGCAACCTTTGAAAGCCAGCACCAGTTTATTGATGCGGTg      >  1:1295920/1‑101 (MQ=255)
                                                                                                 gTTCCCCGCCGTATTGCAACAAGCAATAAAACAGGCACATCCGCTCAAACACGTAGAAGCAACCTTTGAAAGCCAGCACCAGTTTATTGATTCGGGGATaa  <  4:98061/101‑1 (MQ=255)
                                                                                                    |                                                                                                 
AGCAGGGTAATTTGCAATTTATTAATGCCCAGGCGGCGCGGGTCTTGCGCCTTGACGCGACGGCAAGTCAGGGACGGGCAATCACTGAACTCTTAACGTTACCCGCCGTATTGCAACAAGCAATAAAACAGGCACATCCGCTCAAACACGTAGAAGCAACCTTTGAAAGCCAGCACCAGTTTATTGATGCGGTGATAA  >  NC_007779/1253304‑1253501

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