Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273334,2650AG54.8% ‑6.6 / 46.4 31H77R (CAT→CGT) yahGDUF1116 domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (7/7);  new base G (13/4);  total (20/11)
Fisher's exact test for biased strand distribution p-value = 1.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.08e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

ACTCAACTCGACTGGACGCCGCCGGGGCAGGGCAATATGCAGGTGGTGCAGGCGCTGGATAACATTGCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAGGCGATTAACGTGGTGCCGGGCATGACGGCGAAAACCATTCTTCACGCCGGGCCGCCGATCACCTGGGAAAAAATGTGCGGCGCGATGAAAGGCGCGGTCACCGGAGCGCTGG  >  NZ_CP009273/334129‑334402
                                                                                                                                        |                                                                                                                                         
aCTCAACTCGACTGGACGCCGCCGGGGCAGGGCAATATGCAGGTGGTGCAGGCGCTGGATAACATTGCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCGTc                                                                                                                                         >  2:190423/1‑139 (MQ=255)
     aCTCGACTGGACGCCGCCGGGGCAGGGCAATATGCAGGTGGTGCAGGCGCTGGATAACATTGCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCGTCCGGTg                                                                                                                                    <  1:219448/139‑1 (MQ=255)
      cTCGACTGGACGCCGCCGGGGCAGGGCAATATGCAGGTGGTGCAGGCGCTGGATAACATTGCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCGTCCGGTGc                                                                                                                                   >  2:34918/1‑139 (MQ=255)
                     ccgGGGCAGGGCAATATGCAGGTGGTGCAGGCGCTGGATAACATTGCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATc                                                                                                                    <  2:375469/139‑1 (MQ=255)
                      cgGGGCAGGGCAATATGCAGGTGGTGCAGGCGCTGGATAACATTGCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCGTCCGGTGCTGATTGGTTTTGATCa                                                                                                                   <  1:124858/139‑1 (MQ=255)
                            agggcaATATGCAGGTGGTGCAGGCGCTGGATAACATTGCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAGGCGAt                                                                                                             <  1:272357/139‑1 (MQ=255)
                             gggcaATATGCAGGTGGTGCAGGCGCTGGATAACATTGCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCGTCCGGTGCTGATTGGTTTTGATCAGGCGAt                                                                                                             >  1:72352/1‑138 (MQ=255)
                                   tatGCAGGTGGTGCAGGCGCTGGATAACATTGCCGATTCGCCGCTGGCGGACAAAATCGACGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCGTCCGGTGCTGATTGGTTTTGATCAGGCGATTAACgtg                                                                                                      >  1:203006/1‑139 (MQ=255)
                                        aGGTGGTGCAGGCGCTGGATAACATTGCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCGTCCGGTGCTGATTGGTTTTGATCAGGCGATTAACGTGGTGcc                                                                                                 >  1:331072/1‑139 (MQ=255)
                                        aGGTGGTGCAGGCGCTGGATAACATTGCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAGGCGATTAACGTGGTGcc                                                                                                 >  1:172631/1‑139 (MQ=255)
                                             gtgCAGGCGCTGGATAACATTGCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAGGCGATTAACGTGGTGCCGGGCa                                                                                            >  2:302599/1‑139 (MQ=255)
                                              tgCAGGCGCTGGATAACATTGCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCGTCCGGTGCTGATTGGTTTTGATCAGGCGATTAACGTGGTGCCGGGCAt                                                                                           >  2:332418/1‑139 (MQ=255)
                                              tgCAGGCGCTGGATAACATTGCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCGTCCGGTGCTGATTGGTTTTGATCAGGCGATTAACGTGGTGCCGGGCAt                                                                                           >  2:325593/1‑139 (MQ=255)
                                                    cgcTGGATAACATTGCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCGTCCGGTGCTGATTGGTTTTGATCAGGCGATTAACGTGGTGCCGGGCATGACGGc                                                                                     >  1:427908/1‑139 (MQ=255)
                                                        ggATAACATTGCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAGGCGATTAACGTGGTGCCGGGCATGACGGCGaa                                                                                  >  1:365215/1‑138 (MQ=255)
                                                                   ccGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCGTCCGGTGCTGATTGGTTTTGATCAGGCGATTAACGTGGTGCCGGGCATGACGGCGAAAACCATTCTTCa                                                                      >  1:317821/1‑139 (MQ=255)
                                                                   ccGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCGTCCGGTGCTGATTGGTTTTGATCAGGCGATTAACGTGGTGCCGGGCATGACGGCGAAAACCATTCTTCa                                                                      >  2:337728/1‑139 (MQ=255)
                                                                   ccGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCa                                                                                                                   >  2:106542/1‑94 (MQ=255)
                                                                   ccGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCa                                                                                                                   <  1:106542/94‑1 (MQ=255)
                                                                   ccGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAGGCGATTAACGTGGTGCCGGGCATGACGGCGAAAACCATTCTTCa                                                                      >  2:397572/1‑139 (MQ=255)
                                                                      aTTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCGTCCGGTGCTGATTGGTTTTGATCAGGCGATTAACGTGGTGCCGGGCATGACGGCGAAAACCATTCTTCACGc                                                                   >  2:145202/1‑139 (MQ=255)
                                                                                 gCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCGTCCGGTGCTGATTGGTTTTGATCAGGCGATTAACGTGGTGCCGGGCATGACGGCGAAAACCATTCTTCACGCCGGGCCGCCGa                                                        <  2:203006/139‑1 (MQ=255)
                                                                                  cGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAGGCGATTAACGTGGTGCCGGGCATGACGGCGAAAACCATTCTTCACGCCGGGCCGCCGAt                                                       <  1:302599/139‑1 (MQ=255)
                                                                                                                         gTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAGGCGATTAACGTGGTGCCGGGCATGACGGCGAAAACCATTCTTCACGCCGGGCCGCCGATCACCTGGGAAAAAATGTGCGGCGCGATGAAAGGCGCGGt                <  2:170242/139‑1 (MQ=255)
                                                                                                                          tattatCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAGGCGATTAACGTGGTGCCGGGCATGACGGCGAAAACCATTCTTCACGCCGGGCCGCCGATCACCTGGGAAAAAATGTGCGGCGCGATGAAAGgc                     <  2:38969/133‑1 (MQ=255)
                                                                                                                          tattatCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAGGCGATTAACGTGGTGCCGGGCATGACGGCGAAAACCATTCTTCACGCCGGGCCGCCGATCACCTGGGAAAAAATGTGCGGCGCGATGAAAGgc                     >  1:38969/1‑133 (MQ=255)
                                                                                                                               tCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAGGCGATTAACGTGGTGCCGGGCATGACGGCGAAAGCCATTCTTCACGCCGGGCCGCCGATCACCTGGGAAAAAATGTGCGGCGCGATGAAAGGCGCGGTCACCg           <  1:397572/138‑1 (MQ=255)
                                                                                                                                    tCGCGTCCGGTGCTGATTGGTTTTGATCAGGCGATTAACGTGGTGCCGGGCATGACGGCGAAAACCATTCTTCa                                                                      >  1:115568/1‑74 (MQ=255)
                                                                                                                                    tCGCGTCCGGTGCTGATTGGTTTTGATCAGGCGATTAACGTGGTGCCGGGCATGACGGCGAAAACCATTCTTCa                                                                      <  2:115568/74‑1 (MQ=255)
                                                                                                                                    tCGCGTCCGGTGCTGATTGGTTTTGATCAGGCGATTAACGTGGTGCCGGGCATGACGGCGAAAACCATTCTTCACGCCGGGCCGCCGATCACCTGGGAAAAAATGTGCGGCGCGATGAAAGGCGCGGTCACCGGAgcgc     >  2:404695/1‑139 (MQ=255)
                                                                                                                                       cATCCGGTGCTGATTGGTTTTGATCAGGCGATTAACGTGGTGCCGGGCATGACGGCGAAAACCATTCTTCACGCCGGGCCGCCGATCACCTGGGAAAAAATGTGCGGCGCGATGAAAGGCGCGGTCACCGGAGCGCTgg  >  1:235418/1‑139 (MQ=255)
                                                                                                                                        |                                                                                                                                         
ACTCAACTCGACTGGACGCCGCCGGGGCAGGGCAATATGCAGGTGGTGCAGGCGCTGGATAACATTGCCGATTCGCCGCTGGCGGACAAAATCGCCGCCGCTAACCAGCAGGCGCTGGAGCGTATTATCCAGTCGCATCCGGTGCTGATTGGTTTTGATCAGGCGATTAACGTGGTGCCGGGCATGACGGCGAAAACCATTCTTCACGCCGGGCCGCCGATCACCTGGGAAAAAATGTGCGGCGCGATGAAAGGCGCGGTCACCGGAGCGCTGG  >  NZ_CP009273/334129‑334402

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: