Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1222,094,8170TG55.6% ‑3.2 / 16.6 40Q183P (CAG→CCG) CJ019_01995hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (14/2);  new base G (11/9);  total (29/11)
Fisher's exact test for biased strand distribution p-value = 6.71e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.05e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GCCATGCGGTCCAGTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCAAGAAGCCCTCAT  >  ACB122/2094685‑2094946
                                                                                                                                    |                                                                                                                                 
gCCATGCGGTCCAGTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCggg                                                                                                                                 >  3:63924/1‑135 (MQ=255)
  cATGCGGTCCAGTACCTGGGCATTGGCCCGGGCCATCAGTTTCTCAAGATACGTCAGCTGCTATTCTTACACCTTCGGCTCCGGCTGCTTGCGGAAAAACAGACAGCCGCGCACGACGGCGCGGCCCCCCAGGcc                                                                                                                               >  5:83290/1‑135 (MQ=255)
  cATGCGGTCCAGTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGcc                                                                                                                               >  2:414573/1‑135 (MQ=255)
    tGCGGGCCAGTACCGGGGCAGGGGCCGGGGCCAACAGTTTCTCAAGAAACGTCAGCTGCTCTTCAACCACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCCCGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCgg                                                                                                                             <  8:3102/135‑1 (MQ=255)
      cGGTCCAGTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGCAGGCCACGCGCCCCGGGCCCGGc                                                                                                                           >  2:86156/1‑135 (MQ=255)
            aGTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCAGCTTGCGCAAATACTGCACGCCGCGCTCGAACGCCCGCCGCGGCAGGGCCGGCGGGccg                                                                                                                     >  8:63559/1‑134 (MQ=255)
            aGTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCGGCTCTTCATACCCCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGCAGGCCCCGCGCGCCCCGCCCGGCgggggg                                                                                                                     >  5:399690/1‑132 (MQ=255)
                      cATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGCAGGCCAGGCGCGCCGGGCCGGGCGGGCCTTTTTGCAggg                                                                                                           >  4:319107/1‑135 (MQ=255)
                      cATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGGCCCGCGCGCCGGGCCGGGCGGGGCGGTTTGGAggg                                                                                                           >  3:532671/1‑135 (MQ=255)
                      cATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCCGGCGCGCCTGGCCCGGCGGGGCGGTTTGCAggg                                                                                                           >  7:280952/1‑135 (MQ=255)
                      cATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTAGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAACGCCAGGCCCGCCAGGCCGCGCTGGCCGTTTTTGAggg                                                                                                           >  8:47156/1‑135 (MQ=255)
                              cGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCGTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCGGGGCTGGGCGTTTTGAAGGGCGtggtgg                                                                                                   >  8:207174/1‑135 (MQ=255)
                                   cAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCCGGCGCGCCTGGCCGGGCTGGCCGTTTTGCCGGGGCTGTTGGGCGAg                                                                                              >  2:417169/1‑135 (MQ=255)
                                              cAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCACCCGCGCCGTGCCGGGCCGGCCGAGTTTGAGGGCCGGGTGGGCGACGGTGTTgtaga                                                                                   >  3:47573/1‑135 (MQ=255)
                                                atatccGTCAGCGGCTCTTCAACCCCCTGCGGGCCCGGGTGCTCGCTCAATCACTGCCCCCCGCGCGCGCCGGCGCGGCGCGCGGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTgaaga                                                                                   <  8:38489/128‑1 (MQ=255)
                                                         ccctgctATTCAACCACCTCCGGCGCCGGCGGCGGCCCCAAATCCGCCCCGCCGCGCGCGACGGCCCGGCGCGCCGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATgc                                                                        <  2:402843/133‑1 (MQ=255)
                                                          gctgctCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCCCGCGCGCCGGGCCCGGGGGGCCGTTTTTGAGGGCGGGGTGGGCGCGGGTGTTGATGAACACGatattaa                                                                       >  1:119924/1‑132 (MQ=255)
                                                            tgctCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGGCGTTTTGCAGGGCCTGTTGGGCGCGGTTGTTGAAGAACTCTATATgca                                                                       >  3:189667/1‑133 (MQ=255)
                                                              ctctTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCCGGCGCGCCGGGGCGGGCGGGCCGGTTTGCAGGGCCGGGTGGGCGCGGGTGGTGGAGAACTCAATATtgagcag                                                                   >  3:10027/1‑135 (MQ=255)
                                                               tctACAAACACCTTCGGCCCCGGCTGCTGGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATgcagcag                                                                   <  8:147308/134‑1 (MQ=255)
                                                                ctTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCCCGCCTGGCCGGGCTGGCCGTTTTTCAGGGGCTGGTGGCCGAGGTGTTTGTAGAACTCTATATtcagcagcc                                                                 >  2:307344/1‑134 (MQ=255)
                                                                             cGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCGGGGCGGGCTGGCGGTTTTGCAGGGCCGGTTGGCCGAGGTTGTTTAAGAACTCTATATGCAGCAGCCCCACGATATGGCg                                                    >  5:259113/1‑135 (MQ=255)
                                                                              ggCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCCGGCGCGCCGGGCCGGGCGGGGCGTTTTGCAGGGCCTGGTGGCCGAGGGTGGTGAAGAACTCTATAATCAGCAGCCCCCGGATATGGc                                                     >  4:289667/1‑133 (MQ=255)
                                                                              ggCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCGGGCCTTTTTGCCGGGCCGGGTGGCCGAGGTTGTTGAAGAACGAGATATGCAGCAGCCCCAGGATATGGc                                                     >  8:202453/1‑133 (MQ=255)
                                                                                cTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCCGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGAt                                                 >  4:34250/1‑135 (MQ=255)
                                                                                       tgctTGCGCAAATACTGCACGCCGCGCTCGAAGGCCCGGCGCGCCGGGCCGGGCTGGGCGGTTTGGAGGGCCGGTTGGCCGAGGGTGTTGGAGAACTCTATATGGCGCAGCACCAGGCTATGGCGGGTTTGCTTg                                          >  3:205517/1‑135 (MQ=255)
                                                                                          gtGCGAAAATACTGCACGCCGCGCGCGACGGCCCGCCGCGCCGGCCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGa                                         <  2:25605/132‑1 (MQ=255)
                                                                                                           gccgcGCTCGAAGGCCCGGCGCGCCTGGCCGGGCTGGCCGTTTTGCAGGGCCGGTTGGCCGGGGTTGTTGTAGAACTCTATATGCCGCAGCCCCAGGATATGGCGGGGTTGCTTGACCCACTGCTTGGCcgcgcc                      >  3:29712/1‑133 (MQ=255)
                                                                                                           gccgcGCTCGAAGGCCAGGCGCGCCGGGCCGGGGTGGGCGTTTTTCCGGGCCTGTTGGCCGAGGTTGTTGTAGAACTCTATATTCAGCAGCACCACGATATGGGGGATTTCCTTGACCCAATGCTTgggcgcg                        >  2:314242/1‑128 (MQ=255)
                                                                                                           gccgcGCTCGAAGGCCACGCGCGCCTGGCCGGGGTGGCCGTTTTGGAGGGGCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGAc                      >  8:69021/1‑135 (MQ=255)
                                                                                                            ccgcGCGCGCCGGCCCGGCGCGCGGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACa                     <  3:319107/135‑1 (MQ=255)
                                                                                                                gcTCGAAGGCCCGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCGGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCgg                 >  4:521763/1‑135 (MQ=255)
                                                                                                                gcTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCgg                 >  4:124291/1‑135 (MQ=255)
                                                                                                                     ccggcccggcccGCCGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCaag            <  6:256621/133‑1 (MQ=255)
                                                                                                                     ccggCCGGGCGCGCGGCGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCaag            <  3:68887/133‑1 (MQ=255)
                                                                                                                     acggCCCGGCGCGCCGGCCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCaag            <  2:314919/133‑1 (MQ=255)
                                                                                                                          cggggccgcCGGGGCGGGCTGGACTTTTTGCAGGGCCTGTTGGGCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCAAGAAGcc       <  4:47573/129‑1 (MQ=255)
                                                                                                                          cAGGCGCGCCTGGCCCGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTTAAGAACTCAATATGCAGCAGCACCACGAAATGGCGGATTTGCTTGACCCAAAGCTTGGCCTCGACATCGGGCAAGAAGc        >  2:165698/1‑134 (MQ=255)
                                                                                                                            ggCGCGCCGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGGCGAGGTTGTTGAAGAACTCTATATTCAGCAGCACCAGGATATTGCGGGTTTCCTTTACACCAAGCTTGGCCTCGACACCGGGCAAGAAcccc      >  6:244375/1‑134 (MQ=255)
                                                                                                                               gcgcCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTACTTGACCCAAAGCTTGGCCTCGACAGCGGGCAAGAAGCCCTCaa  >  8:498664/1‑134 (MQ=255)
                                                                                                                               gcgcCGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCAAGAAGCCCTCAt  <  8:280952/135‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                 
GCCATGCGGTCCAGTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCAAGAAGCCCTCAT  >  ACB122/2094685‑2094946

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