Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1435,387,6200CG32.4% 48.9 / 12.7 35intergenic (‑1629/+22)CJ019_05070/nptIIhypothetical protein/neomycin phosphotransferase
Reads supporting (aligned to +/- strand):  ref base C (1/22);  new base G (10/1);  total (11/24)
Fisher's exact test for biased strand distribution p-value = 8.88e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GAGGAAGCGGAAGAGCGCCCAATACGCAAACCGCCTCTGGGAGACCAGAAACAAAAAAAGGCCGCGTTAGCGGCCTTCAATAATTGGACCTGGCTCCTAACTGATTTTTAAGGCGACTGATGAGTCGCCTTTTTTTTGTCTAAGAATTCATCAGAAGAACTCGTCAAGAAGGCGATAGAAGGCGATGCGCTGCGAATCGGGAGCGGCGATACCGTAAAGCACGAGGAAGCGGTCAGCCCATTCGCCGCCAAGCTCTTCA  >  ACB143/5387492‑5387750
                                                                                                                                |                                                                                                                                  
gAGGAAGCGGAAGAGCGCCCAATACGCAAACCGCCTCTGGGAGACCAGAAACAAAAAAAGGCCGCGTTAGCGGCCTTCAATAATTGGACCTGGCTCCTAACTGATTTTTAAGGCGACTGATGAGTCGCGgttttt                                                                                                                              >  3:120101/1‑135 (MQ=255)
  ggAAGCGGAAGAGCGCCCAATACGCAAACCGCCTCTGGGAGACCAGAAACAAAAAAAGGCCGCGTTAGCGGCCTTCAATAATTGGACCTGGCTCCTAACTGATTTTTAAGGCGACTGATGAGTCGCGgttttttt                                                                                                                            >  1:358201/1‑135 (MQ=255)
        gtaaGAGCCCCCAATACGCAACCCGCCTCTGGCACCCAAGAAACAAAAAACGGCCGCGTTAGCGGCCTTCAATAATTGGACCTGGCTCCTAACTGATTTTTAAGGCGACTGATGAGTCGCCTTTTTTTTGTCTaa                                                                                                                      <  1:78620/133‑1 (MQ=255)
        ggAAGAGCGCCCAATACGCAAACCGCCTCTGGGAGACCAGAAACAAAAAAAGGCCGCGTTAGCGGCCTTCAATAATTGGACCTGGCTCCTAACTGATTTTTAAGGCGACTGATGAGTCGCGgtttttttttg                                                                                                                         >  3:169155/1‑129 (MQ=255)
                    aaTACGCAAACCGCCTCTGGGAGACCAGAAACAAAAAAAGGCCGCGTTAGCGGCCTTCAATAATTGGACCTGGCTCCTAACTGATTTTTAAGGCGACTGATGATTCGCGTTTTTTTTTTCTAAGAATtaataa                                                                                                            >  5:268132/1‑131 (MQ=255)
                       aCGCAAACCGCCTCTGGGAGACCAGAAACAAAAAAAGGCCGCGTTAGCGGCCTTCAATAATTGGACCTGGCTCCTAACTGATTTTTAAGGCGACTGATGAGTCGCGGTTTTTTTTTTTTATAATTAAACagaaa                                                                                                        >  4:189460/1‑133 (MQ=255)
                              cccccTCTGGTACACAACAAACAAAAAACGGCCGCGTTAGGGGCCTTCAATAATTGGACCTGGCTCCTAACTGATTTTTAAGGCGACTGATGAGTCGCCTTTTTTTTGTCTAAGAATTCATCAGAAGAACTCGTc                                                                                                <  6:307062/132‑1 (MQ=255)
                              ccGCCTCTGGTAGACCAGAAAAAAAAAAAGGCCGCGTTAGCGGCCTTCAATAATTGGACCTGGCTCCTAACTGATTTTTAAGGCGACTGATGAGTCGCCTTTTTTTTGTCTAAGAATTCATCAGAAGAACTCGTc                                                                                                <  7:53041/135‑1 (MQ=255)
                                  ctctGGGAGAACAGAAACAAAAAAAGGCCGCGTTAGCGGCCTTCAATAATTGGACCTGGCTCCTAACTGATTTTTAAGGCGACTGATGAGTCGCGTTTTTTTTGTCTAAGAATTCATCAGAAGAACTCGTCaaga                                                                                            <  8:266290/135‑1 (MQ=255)
                                   tctGGGAGACCAGAAACAAAAAAAGGCCGCGTTAGCGGCCTTCAATAATTGGACCTGGCTCCTAACTGATTTTTAAGGCGACTGATGAGTCGCGTTTTTTTTGCATAAGCTTTACTCAGCGGACATCGTCaagaa                                                                                           >  2:257622/1‑135 (MQ=255)
                                         acaCCAGAAACAAAAAACGGCCGCGTTAGGGGCCTTCAATAATTGGACCTGGCTCCTAACTGATTTTTAAGGCGACTGATGAGTCGCCTTTTTTTTGTCTAAGAATTCATCAGAAGAACTCGTCAagaaggcgat                                                                                     <  1:79730/133‑1 (MQ=255)
                                               acagaaaaaaaaGGCCGCGTTAGCGGCCTTCAATAATTGGACCTGGCTCCTAACTGATTTTTAAGGCGACTGATGAGTCGCCTTTTTTTTGTCTAAGAATTCATCAGAAGAACTCGTCAagaaggcgatagaag                                                                                <  1:117297/129‑1 (MQ=255)
                                                               cggtaaGCGCCCTTCAATAATTGGACCTGGCGCCTAACTGATTTTTAAGGCGACTGATGAGTCGCCTTTTTTTTGTCTAAGAATTCATCAGAAGAACTCGTCAAGAAGGCGATAGAAGGCGATGCGCTGCGAATc                                                               <  1:135935/130‑1 (MQ=255)
                                                                 gtaaGCGGCCTTCAATAATTCGACCTGGCTCCTAACTGATTTTTAAGGCGACTGATGAGTCGCCTTTTTTTTGTCTAAGAATTCATCAGAAGAACTCGTCAAGAAGGCGATAGAAGGCGATGCGCTGCGAATCgg                                                             <  1:154410/132‑1 (MQ=255)
                                                                 gTTAGCGGCCTTCAATAATTGGACCTGGCTCCTAACTGATTTTTAAGGCGACTGATGAGTCGCGGTTTTTTTTTTTAATAATTAATCAGAAGAACTCGTCAATAAAGCGATAGAATGCGATGCGCTGCgcctgt                                                              >  2:117297/1‑129 (MQ=255)
                                                                     gCGGCCTTCAATAATTGGACCTGGCTCCTAACTGATTTTTAAGGCGACTGATGAGTCGCCTTTTTTTTGTCTAAGAATTCATCAGAAGAACTCGTCAAGAAGGCGATAGAAGGCGATGCGCTGCGAATCGGGAg                                                          <  5:124052/134‑1 (MQ=255)
                                                                            tCAATAATTGGACCTGGCTCCTTACTGATTTTTAAGGCGACTGATGAGTCGCCTTTTTTTTGTCTAAGTATTCATCAGAAGTACTAGTCAAGAAGGCGATAGAAGGCGATGCGCTGCGAATCGGGAGCGGCGATa                                                  <  5:178592/135‑1 (MQ=255)
                                                                             cAATAATTGGACCTGGCTCCTAACTGATTTTTAAGGCGACTGATGAGTCGCCGTTTTTTTTTCTAAGAATAGATAAGAAGAAAACGTCAAAAAGTCGATATAAGGCGATGCGCTGCGAATCGGGGGCGGCGATaa                                                 >  3:41619/1‑134 (MQ=255)
                                                                              aataatTGGACCTGGCTCCTAACTGATTTTTAAGGCGACTGATGAGTCGCGGTTTTTTTTTGTTAGTATTCATGAAAAGAACTAGTCAAGAATGCCATATAAGGCGATTCGCTGCGAATCGGGGGCGGCGCTAca                                                >  2:195931/1‑134 (MQ=255)
                                                                                      gACCTGGCTCCTAACTGATTTTTAAGGCGACTGATGAGTCGCGGTTTTTTTTTTTTAGAATTAAAAAGAAAAACACCTCAAAAAAGCGATAGAAGGCGATGCGCTGCGAAACGGGTGCGGCGCAAACGTAAAACa                                        >  4:190903/1‑135 (MQ=255)
                                                                                         cTGGCTCCTAACTGATTTTTAAGGCGACTGATGAGTCGCCTTTTTTTTGTCTAAGAATTCATCAGAAGAACTCGTCAAGAAGGCGATAGAAGGCGATGCGCTGCGAATCGGGAGCGGCGATACCGTAAAGCACGa                                     <  1:72981/135‑1 (MQ=255)
                                                                                          tGGCTCCTAACTGATTTTTAAGGCGACTGATGAGTCGCATTTTTTTTGTCTAAGAATTCATCAGAAGAACTCGTCAAGAAGGCGATAGAAGGCGATGCGCTGCGAATCGGGAGCGGCGATACCGTAAAGCACGAg                                    <  3:269548/135‑1 (MQ=255)
                                                                                           ggCTCCTAACTGATTTTTAAGGCGACTGATGAGTCGCGGGTTTTTTTTATAATAATTTATAATGAGAAATCGTCAACCATGCGCGAGAAGGGGATGCTCGGGGGAACGGGAGGCGCGCTGCACGAAAGCACGAgg                                   >  4:108711/1‑135 (MQ=255)
                                                                                                        tttttAAGGCGACTGATGAGTCGCCTTTTTTTTGTCAAAGAATTCATCAGAAGAACTCGTCAAGAAGGCGATAGACGGCGATGCGCTGCGAATCGGGAGCGGCGATACCGTAAAGCACGAGGAAGCGGTCAGc                        <  7:14679/133‑1 (MQ=255)
                                                                                                         ttttAAGGCGACTGATGAGTCGCCTTTTTTTTGTCTAAGAATTCATCAGAAGAACTCGTCAAGAAGGCGATAGAAGGCGATGCGCTGCGAATCGGGAGCGGCGATACCGTAAAGCACGAGGAAGCGGTCAGccc                      <  4:194162/134‑1 (MQ=255)
                                                                                                                gctaCTGATGAGTCTCCTTTTTTTTGTCTAAGAATTCATCAGAAGAACTCGTCAAGAAGGCGATAGAAGGCGATGCGCTGCGAATCGGGAGCGGCGATACCGTAAAGCACGAGGAAGCGGTCAGCCCATTcgccg              <  2:329380/132‑1 (MQ=255)
                                                                                                                gCGACTGATGAGTCGCCTTTTTTTTGTCTAAGAATTCATCAGAAGAACTCGTCAAGAAGGCGATAGAAGGCGATGCGCTGCGAATCGGGAGCGGCGATACCGTAAAGCACGAGGAAGCGGTCAGCCCATTcgccg              <  2:250798/135‑1 (MQ=255)
                                                                                                                gCGACTGATGAGTCGCCTTTTTTTTGTCTAAGAATTCATCACAAGAACTCGTCAAGAAGGCGATAGAAGGCGATGCGCTGCGAATCGGGAGCGGCGATACCGTAAAGCACGAGGAAGCTGTCAGCCCATTcgccg              <  7:340724/135‑1 (MQ=255)
                                                                                                                  gccTGATGAGTCGCCTTTTTTTTGTCTAAGACTTCATCAGAAGAACTCGTCAAGAAGGCGATAGAAGGCGATGCGCTGCGAATCGGGAGCGGCGATACCGTAAAGCACGAGGAAGCGGTCAGCCCATTCgccgcc            <  6:100110/133‑1 (MQ=255)
                                                                                                                  gACTGATGAGTCGCCTTTTTTTTGTCTAAGAATTCATCAGACGAACTCGTCAAGAAGGCGATAGAAGGCGATGCGCTGCGAATCGGGAGCGGCGATACCGTAAAGCACGAGGAAGCGGTCAGCCCATTCgccgcc            <  8:125083/135‑1 (MQ=255)
                                                                                                                  gACTGATGAGTCGCCTTTTTTTTGTCTAAGAATTCATCAGAAGAACTCGTCAAGAAGGCGATAGAAGGCGATGCGCTGCGAATCGGGAGCGGCGATACCGTAAAGCACGAGGAAGCGGTCAGCCCATTCgccgcc            <  4:146798/135‑1 (MQ=255)
                                                                                                                  gACTGATGAGTCGCCTTTTTTTTGTCTAAGAATTCATCAGAAGAACTCGGCAAGAAGGCGATAGAAGGCGATGCGCTGCGAATCGGGAGCGGCGATACCGTAAAGCACGAGGAAGCGGTCAGCCCATTCgccgcc            <  7:52326/135‑1 (MQ=255)
                                                                                                                        tgaGGCGCC‑TTTTTTTGTCTAAGAATTCATCAGAAGAACTCGTCAAGAAGGCGATAGAAGGCGATGCGCTGCGAATCGGGAGCGGCGATACCGTAAAGCACGAGGAAGCGGTCAGCCCATTCGCCGCCAAGct       <  8:207344/133‑1 (MQ=255)
                                                                                                                           gTCGCCTTTTTTTTGTCTAAGAATTCATCAGAAGAACTCGTCAAGAAGGCGATAGAAGGCGATGCGCTGCGAATCGGGAGCGGCGATACCGTAAAGCACGAGGAAGCGGTCAGCCCATTCGCCGCCAAGCTCTTc   <  5:292572/135‑1 (MQ=255)
                                                                                                                            tCGCCTTTTTTTTGTCTAAGAATTCATCAGAAGAACTCGTCAAGAAGGCGATAGAAGGCGATGCGCTGCGAATCGGGAGCGGCGATACCGTAAAGCACGAGGAAGCGGTCAGCCCATTCGCCGCCAAGCTCTTCa  <  3:14397/135‑1 (MQ=255)
                                                                                                                                |                                                                                                                                  
GAGGAAGCGGAAGAGCGCCCAATACGCAAACCGCCTCTGGGAGACCAGAAACAAAAAAAGGCCGCGTTAGCGGCCTTCAATAATTGGACCTGGCTCCTAACTGATTTTTAAGGCGACTGATGAGTCGCCTTTTTTTTGTCTAAGAATTCATCAGAAGAACTCGTCAAGAAGGCGATAGAAGGCGATGCGCTGCGAATCGGGAGCGGCGATACCGTAAAGCACGAGGAAGCGGTCAGCCCATTCGCCGCCAAGCTCTTCA  >  ACB143/5387492‑5387750

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