Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1221,348,4200CG32.0% 67.2 / 14.3 50D4H (GAC→CAC) CJ019_01261hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (9/25);  new base G (15/1);  total (24/26)
Fisher's exact test for biased strand distribution p-value = 8.13e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ATCTTGTCGTCGGTGTAAGGGACGACAGTGCCGTTACGCTTGATCACGCGCAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCT‑CCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCA‑CCGTTCCGTCCACTTGCTCGATAACCTCGCGGGGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGA  >  ACB122/1348284‑1348555
                                                                                                                                        |                                                                                                                                         
aTCTTGTCGTCGGTGTAAGGGACGACAGTGCCGTTACGCTTGATCACGCGCAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTc                                                                                                                                           <  6:395281/137‑1 (MQ=255)
    tGTCGTCGGTGTAAGGGACGACAGTGCCGTTACGCTTGATCACGCGCAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGtt                                                                                                                                       <  8:211274/137‑1 (MQ=255)
    tGTCGTCGGTGTAAGGGACGACAGTGCCGTTACGCTTGATCACGCGCAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGtt                                                                                                                                       <  7:176734/137‑1 (MQ=255)
     gtcgtcGGTGTAAGGGACGACAGTGCCGTTACGCTTGATCACGCGCAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTGGGgtt                                                                                                                                      >  5:227414/1‑137 (MQ=255)
             tgtAAGGGACGACAGTGCCGTTACGCTTGATCACGCGCAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCAtgggt                                                                                                                              <  1:136507/137‑1 (MQ=255)
             tgtAAGGGACGACAGTGCCGTTACGCTTGATCACGCGCAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGGGGGGTTGCAtgggg                                                                                                                              >  8:25412/1‑136 (MQ=255)
              gtAAGGGACGACAGTGCCGTTACGCTTGATCACGCGCAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCAtgggtg                                                                                                                             <  7:291526/137‑1 (MQ=255)
                      cgacAGTGCCGTTACGCTTGATCACGCGCAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGtct‑c                                                                                                                    <  7:187165/137‑1 (MQ=255)
                      cgacAGTGCCGTTACGCTTGATCACGCGCAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGtct‑c                                                                                                                    <  3:227657/137‑1 (MQ=255)
                      cgacAGTGCCGTTACGCTTGATCACGCGCAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGGTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGtct‑c                                                                                                                    <  5:143549/137‑1 (MQ=255)
                          aGTGCCGTTACGCTTGATCACGCGCAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGGGGGTGTGGGTGGGGGGGTGtcggcaaa                                                                                                                 >  1:205779/1‑131 (MQ=255)
                          aGTGCCGTTACGCTTGATCACGCGCAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGGGGGGTTGGATGGGGGGGTGTGT‑Cccca                                                                                                                >  3:249739/1‑137 (MQ=255)
                           gTGCCGTTACGCTTGATCACGCGCAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCT‑CCACAg                                                                                                               <  4:181202/137‑1 (MQ=255)
                           gTGCCGTTACGCTTGATCACGCGCAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCT‑CCACAg                                                                                                               <  5:131071/137‑1 (MQ=255)
                              ccGTTACGCTTGATCACGCGCAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGGTTGAATAGGTGGGTGTGTGCCACAGtt                                                                                                             >  2:314530/1‑137 (MQ=255)
                                 ttACGCTTGATCACGCGCAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGGATGGGTTGGTGTAT‑CAACAGGTTGta                                                                                                         >  4:289241/1‑137 (MQ=255)
                                 ttACGCTTGATCACGCGCAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCT‑CCCCAGtttatt                                                                                                         >  2:313573/1‑134 (MQ=255)
                                      cTTGATCACGCGCAGTTGGACGGGTGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGGCGGGTTGCATTGGTGGGTGTCT‑CCCCAGTTTAGATAttt                                                                                                    >  6:342396/1‑136 (MQ=255)
                                         taTCACGCGCAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCT‑CCACAGTTTCTATATTGTTc                                                                                                 <  3:289241/136‑1 (MQ=255)
                                         gATGACGCGCAGTGGGCCGGGGGCGGTGGCTGCCAGATCCTGGTTGGAATGGGCGGCCTGCGGCACCTGGGCCTGCGTGTTCTCGCGAGTCGTGTCGGATTGCATGGGTGGGTGTCT‑CCACAGTTTCTATATTGTTc                                                                                                 <  8:70622/137‑1 (MQ=255)
                                         gATCACGCGCAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCT‑CCACAGTTTCTATATTGTTc                                                                                                 <  3:365792/137‑1 (MQ=255)
                                         gATCACGCGCAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCT‑CCACAGTTTCTATATTGTTc                                                                                                 <  1:240548/137‑1 (MQ=255)
                                         gATCACGCGCAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCT‑CCACAGTTTCTATATTGTTc                                                                                                 <  7:347681/137‑1 (MQ=255)
                                         gATCACGCGCAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGAGGGTTCTCGCGAGTAGTGTGTGTTTTCATTGGTTTGTGTGT‑CCACAGTTTATTTATttttt                                                                                                 >  5:273013/1‑136 (MQ=255)
                                           tCACGCGCAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGGGGGGTTGGATGGGTGGGTGTGT‑GGCCAGTTTGGATATTTTTTAg                                                                                               >  2:183490/1‑137 (MQ=255)
                                            cACGCGCAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGTTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCT‑CCACAGTTTCTATATTGTTAAgg                                                                                              >  8:384772/1‑137 (MQ=255)
                                            cACGCGCAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTGGGTTTGTGTGGGTGGGGGGCT‑GGACACTTTGTATAATTTTTAgg                                                                                              >  3:12521/1‑137 (MQ=255)
                                            cACGCGCAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGGGGGGTTTGATGGGGGGGGGTGT‑CCACACTTTGTATATTGTTAAgg                                                                                              >  7:174307/1‑137 (MQ=255)
                                                 gcAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGGGGGTTTGGATGGGTGGGTGTGT‑GCACAGATTTGTTATTTTTTATGCGccg                                                                                         >  3:331550/1‑136 (MQ=255)
                                                  cAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGGTTGCATGGGTGGGTGTCT‑CCACAGTTTCGATATTGTTAAGGCGCCtt                                                                                        >  6:189360/1‑137 (MQ=255)
                                                   aGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGGGGGGTTGGATGGGTGGGGGTGT‑GGCCAGTTTTTATTTTGTTAAGgcggcgtt                                                                                       >  1:316565/1‑132 (MQ=255)
                                                        gCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGGCGGTTTGGATGGGTGGGTGTCT‑CCACAGTTTATATATTGTTAAGGCGCCTTATgggg                                                                                  >  6:207287/1‑136 (MQ=255)
                                                                cGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGGCGGGTTGGATGGGTGGGTGGGT‑CCACACGTTCTATATTTTTCATGGGCCGTTTTGGCGGCcc‑ccc                                                                         >  4:78856/1‑136 (MQ=255)
                                                                cGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGGCGGGTTGGATGGGGGGGTGGGT‑GCACAGGTTGTATATTGTTGAGGGGGCGTTTTGGCGCCCA‑ccc                                                                         >  4:231536/1‑136 (MQ=255)
                                                                  gTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGGGGGTTTGGATGGGTGGGGGTGT‑GCACAGTTTGTATTTTTTTCAGGCGCCTTATTGGCGGCCC‑CCCtt                                                                       >  6:273420/1‑137 (MQ=255)
                                                                      cAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGGGGGTTTGAATGGGTGGGTGTGT‑CCACAGTTTGTATATTTTTCAGGCGGCGTATTGGCGCCCC‑CCCTGCCgg                                                                   >  6:124426/1‑136 (MQ=255)
                                                                        gCCAGATCCTGGTTGGAATCGGCGTCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCT‑CCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCA‑CCGTTCCGTcc                                                                 <  7:295494/137‑1 (MQ=255)
                                                                               ccTGGTTGTAATCGGCGGCGGGGGGCCCGTGGGCGTGGTGGTTCTCGCTAGTGGTGTGGGTTGGCATGGGTGGGTGTCT‑CCACAGTTTCTATATGTTTCAGGCGCCTTATGGGCGCCCA‑CCGTTCCGTCCACTTGCt                                                          <  8:255352/137‑1 (MQ=255)
                                                                                cTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTGGGTTTTGATGGGTGGGGGTGT‑CCACAGTTTCTATTTTGTTTAGGGGGCGTATGGGCGCCCC‑CCGGTGCGGGCACTTGcgg                                                         >  3:87057/1‑135 (MQ=255)
                                                                                cTGGTGGGACTCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCT‑CCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCA‑CCGTTCCGTCCACTTGCTc                                                         <  6:104248/137‑1 (MQ=255)
                                                                                  ggTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCT‑CCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCA‑CCGTTCCGTCCACTTGCTCGa                                                       <  5:180584/137‑1 (MQ=255)
                                                                                       taaTCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGGGTCT‑CCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCA‑CCGTTCCGTCCACTTGCTCGATAAcc                                                  <  3:67821/136‑1 (MQ=255)
                                                                                         aTCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCT‑CCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCA‑CCGTTCCGTCCACTTGCTCGATAACCTc                                                <  5:280641/137‑1 (MQ=255)
                                                                                         aTCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCT‑CCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCA‑CCGTTCCGTCCACTTGCTCGATAACCTc                                                <  7:59759/137‑1 (MQ=255)
                                                                                           cggcggCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGGGGGGTTGGATGGGGGGGGGTGT‑CCACCCTTTATATATTTTTAAGGCGCCTTTTGGGCGCCCACCCGTTCGGACC‑CTTGCTCGATAAccaacc                                              >  4:143585/1‑133 (MQ=255)
                                                                                                 cgtGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCT‑CCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCA‑CCGTTCCGTCCACTTGCTCGATAACCTCGCGGGGAt                                        <  4:58396/135‑1 (MQ=255)
                                                                                                 ccggCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGTTTTGCATGGGTGGGTGTCT‑CCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCA‑CCGTTCCGTCCACTTGCTCGATAACCTCGCGGGGAt                                        <  8:244079/134‑1 (MQ=255)
                                                                                                      ggCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTGGGTTTGGTTGGGGGGGTGTGT‑CCCCAGTGTGTATTTTGTTGAGGCGCGTGATGGGGGGCCC‑CCGGTCCGGCCCCTTTCGCGATAACCTCGCGGGGGTGGGCa                                   >  4:232479/1‑137 (MQ=255)
                                                                                                                        ttCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCT‑CCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCA‑CCGTTCCGTCCACTTGCTCGATAACCTCGCGGGGATGCGCAATGCATGCGCATccgccc                 <  5:189360/137‑1 (MQ=255)
                                                                                                                        ttCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCT‑CCACAGTTTCTATATTGTGCAGGCGCCTTATGGGCGCCCA‑CCGTTCCGTCCACTTGCTCGATAACCTCGCGGGGATGCGCAATGCATGCGCATccgccc                 <  7:25412/137‑1 (MQ=255)
                                                                                                                                aGTCGTGTCGGTTTGCATGGGTGGGTGTCT‑CCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCA‑CCGTTCCGTCCACTTGCTCGATAACCTCGCGGGGATGCGCAATGCATGCGCATCCGCCCGCGCGGGc         <  8:174307/137‑1 (MQ=255)
                                                                                                                                       gcGGTTTGCATGGGTGGGTGGCT‑CCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCA‑CCGTTCCGTCCACTTGCTCGATAACCTCGCGGGGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACTGa  <  1:313573/136‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                         
ATCTTGTCGTCGGTGTAAGGGACGACAGTGCCGTTACGCTTGATCACGCGCAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCT‑CCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCA‑CCGTTCCGTCCACTTGCTCGATAACCTCGCGGGGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGA  >  ACB122/1348284‑1348555

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: ‑