Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1221,740,4070TG46.5% 4.6 / 15.6 43E53A (GAG→GCG) CJ019_01643hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/23);  new base G (1/19);  total (1/42)
Fisher's exact test for biased strand distribution p-value = 4.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TGTCCACCTTCTCTTGGAAGGCCGGGTGCATTGTTGTCTGCATGGATGATCCTCGGGTGGGTCAGGCGTGGTATTCGAAGGCCTCGGCCTTACGAACGATTCGAACTTGGGCGGTGCGCCGCTCCGGCGCGCGGCGGTCGCGGCGCATGGGGTCGCTGTCGTTGATGGCGGCATGCATTGCCACAAGACAGACCAGGGCGAGGCAGAGCGGAGTGATGATCTTCTGGTATGCCGCCTCTGCTACGAGAGCAGT  >  ACB122/1740285‑1740537
                                                                                                                           |                                                                                                                                  
tgtCCACCTTCTCTTGGAAGGCCGGGTGCATTGTTGTCTGCATGGATGATCCTCGGGTGGGTCAGGCGTGGTATTCGAAGGCCTCGGCCTTACGAACGATTCGAACTTGGGCGGTGCGCCGCGCCGGCGCGcggcgg                                                                                                                      <  3:321648/137‑1 (MQ=255)
tgtCCACCTTCTCTTGGAAGGCCGGGTGCATTGTTGTCTGCATGGATGATCCTCGGGTGGGTCAGGCGGGGTATTCGACGGCCTCGGCCTTACGAACGATTCCAACTTGGGCGGGGCGCCGCTCCGGCGCGcggcgg                                                                                                                      <  2:251014/137‑1 (MQ=255)
ggtCCACCTTCTCTTGGATGGCCGGGTGCTTTGTTGTCTGAATGTATGATCCTCGGGTGGGTCAGGCGTGGTATTCGAAGGCCTCGGCCTTACGAACGATTCGAACTTGGGCGGTGCGCCGCTCGGGCGCGcggcgg                                                                                                                      <  8:129684/136‑1 (MQ=21)
  cccACCTTCTCTTGCAAGGCCGGGTGCTTGGGTGGCGGCATGGCTGATCCTCGGGTGGGGCAGGCGTGGTATTCGAAGGCCGCGGCCTTACGAACGATTCGAACTTGGGCGGGGCGCCGCTCCGGCGCGCGGCGGTc                                                                                                                    <  2:108765/136‑1 (MQ=18)
                         gTGCATTGTTGTCTGCATGGATGATCCTCGGGTGGGTCAGGCGTGGTATTCGAAGGCCTCGGCCTTACGAACGATTCGAACTTGGGCGGTGGGCCGCGCCGGCCCCCCGCGGGCGCGCCGCCAGGGGgcgcaggggg                                                                                             >  8:198876/1‑131 (MQ=17)
                           gCATTGTTGGCTGCATGGATGATCCTCGGGGGGGGCAGGCGTGGTATTCGAAGGCCTCGGCCTTACTAACGATTCGAACTTGGGCGGTGCGCCGCTCCGGCGCGCGGCGGTCGCGGCGCATGGGGTCGCTGTCGTTg                                                                                           <  6:120966/137‑1 (MQ=255)
                                      tccATGGATGATCCTCGGGTGGGTCAGGCGTGGTATTCGAAGGCCTCGGCCTTACGAACGATTCGAACTTGGGCGGTGCGCCGCGCCGGCGCGCGGCGGTCGCGGCGCATGGGGTCGCTGTCGTTGATGGCGgcatg                                                                                <  4:411378/135‑1 (MQ=255)
                                      tGCTTGGATGATCGTCGGGGGGGGCAGGCGTGGTATTCGAAGGCCTCGGCCTTACGAACGATTCGAACTTGGGCGGTGCGCCGCTCCGGCGCGCGGCGGTCGCGGCGCATGGGGTCGCTGTCGTTGATGGCGgcatg                                                                                <  6:431294/137‑1 (MQ=255)
                                      tGCATGGATGATCCTGGGGTGGGTCAGCGGGTGTATGCGAAGGCCGCGGCCTTACGATCGATTCTAACTTGGGCGGTGCGCCGCGCCGGCGCGCGGAGGTCGCGGCGCATGGGGTCGCTGTCGTTGATGGCGgcatg                                                                                <  8:34019/137‑1 (MQ=18)
                                      tGCATGGATGATCCTCGGGTGGGGCAGGCGGGGTATTCGAAGGCCTCGGCCTTACGAACGATTCGAACTTGGGCGGTGCGCCGCTCCGGCGCGCGGCGGTCGCGGCGCATGGGGTCGCTGTCGTTGATGGCGgcatg                                                                                <  8:99321/137‑1 (MQ=255)
                                      ggcttGGATGACCCTGGGGTGGGTCTGGCGGTGTATTCGACGGCCTGGGCGTTACGAACGATTCGAAGTTGGGCGGTGCGCCGCGCCGGCGCGCGGCGGTCGCGGCGCATGGGGTCGCTGTCGTTGATGGCGgcatg                                                                                <  2:405521/133‑1 (MQ=17)
                                           ggattctCCTCGGGTGGGTCAGGCGTGGTATTCGAAGGCCTCGGCCTTACGAACGATTCGAACTTGGGCGGTGCGCCGCTCCGGCGCGCGGCGGTCGCGGCGCATGGGGTCGCTGTCGTTGATGGCGGCATGCATTg                                                                           <  1:250924/131‑1 (MQ=255)
                                           ggaggaTCCTGGGGGGGGTCAGGCGTGGTATTCGAAGGCCTCGGCCTTACGAACGATTCTAACTTGGGCGGTGCGCCGCTCCGGCGCGCGGCGGTCGCGGCGCATGGGGTCGCTGTCGTTGATGGCGGCATGCATTg                                                                           <  1:32192/137‑1 (MQ=255)
                                           ggATGATCCTCGGTGGGGTCAGGCGTGGTATTCGAAGGCCTCGGCCTTACGAACGATTCGAACTTGGGCGGTGCGCCGCTCCGGCGCGCGGCGGTCGCGGCGCATGGGGTCGCTGTCGTTGATGGCGGCATGCATTg                                                                           <  8:18768/137‑1 (MQ=255)
                                           ggATGATCCTCGGGTGGGTCAGGCGTGGTATTCGAAGGCCTCGGCCTTACGAACGATTCGAACTTGGGCGGTGCGCCGCTCCGGCGCGCGGCGGTCGCGGCGCATGGGGTCGCTGTCGTTGATGGCGGCATGCATTg                                                                           <  6:87340/137‑1 (MQ=255)
                                           ggATGATCCGCGGGTGGGTCAGGCGTGGTATTCGAAGGCCTCGGCCTTACGAACGATTCGAACTTGGGCGGTGCGCCGCGCCGGCGCGCGGCGGTCGCGGCGCATGGGGTCGCTGTCGTTGATGGCGGCATGCATTg                                                                           <  7:8363/137‑1 (MQ=255)
                                                ctCCTGGGGTGGGTCAGGCGTGGTATTCGACGGCCTCGGCCTTACGAACGATTCGAACTTGGGCGGTGCGCCGCGCCGGCGCGCGGCGGTCGCGGCGCATGGGGTCGCTGTCGTTGATGGCGGCATGCATTGCcaca                                                                      <  4:221002/136‑1 (MQ=255)
                                                cgcgtCGGGGGGGTCCGCCGTGGTATACGAAGGCCACGGCCTTACGAACGATTCGACCTTGGGCGGGGCGCCGCGCCGGCGCGCGGCGGTCGCGGCGCATGGGGTCGCTGTCGTTGATGGCGGCATGCATTGCcaca                                                                      <  6:370425/133‑1 (MQ=17)
                                                     ggggtgggtCAGGCGTGTTATTCGAAGGCCTCGGCCTTACGAACTATTCGAACTTGGGCGGGGCGCCGCTCCGGCGCGCGGCGGTCGCGGCGCATGGGGTCGCTGTCGTTGATGGCGGCATGCATTGCCACAagaca                                                                 <  1:277326/136‑1 (MQ=255)
                                                     ggggtgggtCAGGCGTGGTATTCGAAGGCCTCGGCCTTACGAACGATTCGAACTTGGGCGGTGCGCCGCTCCGGCGCGCGGCGGTCGCGGCGCATGGGGTCGCTGTCGTTGATGGCGGCATGCATTGCCACAagaca                                                                 <  7:99285/136‑1 (MQ=255)
                                                     ggggggggtCTGGCGTGGTATTCGAAGGCCTCGGCCTTACGAACGATTCGAACTTGGGCGGGGCGCCGCGCCGGCGCGCGGCGGTCGCGGCGCATGGGGTCGCTGTCGTTGATGGCGGCATGCATTGCCACAagaca                                                                 <  7:53580/136‑1 (MQ=255)
                                                          gggtCAGGCGTGGTATGCGACGGCCTCGGCCTTACTAACTATTCCAACTTGGGCGGGGCGCCGCGCCGGCGCGCGGCGGTCGCGGCGCATGGGGTCGCTGTCGTTGATGGCGGCATGCATTGCCACAAGACAGACCa                                                            <  2:183356/137‑1 (MQ=18)
                                                            ttCAGGGGTGGATTTCGAAGGCCTCGGCCTAACGATCTATTCCAAGTTGGGCGGTGCGCCGCTCCGGCGCGCGGCGGTCGCGGCGCATGGGGTCGCTGTCGTTGATGGCGGCATGCATTGCCACAAGACAGACCAgg                                                          <  7:378999/136‑1 (MQ=21)
                                                            ggctggCGTTGTATTCGAAGGCCTCGGCCTTACGAACGATGCGACCTTGGGCGGGGCGCCGCGCCGGCGCGCGGCGGTCGCGGCGCATGGGGTCGCTGTCGTTGATGGCGGCATGCATTGCCACAAGACAGACCAgg                                                          <  2:350478/133‑1 (MQ=21)
                                                               cggCGTGGTATTCGAAGGCCTCGGCCTTACGAACGATTCGAACTTGGGCGGTGCGCCGCTCCGGCGCGCGGCGGTCGCGGCGCATGGGGTCGCTGTCGTTGATGGCGGCATGCATTGCCACAAGACAGACCAGGGCg                                                       <  5:394889/136‑1 (MQ=255)
                                                               cggCGTGGTATTCGAAGGCCTCGCCATAACCAACGATGCAAACTGGGGCGGGTGCCGGCGCCGGCGCGCGGCGGTCGCGGCGCATGGGGTCGCTGTCGTTGATGGCGGCATGCATTGCCACAAGACAGACCAGGGCg                                                       <  2:404201/136‑1 (MQ=21)
                                                                    gggTATTCGAAGGCCTCGGCCTTACGAACGATTCGACCTTGGGCGGGGCGCCGCGCCGGCGCGCGGCGGTCGCGGCGCATGGGGTCGCTGTCGTTGATGGCGGCATGCATTGCCACAAGACAGACCAGGGCGAGGCa                                                  <  3:158338/136‑1 (MQ=255)
                                                                      ttATTGGAAGGCCTGGGCCTTACGAACGATTCGAAGTTGGGGGGGGCGCCGCGCGGGCGCGGGGCGGTCGCGGCGCATGGGGTCGCTGTCGTTGATGGCGGCATGCATTGCCACAAGACAGACCAGGGCGAGGCaga                                                <  3:57573/136‑1 (MQ=25)
                                                                            gAAGGCCTCGGCCTTACGAACTATTCGATCTGGGGCGGTGCGCCGCTCCGGCGCGCGGCGGTCGCGGCGCATGGGGTCGCTGTCGTTGATGGCGGCATGCATTGCCACAAGACAGACCAGGGCGAGGCAGAGCGGAg                                          <  6:320933/137‑1 (MQ=255)
                                                                            gAAGGCCTCGGCCTTACGAACGATTCGAACTTGGGCGGTGCGCCGCTCCGGCGCGCGGCGGTCGCGGCGCATGGGGTCGCTGTCGTTGATGGCGGCATGCATTGCCACAAGACAGACCAGGGCGAGGCAGAGCGGAg                                          <  6:113947/137‑1 (MQ=255)
                                                                            gAAGGCCTCGGCCTTACGAACGATTCGAACTTGGGCGGTGCGCCGCGCCGGCGCGCGGCGGTCGCGGCGCATGGGGTCGCTGTCGTTGATGGCGGCATGCATTGCCACAAGACAGACCAGGGCGAGGCAGAGCGGAg                                          <  3:148262/137‑1 (MQ=255)
                                                                              gggacgcGGCCTTACGAACTATTCGAACTTGGGCGGGGCGCCGCTCCGGCGCGCGGCGGTCGCGGCGCATGGGGTCGCTGTCGTTGATGGCGGCATGCATTGCCACAAGACAGACCAGGGCGAGGCAGAGCGGAGtg                                        <  3:138577/131‑1 (MQ=255)
                                                                              cggCCTGGGACTTACGAACGATGCGAATTTGGGCGGTGCGCCGCTCCGGCGCGCGGCGGTCGCGGCGCATGGGGTCGCTGTCGTTGATGGCGGCATGCATTGCCACAAGACAGACCAGGGCGAGGCAGAGCGGAGtg                                        <  7:25104/136‑1 (MQ=255)
                                                                              cggCCGCGGCCTTACTAACTATTCGACCTTGGGCGGGCCGCCGCTCGGGCGCGCGGCGGTCGCGGCGCATGGGGTCGCTGTCGTTGATGGCGGCATGCATTGCCACAAGACAGACCAGGGCGAGGCAGAGCGGAGtg                                        <  2:403656/136‑1 (MQ=25)
                                                                              aGGCCTCGGCCTTACGAACGATTCGAACTTGGGCGGTGGGCAGCTCGGGCGCGCGGCGGTGGCGGCGCAGGGGGTCGCTGTCGTTGATGGCGGGATGCATTGCCACAAGACAGACCAGGGCGAGGCAGAGCGGAGtg                                        <  2:274095/137‑1 (MQ=25)
                                                                                 ccTCGGCCTTACGAACGATTCGAACTTGGGCGGTGCGCCGCTCCGGCGCGCGGCGGTCGCGGCGCATGGGGTCGCTGTCGTTGATGGCGGCATGCATTGCCACAAGACAGACCAGGGCGAGGCAGAGCGGAGtgatg                                     <  5:102048/137‑1 (MQ=255)
                                                                                 ccTCGGCCTTACGAACGATTCGAACTTGGGCGGTGCGCCGCTCCGGCGCGCGGCGGTCGCGGCGCATGGGGGCGCTGTCGTTGATGGCGGCATGCATTGCCACAAGACAGACCAGGGCGAGGCAGAGCGGAGtgatg                                     <  8:339808/137‑1 (MQ=255)
                                                                                    cGGCCTTACGAACGATGCTAACTTGGGCGGTGCGCCGCGCCGGCGCGCGGCGGTCGCGGCGCATGGGGTCGCTGTCGTTGATGGCGGCATGCATTGCCACAAGACAGACCAGGGCGAGGCAGAGCGGAGTGATGAtc                                  <  4:414774/137‑1 (MQ=255)
                                                                                                cgaTTCGAACTTGGGCGGGGCGCCGCGCCGGCGCGCGGCGGTCGCGGCGCATGGGGTCGCTGTCGTTGATGGCGGCATGCATTGCCACAAGACAGACCAGGGCGAGGCAGAGCGGAGTGATGATCTTCTGGTATgcc                      <  2:141083/137‑1 (MQ=255)
                                                                                                 tatccGAAGTGGGGCGGGCCGCCGCGCCGGCGCGCGGCGGGCGCGGCGCATGGGGTCGCTGTCGTTGATGGCGGCATGCATTGCCACAAGACAGACCAGGGCGAGGCAGAGCGGAGTGATGATCTTCTGGTATgccg                     <  1:342149/133‑1 (MQ=38)
                                                                                                           ggggCGGTGCGCCGCTCCGGCGCGCGGCGGTCGCGGCGCATGGGGTCGCTGTCGTTGATGGCGGCATGCATTGCCACAAGACAGACCAGGGCGAGGCAGAGCGGAGTGATGATCTTCTGGTATGCCGCCTCTGCTAc           <  2:308304/136‑1 (MQ=255)
                                                                                                             ggCGGGGCGCCGCGCCGGCGCGCGGCGGGCGCGGCGCATGGGGTCGCTGTCGTTGATGGCGGCATGCATTGCCACAAGACAGACCAGGGCGAGGCAGAGCGGAGTGATGATCTTCTGGTATGCCGCCTCTGCTACga         <  2:384177/137‑1 (MQ=255)
                                                                                                              ggcgtggcgcCGCGCGGGCGCGCGGCGGGCGCGGCGAATGGGGTCGCTGTCGTTGATGGCGGCATGCATTGCCACAAGACAGACCAGGGCGAGGCAGAGCGGAGTGATGATCTTCTGGTATGCCGCCTCTGCTACga         <  6:235683/131‑1 (MQ=255)
                                                                                                                     cgccgcgccgGCGCGCGGCGGTCGCGGCGCATGGGGTCGCTGTCGTTGATGGCGGCATGCATTGCCACAAGACAGACCAGGGCGAGGCAGAGCGGAGTGATGATCTTCTGGTATGCCGCCTCTGCTACGAGAGCAGt  <  5:415964/137‑1 (MQ=255)
                                                                                                                     cgccgcgccgGCGCGCGGCGGTCGCGGCGCATGGGGTCGCTGTCGTTGATGGCGGCATGCATTGCCACAAGACAGACCAGGGCGAGGCAGAGCGGAGTGATGATCTTCTGGTATGCCGCCTCTGCTACGAGAGCAGt  <  5:140874/137‑1 (MQ=255)
                                                                                                                           |                                                                                                                                  
TGTCCACCTTCTCTTGGAAGGCCGGGTGCATTGTTGTCTGCATGGATGATCCTCGGGTGGGTCAGGCGTGGTATTCGAAGGCCTCGGCCTTACGAACGATTCGAACTTGGGCGGTGCGCCGCTCCGGCGCGCGGCGGTCGCGGCGCATGGGGTCGCTGTCGTTGATGGCGGCATGCATTGCCACAAGACAGACCAGGGCGAGGCAGAGCGGAGTGATGATCTTCTGGTATGCCGCCTCTGCTACGAGAGCAGT  >  ACB122/1740285‑1740537

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: