Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1435,387,6200CG31.6% 49.9 / 19.0 38intergenic (‑1629/+22)CJ019_05070/nptIIhypothetical protein/neomycin phosphotransferase
Reads supporting (aligned to +/- strand):  ref base C (1/25);  new base G (12/0);  total (13/25)
Fisher's exact test for biased strand distribution p-value = 4.80e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TGAGCGAGGAAGCGGAAGAGCGCCCAATACGCAAACCGCCTCTGGGAGACCAGAAACAAAAAAAGGCCGCGTTAGCGGCCTTCAATAATTGGACCTGGCTCCTAACTGATTTTTAAGGCGACTGATGAGTCGCC‑TTTTTTTTGTCTAAGAATTCATCAGAAGAACTCGTCAAGAAGGCGATAGAAGGCGATGCGCTGCGAATCGGGAGCGGCGATACCGTAAAGCACGAGGAAGCGGTCAGCCCATTCGCCGCCAAGCTCTTCA  >  ACB143/5387487‑5387750
                                                                                                                                     |                                                                                                                                   
tGAGCGAGGACGCGGAAGAGCGCCCAATACGCAAACCGCCTCTGGGACACCAGAAACAAAAAAAGGCCGCGTTAGCGGCCTTCAATAATTGGACCTGGCTCCTAACTGATTTTTAAGGCGACTGATGAGTCGCC‑ttt                                                                                                                                 <  8:39393/137‑1 (MQ=255)
       ggAAGCGGAAGAGCGCCCAATACGCAAACCGCCTCTGGGAGACCAGAAACAAAAAAAGGCCGCGTTAGCGGCCTTCAATAATTGGACCTGGCTCCTAACTGATTTTTAAGGCGACTGATGAGTCGCG‑gttttttttt                                                                                                                          >  1:71986/1‑135 (MQ=255)
         aaGCGGAAGAGCGCCCAATACGCAAACCGCCTCTGGGAGACCAGAAACAAAAAAAGGCCGCGTTAGCGGCCTTCAATAATTGGACCTGGCTCCTAACTGATTTTTAAGGCGACTGATGAGTCGCC‑ttttttttttat                                                                                                                        >  3:48644/1‑133 (MQ=255)
           gCGTAAGAGCCCCAAATCCCCAACCCCCCAGGGGCACACCACAAAAAAAAAAAGGCCGCGTTAGGGGCCTTCAATAATTGGACCTGGCTCCTAACTGATTTTTAAGGCGACTGATGAGTCGCC‑TTTTTTTTGTCTaa                                                                                                                      <  3:322094/137‑1 (MQ=255)
           gCGGAAGAGCGCCCAATACGCAAACCGCCTCTGGGAGACCAGAAACAAAAAAAGGCCGCGTTAGCGGCCTTCAATAATTGGACCTGGCTCCTAACTGATTTTTAAGGCGACTGATGAGTCGCG‑TGTTTTTTTTGTaa                                                                                                                      >  8:172870/1‑137 (MQ=255)
           ccgcaaGAGCCCCAAATACCCAACCCGCCTCTGGGACACCACAAACAAAAAAAGGCCGCGTTAGCGGCCTTCAATAATTGGACCTGGCTCCTAACTGATTTTTAAGGCGACTGATGAGTCGCC‑TTTTTTTTGTCTaa                                                                                                                      <  3:35937/133‑1 (MQ=255)
               aaGAGCGCCCAATACGCAAACCGCCTCTGGGAGACCAGAAACAAAAAAAGGCCGCGTTAGCGGCCTTCAATAATTGGACCTGGCTCCTAACTGATTTTTAAGGCGACTGATGAGTCGCG‑GTTTTTTTGTGTAAGAAt                                                                                                                  >  4:205432/1‑137 (MQ=255)
               aaGAGCGCCCAATACGCAAACCGCCTCTGGGAGACCAGAAACAAAAAAAGGCCGCGTTAGCGGCCTTCAATAATTGGACCTGGCTCCTAACTGATTTTTAAGGCGACTGATGAGTCGCG‑GGTTTTTTTTGTaattat                                                                                                                  >  2:249514/1‑133 (MQ=255)
                            aCGCAAACCGCCTCTGGGAGACCAGAAACAAAAAAAGGCCGCGTTAGCGGCCTTCAATAATTGGACCTGGCTCCTAACTGATTTTTAAGGCGACTGATGAGTCGCG‑GTTTTTTTGTCTTAGGATTGAATAGAATAAc                                                                                                     >  5:492260/1‑137 (MQ=255)
                                aaaCCGCCTCTGGGAGCCCAGAAAAAAAAAAAGGCCGCTTTAGCGGCCTTCAATAATTGGACCTGGCTCCTAACTGATTTTTAAGGCGACTGATGAGTCGCC‑TTTTTTTTGTCTAAGAATTCATCAGAAGAACTCGt                                                                                                 <  8:126232/137‑1 (MQ=255)
                                    cGCCTCTGGGAGACCAGAAACAAAAAAAGGCCGCGTTAGCGGCCTTCAATAATTGGACCTGGCTCCTAACTGATTTTTAAGGCGACTGATGAGTCGCG‑TGTTTTTTTTTTTTGAAATCATCAGAAAAAATCGTCaag                                                                                             >  1:213769/1‑137 (MQ=255)
                                     gCCTCTGGGACACCAGAAAAAAAAAAAGGCCGGGTTAGCGGCCTTCAATAATTGGACCTGGCTCCTAACTGATTTTTAAGGCGACTGATGAGTCGCC‑TTTTTTTTGTCTAAGAATTCATCAGAAGAACTCGTCaaga                                                                                            <  4:48644/137‑1 (MQ=255)
                                                   aGAAACAAAAAAAGGCCGCGTTAGCGGCCTTCAATAATTGGACCTGGCTCCTAACTGATTTTTAAGGCGACTGATGAGTCGCG‑GTTTTTTTTTGTTATAATTTATGAGAAGAAAAAGGCAaaaaggcgctataaagc                                                                              >  5:337642/1‑137 (MQ=255)
                                                     aaaCAAAAAAAGGCCGCTTTAGGGGCCTTCAATAATTGGACCTGGCTCCTAACTGATTTTTAAGGCGACTGATGAGTCGCCTTTTTTTTTGTCTAAGAATTCATCAGAAGAACTCGTCAagaaggcgatagaaggcg                                                                             <  2:57981/137‑1 (MQ=255)
                                                       aaaaaaaacGGCCGCTTTAGGGGCCTTCAATAATTGGACCTGGCTCCTAACTGATTTTTAAGGCGACTGATGAGTCGCC‑TTTTTTTTGTCTAAGAATTCATCAGAAGAACTCGTCAAGAAGGCGATAGAAGGCGATg                                                                          <  7:185157/135‑1 (MQ=255)
                                                       aaaaaaaacGGCCGCGTTAGCGGCCTTCAATAATTGGACCTGGCTCCTAACTGATTTTTAAGGCGACTGATGAGTCGCC‑TTTTTTTTGTCTAAGAATTCATCAGAAGAACTCGTCAAGAAGGCGATAGAAGGCGATg                                                                          <  3:211516/135‑1 (MQ=255)
                                                              aaGGCCGCGTTAGCGGCCTTCAATAATTGGACCTGGCTCCTAACTGATTTTTAAGGCGACTGATGAGTCGCG‑GTTTTTTTGTGATATAAATAATAACAAGAAATCGGGAAGAAGGGCTTAGAAGGCGATGGgcggcg                                                                   >  8:123864/1‑137 (MQ=255)
                                                                          gCGTCCTTCAATAATTGGACCTGGCTCATAACTGATTTTTAATGCGACTGATGAGTCGCC‑TTTTTTTTGTCTAAGAATTCATCAGAAGAACTCGTCAAGAAGGCGATAGAAGGCGATGCGCTGCGAATCGGGAgcgg                                                       <  7:146105/137‑1 (MQ=255)
                                                                          gCGGCCTTCAATAATTGGACCTGGCTCCTAACTGATTTTTAAGGCGACTGATGAGTCGCC‑TTTTTTTTGTCTAGGAATTCATCAGAAGAACTCGTCAAGAAGGCGATAGAAGGCGATGCGCTGCGAATCGGGAgcgg                                                       <  2:276839/137‑1 (MQ=255)
                                                                          gCGGCCTTCAATAATTGGACCTGGCTCCTAACTGATTTTTAAGGCGACTGATGAGTCGCC‑TTTTTTTTGTCTAAGAATTCATCAGAAGAACTCGTCAAGAAGGCGATAGAAGGCGATGCGCTGCGAATCGGGAgcgg                                                       <  8:189612/137‑1 (MQ=255)
                                                                          gCGGCCTTCAATAATTGGACCTGGCTCCTAACTGATTTTTAAGGCGACTGATGAGTCGCC‑TTTTTTTTGTCTAAGAATTCATCAGAAGAACTCGTCAAGAAGGCGATAGAAGGCGATGCGCTGCGAATCGGGAgcgg                                                       <  6:488836/137‑1 (MQ=255)
                                                                          gCGGCCTTCAATAATTGGACCTGGCTCCTAACTGATTTTTAAGGCGACTGATGAGTCGCC‑TTTTTTTTGTCTAAGAATTCATCAGAAGAACTCGTCAAGAAGGCGATAGAAGGCGATGCGCTGCGAATCGGGAgcgg                                                       <  4:135724/137‑1 (MQ=255)
                                                                               cTTCAATAATTGGACCTGGCTCCTAACTGATTTTTAAGGCGACTGATGAGTCGCC‑TTTTTTTTGTCTAAGAATTCATCAGAAGAACTCGTCAAGAAGGCGATAGAAGGCGATGCGCTGCGAATCGGGAGCGGCtata                                                  <  4:455222/137‑1 (MQ=255)
                                                                                      aatTGGACCTGGCTCCTAACTGATTTTTAAGGCGACTGATGAGTCGCC‑TTTTTTTTGTCTAAGAATTCATCAGAAGAACTCGTCAAGAAGGCGATAGAAGGCGATGCGCTGCGAATCGGGAGCGGCGATACCGTaaa                                           <  8:311556/137‑1 (MQ=255)
                                                                                           gACCTGGCTCCTAACTGATTTTTAAGGCGACTGATGAGTCGCG‑GTTTTTTTTTCGTATAATTAAACAAAAGAAATACTCAAGAAGGCGCTAGAAGGCGCTGCGGGGGGAAAGGGGGGCGGCGATACCGGAAAGCAcc                                      >  4:171195/1‑136 (MQ=255)
                                                                                             cctttcTCATAACTGATTTTTAAGGCGACTGATGAGTCGGC‑TTTTTTTTGTCTAATAATTCATCAGAAGAACTCGTCAAGAAGGCGATAGAAGGCGATGCGCTGCGAATCGGGAGCGGCGATACCGTAAAGCACGAg                                    <  2:437681/132‑1 (MQ=255)
                                                                                             ccTGGCTCCTAACTGATTTTTAAGTCGACTGATGAATCTCC‑TTTTTTTTGTCTAAGAATTCATCATAAGAACTCGTCAAGAAGGCGATAGAAGGCGATGCGCTGCGAATCGGGAGCGGCGATACCGTAAAGCACGAg                                    <  4:124608/137‑1 (MQ=255)
                                                                                             ccTGGCTCCTAACTGATTTTTAAGGCGACTGATGAGTCGCC‑TTTTTTTTGTCTAAGAATTCATCAGAAGAACTCGTCAAGAAGGCGATAGAAGGCGATGCGCTGCGAATCGGGAGCGGCGATACCGTAAAGCACGAg                                    <  2:54276/137‑1 (MQ=255)
                                                                                                   tCCTAACTGATTTTTAAGGCGACTGATGAGTCGCG‑GTTTTTTTGTGTTATAATTAATAAGAAGAAATCGTCAAGAAGGCGATAGAAGGCGCTGCGCTGCGAATCGGGGGCGGCGATTCCGGAAATCACCAGGCAGCg                              >  8:177263/1‑137 (MQ=255)
                                                                                                          tGATTTTTAAGGCGACTGATGAGTCGCG‑GTTTTTTTTTGTTATAAATAATCAGAAAAAATAGGGAAGAAGGGGAAAGAATTCGATGCGCAGCGAATCGGGAGCGGGGCTACCGTTAAGCCCGAGGACGCGGTCAGcc                       >  4:82575/1‑137 (MQ=255)
                                                                                                           taTTTTTAAGGCGACTGATGAGTCGCC‑TTTTTTTTGTCTAAGAATTCATCAGAAGAACTCGTCAAGAAGGCGATAGAAGGCGAAGCGCTGCGAATCGGGAGCGGCGATACCGTAAAGCACGAGGAAGCGGTCAGccc                      <  1:447169/136‑1 (MQ=255)
                                                                                                               tttAAGGCGACTGATGAGTCGCG‑GTTTTTTTGTGTAAGAATTGATCAGAAGAACTCGGGAAGAAGGCGATAGAAGGCGCTGCGCTGCGAATCGGGGGCGGCGCTACCGTAAAGCACGAGGAAGCGGGCAGCCCAtta                  >  6:470068/1‑136 (MQ=255)
                                                                                                                   aGGCGACTGATGAGTCGCC‑TTTTTTTTGTCTAAGAATTCATCAGAAGAACTCGTCAAGAAGGCGATAGAAGGCGATGCGCTGCGAATCGGGAGCGGCGATACCGTAAAGCACGAGGAAGCGGTCAGCCCATTcgccg              <  3:200174/137‑1 (MQ=255)
                                                                                                                     gCGACTGTTGAGTCGCC‑TTTTTTTTGTCTAAGAATTCATCAGAAGAACTCGTCAAGAAGGCGATAGAAGGCGATGCGCTGCGAATCGGGAGCGGCGATACCGTAAAGCACGAGGAATCGGTCAGCCCATTCgccgcc            <  3:381358/137‑1 (MQ=255)
                                                                                                                     gCGACTGATGAGTCGCC‑TTTTTTTTGTCTAAGAATTCATCAGAAGAACTCGTCAAGAAGGCGATAGAAGGCGATGCGCTGCGAATCGGGAGCGGCGATACCGTAAAGCACGAGGAAGCGGTCAGCCCATTCgccgcc            <  7:182847/137‑1 (MQ=255)
                                                                                                                          tgatgaGTCGCC‑TTTTTTTTGTCTAAGAATTCATCAGAAGAACTCGTCAAGAAGGCGATAGAAGGCGATGCGCTGCGAATCGGGAGCGGCGATACCGTAAAGCACGAGGAAGCGGTCAGCCCATTCGCCGCCAAGct       <  1:104105/137‑1 (MQ=255)
                                                                                                                             tgaGTCGCC‑TTTTTTTTGTCTAAGAATTCATCAGAAGAACTCGTCAAGAAGGCGATAGAAGGCGATGCGCTGCGAATCGGGAGCGGCGATACCGTAAAGCACGAGGAAGCGGTCAGCCCATTCGCCGCCAAGCTCtt    <  4:8530/137‑1 (MQ=255)
                                                                                                                               aGTCGCC‑TTTTTTTTGTCTAAGAATTCATCAGAAGAACTCGTCAAGAAGGCGATAGAAGGCGATGCGCTGCGAATCGGGAGCGGCGATACCGTAAAGCACGAGGAAGCGGTCAGCCCATTCGCCGCCAAGCTCTTCa  <  8:271721/137‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                   
TGAGCGAGGAAGCGGAAGAGCGCCCAATACGCAAACCGCCTCTGGGAGACCAGAAACAAAAAAAGGCCGCGTTAGCGGCCTTCAATAATTGGACCTGGCTCCTAACTGATTTTTAAGGCGACTGATGAGTCGCC‑TTTTTTTTGTCTAAGAATTCATCAGAAGAACTCGTCAAGAAGGCGATAGAAGGCGATGCGCTGCGAATCGGGAGCGGCGATACCGTAAAGCACGAGGAAGCGGTCAGCCCATTCGCCGCCAAGCTCTTCA  >  ACB143/5387487‑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: ‑