Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1434,519,9650AT50.0% ‑0.5 / 26.1 28intergenic (+43/‑39)purB/roxAAdenylosuccinate lyase/50S ribosomal protein L16 3‑hydroxylase
Reads supporting (aligned to +/- strand):  ref base A (4/10);  new base T (14/0);  total (18/10)
Fisher's exact test for biased strand distribution p-value = 1.53e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

AGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGACTACTGGCAGAAGAAGCCGCTGCTGGTACGCC  >  ACB143/4519841‑4520097
                                                                                                                              |                                                                                                                                     
aGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGaa                                                                                                                           <  8:124328/137‑1 (MQ=255)
     aaGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGGTTTTTTTTGCCGTAGAAAAttt                                                                                                                      >  3:155104/1‑136 (MQ=255)
     aaGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAAAGCCCCGCCTGGCCGGGCGGTTTTTTTTGCCGTAGAAAAttt                                                                                                                      >  2:84627/1‑136 (MQ=255)
            ggCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGttttt                                                                                                               <  3:354386/137‑1 (MQ=255)
              ccGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTc                                                                                                             <  6:364354/137‑1 (MQ=255)
               cGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTTTTTCCGTAAGGAAATTAAGttttttga                                                                                                            >  6:112183/1‑135 (MQ=255)
                   cTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGAGTACAACGCCCGGCCAGGCCGGGCGTGTTTTTTTCCGTGAGAAAATTACGTTTTTTAAATat                                                                                                        >  6:183737/1‑136 (MQ=255)
                    tCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGGCGGGCGTTTTTATTTCCCGACGAAAATTAAGTTTTTTTaaaaag                                                                                                       >  7:369768/1‑133 (MQ=255)
                    tCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTTTTACCCGACGAAAAATAAGGTTTTTGAATatg                                                                                                       >  2:39760/1‑135 (MQ=255)
                    tCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAgg                                                                                                       <  6:204073/137‑1 (MQ=255)
                    tCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCGGGCCGGGCGGTGTTTTTTCCGGAGCGAAATTAAGGTTTTTTAATTgg                                                                                                       >  1:258221/1‑137 (MQ=255)
                    tCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTAAAAAGCCCGGCCGGGGCGGGCGGTTTTTTTTCCGGAACAAAATTTAATTTTTTGAATatg                                                                                                       >  7:338550/1‑135 (MQ=255)
                    tCAAGCAGCTGACCCCCGCTACCAACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACCCCCGGGCTGGCCGGGCGTGTTTTTTAAGGGACGAAAATTAGGTTTTTTCaaacag                                                                                                       >  6:68543/1‑133 (MQ=255)
                               aCCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCGGGCCGGGCGGTTTTTTTTCCGGACGCAAAATACGGTTTTTTAATATGTTGGACATGaa                                                                                            >  7:243084/1‑137 (MQ=255)
                                 ccccGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTACCGGACGAAAATTAAGGTTTTTAAATATGTTTAAAAAGAAta                                                                                          >  7:107524/1‑136 (MQ=255)
                                  cccGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATcc                                                                                         <  3:223479/137‑1 (MQ=255)
                                    cGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTg                                                                                       <  3:108528/137‑1 (MQ=255)
                                         cctacATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTGTTTTTTTCCGGACGCAAATTTAGGTTTTTTAAAATGTTGTACAAGAATACGGATaaa                                                                                  >  4:223479/1‑135 (MQ=255)
                                                cGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTACCGTTCGAAAATTAAGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGc                                                                           >  2:202753/1‑137 (MQ=255)
                                                 ggTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTTTTGCCGGTCGAAAATTAAGGTTTTTGAATTTGTTGAACAAGAATACGGATACACCACTGcc                                                                          >  4:174192/1‑136 (MQ=255)
                                                               ccAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCa                                                            <  8:303247/137‑1 (MQ=255)
                                                                      aaaCGCATCTAAGCTGAGCGTCGCGTACAACGCCCCGCCTGGCCGGGCGGTTTTTTTTCCGGACGAAAATTTAGTTTTTTGAAAAAGGTTTACATGAATCCGTTAAAACCACAGCCCCTGCTCGGCGGCAtcgaggc                                                     >  8:226636/1‑132 (MQ=255)
                                                                                  gCTGAGCGTCGAGTACAACGCCCGGCCTGGCCGGGAGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCTCGAATTCAt                                        <  3:312102/137‑1 (MQ=255)
                                                                                         gTCGCGTACAACGCCCGGCCTGGCCGGGCGTGTTTTTTAGCGGAACCAAAATAAGGTTTTTTAATTGGTTGGAAAAGAATACGGCTAAGACCCATCAGCAGCGCGGCGGCATCATCGGCCCCGGATTGATTcgcggc                                 >  8:106476/1‑135 (MQ=255)
                                                                                                          gCCTGGCCGGGCGTTTTTTTTCCCGTAAGAAAATTTAGGTTTTTCAAAAGGTTGGAAATTAATCCTGCTACTCCACAGAAGCTGCTCGGCGGGATCACGGCCCGCCAATTCAAGCGGGCAATCTGGCAGAATAAGcc                >  6:162758/1‑137 (MQ=255)
                                                                                                          gCCTGGCCGGGCGTTTTTATTGCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGTATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGACTACTGGCAGAAGAAGcc                >  6:335428/1‑137 (MQ=255)
                                                                                                               gCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGACTACTGGCAGAAGAAGCCgctgc           <  4:155104/137‑1 (MQ=255)
                                                                                                                          ttttATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGACTACTGGCAGAAGAAGCCGCTGCTGGTACGcc  <  1:184467/137‑1 (MQ=255)
                                                                                                                              |                                                                                                                                     
AGCCGAAGCCAAGGCCGAGCTCAAGCAGCTGACCCCCGCTACCTACATCGGTAACGCGGCAGCCCAGGCCAAACGCATCTAAGCTGAGCGTCGCGTACAACGCCCGGCCTGGCCGGGCGTTTTTATTCCCGTACGAAAATTACGTTTTTTCAATAGGTTGAACATGAATCCTGATACTCCACTGCAGCTGCTCGGCGGCATCTCGGCCCGCGAATTCATGCGCGACTACTGGCAGAAGAAGCCGCTGCTGGTACGCC  >  ACB143/4519841‑4520097

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: