Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1221,721,8690GC26.0% 87.0 / 11.8 50V98V (GTC→GTGCJ019_01617hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (30/7);  new base C (0/13);  total (30/20)
Fisher's exact test for biased strand distribution p-value = 2.18e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCACCCGGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACAGGCTGTTCCATTTCCCTACAGCCT  >  ACB122/1721733‑1721999
                                                                                                                                         |                                                                                                                                  
gCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCg                                                                                                                                    <  5:156023/137‑1 (MQ=255)
     tATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCa                                                                                                                               <  2:222039/137‑1 (MQ=255)
        tCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGcc                                                                                                                            >  5:120458/1‑137 (MQ=255)
        tCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGcc                                                                                                                            >  5:149211/1‑137 (MQ=255)
        tCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGcc                                                                                                                            >  5:170455/1‑137 (MQ=255)
        tCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGcc                                                                                                                            >  5:363104/1‑137 (MQ=255)
        tCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCACCGTCTTACTCGACCCAGcc                                                                                                                            >  8:139715/1‑137 (MQ=255)
        tCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGGCAGCATCTTATTCCCCGTCTGACTCGACCCCGcc                                                                                                                            >  7:127463/1‑137 (MQ=255)
        tCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGAAGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGcc                                                                                                                            >  5:305284/1‑137 (MQ=255)
             ggTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCAcccc                                                                                                                       >  7:94314/1‑136 (MQ=255)
              gTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATTATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGCCCCAGCGGTCAGCATCTTCTTACCCCTCTTCCTCGACCCCGCCAccccg                                                                                                                      >  3:16240/1‑135 (MQ=255)
              gTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCACCCgg                                                                                                                      >  3:364127/1‑137 (MQ=255)
               aaaaCCAATTCATACAGAAACTAGAAAGCCTGCGTTTTAATGTTGTACGCCAGGTGTTTTCCTGTTCACGGTCGGCGCTGATGCGGACCCTGCGGTCAGCTTCTTCTTCCCCGTCTTCCTCGACCCCGCCACCCggg                                                                                                                     <  4:62787/136‑2 (MQ=255)
                  aCCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCACCCGgcgcg                                                                                                                  >  5:326921/1‑137 (MQ=255)
                         ccTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCCGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGCCCCCGCCCCCCGGcgcgccgcgcc                                                                                                           >  2:380310/1‑136 (MQ=255)
                           tACAGAAACGAGCAAGCCAGCCGTTTCTTGATGTAAGCCAGGTGAATTCCTCTTCACCGTCGGCGCAGATTCGGACCCAGCCGTGAGGAACTTCTTCCCCCTCTTACTCGACCCATCCACCCGGCCCGCCGCCCAcg                                                                                                         >  5:254092/1‑136 (MQ=255)
                                aaaCGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCACCCGGCGCGCAGCGCACTTCcac                                                                                                    >  1:52263/1‑137 (MQ=255)
                                     aGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCACCCGGCGCGCAGCGCACTTCcccaccca                                                                                               >  8:40435/1‑137 (MQ=255)
                                     aGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCATCGACACATCCACCCGCCGCGCCGCGCACTTAAACAccaa                                                                                               >  1:43162/1‑135 (MQ=255)
                                     aGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCAACTTCTTACCCGTCTTCCTCGACCCCGCCACCCCGCGCGCAGCGCACTTACACAccag                                                                                               >  1:330009/1‑135 (MQ=255)
                                        aaGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGCCCCAGCCACCCGGCGCGCAGCGCCCTTCCACACCCAgcg                                                                                            >  2:150091/1‑137 (MQ=255)
                                        aaGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCCGCCACCCGGCGCGCAGCGCACTTCCACACCCAgcc                                                                                            >  5:286343/1‑136 (MQ=255)
                                                       tgatgTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTACCCCGTCTTCCTCGACCCAGCCACCCGGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGgcgc                                                                             >  2:89294/1‑137 (MQ=255)
                                                               gCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCACCCGGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAg                                                                     >  1:44487/1‑137 (MQ=255)
                                                                cccggTCATTCCCTCTCCCCGGGCGGGGATGATGCGGACCCAGGGGTCAGCATCTTCTGCCCCGTCTACCTCGCCCCACCCACCCGGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAgg                                                                    <  1:250437/134‑1 (MQ=255)
                                                                      taTTCCCTCTTCACGGGGGGCGGTGTGGCGCCCCCAGGGGTCACATGGTTCTCCCCCGCCTCCCTCCCCCCAGCCCCCCGGCGCCCCGCCCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGt                                                              <  4:88838/136‑1 (MQ=255)
                                                                           ccTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCACCCGGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCg                                                         >  5:394115/1‑137 (MQ=255)
                                                                                 cACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCACCCGGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCAc                                                   >  6:372819/1‑137 (MQ=255)
                                                                                    cGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCACCCGGCGCGCCGCGCACTTCCACCCCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACggg                                                >  8:59284/1‑136 (MQ=255)
                                                                                    cGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCACCCGGCGCGCAGCGCACTTCCACACCCAGCCCGGCAAATGCGGCGCGGGCGCCGGCAACCTCATCGGCCCACCGc                                                >  1:171612/1‑137 (MQ=255)
                                                                                           gcTGATGCGGACCCCGCGGTCAGCATCTTCTTCCCCGTCTTCCTCCCCCCCCCCCCCCGGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGt                                         <  5:386319/137‑1 (MQ=255)
                                                                                               ttggggACCCAGCGGTCTGCTTCTTATTCCCCGTCTTCCTCCACCCAGCCCCCCGGCGCGCCGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCt                                     <  6:326921/133‑1 (MQ=255)
                                                                                               tggcgcACCCCGCGGTCAGATTGTTCTACCCGGTCTTCCTCCACCCCCCCCCCCGGCGCGCCGCCCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCt                                     <  4:189800/135‑1 (MQ=255)
                                                                                               tgccGGCCCCCGCGGTCAGCATGTTCTCCCCCGTCTCCCCCCCCCCCCCCCCCCGGCGCGCCGCCCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCt                                     <  3:359209/134‑1 (MQ=255)
                                                                                               ctGCGGACCCCGCGGTCATCATCTTCTGCCCCGTCTTCCTCGCCCCCGCCCCCCGGCGCGCCCCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCt                                     <  6:316923/136‑1 (MQ=255)
                                                                                               aggCGGCCCCCGGGGCCACCATTTTCTCCCCCGGCTCCCCCCCCCCCCCCCCCCGGCGCGCACCCCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCt                                     <  3:353543/135‑1 (MQ=255)
                                                                                                 gCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCACCCGGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCttt                                   >  7:69017/1‑137 (MQ=255)
                                                                                                            ggTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCACCCGGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCCCGGCGTCTGGTCGCTGTCCAGATACCgg                        >  6:337482/1‑137 (MQ=255)
                                                                                                            ggTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCACCCGGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCATATACAgg                        >  3:133294/1‑137 (MQ=255)
                                                                                                            ggTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCACCCGGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACAgg                        >  3:139250/1‑137 (MQ=255)
                                                                                                            ggTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCACCCGGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGAGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCGTTCCAGATACAgg                        >  6:332736/1‑137 (MQ=255)
                                                                                                                 gaaTCTTCTTCCCCGTCTTCCTCCACCCACCCACCCGGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACAGGCTGtt                   <  8:94314/135‑1 (MQ=255)
                                                                                                                  attcttctTCCCGTTCTCCCTCGCCCCCCCCCCCGGGGGCGCCGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACAGGCTGTTc                  <  5:302884/135‑1 (MQ=255)
                                                                                                                   aTCTGCTCCCCCGTCTACCTCCCCCCCCCCCCCCGGCGCGCCGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACAGGCTGTTcc                 <  8:42200/137‑1 (MQ=255)
                                                                                                                    tcttctCCCCCGTCTCCCCCCCCCCCCCCCCCCGGCGCGCCGCCCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACAGGCTGTTCCa                <  4:259714/137‑1 (MQ=255)
                                                                                                                         cttcCCCGTCTTCCTCGACCCAGCCACCCGGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACAGGCTGTTCCATTTc            <  6:363104/137‑1 (MQ=255)
                                                                                                                          ttcCCCTTCTTCCTCCACCCCCCCCCCCGGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACAGGCTGTTCCATTTcc           <  4:170503/137‑1 (MQ=255)
                                                                                                                          ttcCCCGTCTTCCTCGACCCAGCCCCCCGGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACAGGCTGTTCCATTTcc           <  2:44487/137‑1 (MQ=255)
                                                                                                                          taccccGCCTCCCACCACCCCGCCACCCGGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACAGGCTGTTCCATTTcc           <  6:120458/135‑1 (MQ=255)
                                                                                                                               ggTCTTCCTCCCCCCACCCCCCGGGCGGCCCGCCCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACAGGCTGTTCCATTTCCCTACa      <  4:16240/136‑1 (MQ=255)
                                                                                                                                   ttCCTCCACCCCCCCACCCGGCGCCCCGCCCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACAGGCTGTTCCATTTCCCTACAGCCt  <  1:286626/137‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                  
GCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCACCCGGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACAGGCTGTTCCATTTCCCTACAGCCT  >  ACB122/1721733‑1721999

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: