Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273289,2560TC54.5% ‑1.2 / 28.9 22D125G (GAC→GGC) yagLrecombinase family protein
Reads supporting (aligned to +/- strand):  ref base T (4/6);  new base C (3/9);  total (7/15)
Fisher's exact test for biased strand distribution p-value = 6.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

ATGTCAATCGTATCCTGACGAGCAATCAAAAAAGACAAGCAACTTTCAGCAGAAGGCTCTATCCTTGAAGCTGGATGGTAGCTGTAGATAAACATTTCTTTTTGAAGGCAAAAAAATAAAATATTCTGTAAGTCCTCCACAGTTCTGCCAAGACAGCTTAGATCTGGAATGACAAGCGTGTCCTGCCGGTTCATTTCATGATGAATCAATTGTCTGAGTAAACGTTGCTCTGAAACAAC  >  NZ_CP009273/289124‑289362
                                                                                                                                    |                                                                                                          
aTGTCAATCGTATCCTGACGAGCAATCAAAAAAGACAAGCAACTTTCAGCAGAAGGCTCTATCCTTGAAGCTGGATGGTAGCTGTAGATAAACATTTCTTTTTGAAGGCAAAAAAATAAAATATTCTGTAAGTCCTCca                                                                                                      <  2:217459/139‑1 (MQ=255)
         gTATCCTGACGAGCAATCAAAAAAGACAAGCAACTTTCAGCAGAAGGCTCTATCCTTGAAGCTGGATGGTAGCTGTAGATAAACATTTCTTTTTGAAGGCAAAAAAATAAAATATTCTGTAAGCCCTCCACAGTTCTGc                                                                                             >  1:79606/1‑139 (MQ=255)
              cTGACGAGCAATCAAAAAAGACAAGCAACTTTCAGCAGAAGGCTCTATCCTTGAAGCTGGATGGTAGCTGTAGATAAACATTTCTTTTTGAAGGCAAAAAAATAAAATATTCTGTAAGCCCTCCACAGTTCTGCCAAGa                                                                                        <  2:79606/139‑1 (MQ=255)
              cTGACGAGCAATCAAAAAAGACAAGCAACTTTCAGCAGAAGGCTCTATCCTTGAAGCTGGATGGTAGCTGTAGATAAACATTTCTTTTTGAAGGCAAAAAAATAAAATATTCTGTAAGCCCTCCACAGTTCTGCCAAGa                                                                                        <  1:356837/139‑1 (MQ=255)
                 aCGAGCAATCAAAAAAGACAAGCAACTTTCAGCAGAAGGCTCTATCCTTGAAGCTGGATGGTAGCTGTAGATAAACATTTCTTTTTGAAGGCAAAAAAATAAAATATTCTGTAAGTCCTCCACAGTTCTGCCAAGACAg                                                                                     >  1:222296/1‑139 (MQ=255)
                 aCGAGCAATCAAAAAAGACAAGCAACTTTCAGCAGAAGGCTCTATCCTTGAAGCTGGATGGTAGCTGTAGATAAACATTTCTTTTTGAAGGCAAAAAAATAAAATATTCTGTAAGCCCTCCACAGTTCTGCCAAGACAg                                                                                     >  1:486300/1‑139 (MQ=255)
                        aTCAAAAAAGACAAGCAACTTTCAGCAGAAGGCTCTATCCTTGAAGCTGGATGGTAGCTGTAGATAAACATTTCTTTTTGAAGGCAAAAAAATAAAATATTCTGTAAGTCCTCCACAGTTCTGCCAAGACAGCt                                                                                   >  2:96419/1‑134 (MQ=255)
                               aaGACAAGCAACTTTCAGCAGAAGGCTCTATCCTTGAAGCTGGATGGTAGCTGTAGATAAACATTTCTTTTTGAAGGCAAAAAAATAAAATATTCTGTAAGTCCTCCACAGTTCTGCCAAGACAGCTTAGATCTGGAAt                                                                       <  2:222296/139‑1 (MQ=255)
                                     aGCAACTTTCAGCAGAAGGCTCTATCCTTGAAGCTGGATGGTAGCTGTAGATAAACATTTCTTTTTGAAGGCAAAAAAATAAAATATTCTGTAAGTCCTCCACAGTTCTGCCAAGACAGCTTAGATCTGGAATGACAAg                                                                 <  2:279021/139‑1 (MQ=255)
                                                   gAAGGCTCTATCCTTGAAGCTGGATGGTAGCTGTAGATAAACATTTCTTTTTGAAGGCAAAAAAATAAAATATTCTGTAAGCCCTCCACAGTTCTGCCAAGACAGCTTAGATCTGGAATGACAAGCGTGTCCTGCCGGt                                                   <  1:421473/139‑1 (MQ=255)
                                                               cTTGAAGCTGGATGGTAGCTGTAGATAAACATTTCTTTTTGAAGGCAAAAAAATAAAATATTCTGTAAGTCCTCCACAGTTCTGCCAAGACAGCTTAGATCTGGAATGACAAGCGTGTCCTGCCGGTTCATTTCatga                                        <  1:96419/138‑1 (MQ=255)
                                                                      cTGGATGGTAGCTGTAGATAAACATTTCTTTTTGAAGGCAAAAAAATAAAATATTCTGTAAGCCCTCCACAGTTCTGCCAAGACAGCTTAGATCTGGAATGACAAGCGTGTCCTGCCGGTTCATTTCATGATGAATCaa                                <  2:264834/139‑1 (MQ=255)
                                                                          aTGGTAGCTGTAGATAAACATTTCTTTTTGAAGGCAAAAAAATAAAATATTCTGTAAGCCCTCCACAGTTCTGCCAAGACAGCTTAGATCTGGAATGACAAGCGTGTCCTGCCGGTTCATTTCATGATGAATCAATTGt                            <  1:158496/139‑1 (MQ=255)
                                                                                          aaCATTTCTTTTTGAAGGCAAAAAAATAAAATATTCTGTAAGCCCTCCACAGTTCTGCCAAGACAGCTTAGATCTGGAATGACAAGCGTGTCCTGCCGGTTCATTTCatgatg                                      <  1:349724/113‑1 (MQ=255)
                                                                                          aaCATTTCTTTTTGAAGGCAAAAAAATAAAATATTCTGTAAGCCCTCCACAGTTCTGCCAAGACAGCTTAGATCTGGAATGACAAGCGTGTCCTGCCGGTTCATTTCatgatg                                      >  2:349724/1‑113 (MQ=255)
                                                                                          aaCATTTCTTTTTGAAGGCAAAAAAATAAAATATTCTGTAAGCCCTCCACAGTTCTGCCAAGACAGCTTAGATCTGGAATGACAAGCGTGTCCTGCCGGTTCATTTCATGATGAATCAATTGTCTGAGTAAACGTTGct            <  2:160247/139‑1 (MQ=255)
                                                                                            cATTTCTTTTTGAAGGCAAAAAAATAAAATATTCTGTAAGTCCTCCACAGTTCTGCCAAGACAGCTTAGATCTGGAATGACAAGCGTGTCCTGc                                                       >  1:370559/1‑94 (MQ=255)
                                                                                            cATTTCTTTTTGAAGGCAAAAAAATAAAATATTCTGTAAGTCCTCCACAGTTCTGCCAAGACAGCTTAGATCTGGAATGACAAGCGTGTCCTGc                                                       <  2:370559/94‑1 (MQ=255)
                                                                                              tttCTTTTTGAAGGCAAAAAAATAAAATATTCTGTAAGTCCTCCACAGTTCTGCCAAGACAGCTTAGATCTGGAATGACAAGCGTGTCCTGCCGGTTCATTTCATGATGAATCAATTGTCTGAGTAAACGTTGCTCTg         >  1:282310/1‑138 (MQ=255)
                                                                                                     ttGAAGGCAAAAAAATAAAATATTCTGTAAGTCCTCCACAGTTCTGCCAAGACAGCTTAGATCTGGAATGACAAGCGTGTCCTGCCGGTTCATTTCATGATGAATCAATTGTCTGAGTAAACGTTGCTCTGAaacaac  <  1:307307/138‑1 (MQ=255)
                                                                                                     ttGAAGGCAAAAAAATAAAATATTCTGTAAGCCCTCCACAGTTCTGCCAAGACAGCTTAGATCTGGAATGACAAGCGTGTCCTGCCGGTTCATTTCATGATGAATCAATTGTCTGAGTAAACGTTGCTCTGAaacaac  <  1:458590/138‑1 (MQ=255)
                                                                                                     ttGAAGGCAAAAAAATAAAATATTCTGTAAGCCCTCCACAGTTCTGCCAAGACAGCTTAGATCTGGAATGACAAGCGTGTCCTGCCGGTTCATTTCATGATGAATCAATTGTCTGAGTAAACGTTGCTCTGAaacaac  <  2:486300/138‑1 (MQ=255)
                                                                                                                                    |                                                                                                          
ATGTCAATCGTATCCTGACGAGCAATCAAAAAAGACAAGCAACTTTCAGCAGAAGGCTCTATCCTTGAAGCTGGATGGTAGCTGTAGATAAACATTTCTTTTTGAAGGCAAAAAAATAAAATATTCTGTAAGTCCTCCACAGTTCTGCCAAGACAGCTTAGATCTGGAATGACAAGCGTGTCCTGCCGGTTCATTTCATGATGAATCAATTGTCTGAGTAAACGTTGCTCTGAAACAAC  >  NZ_CP009273/289124‑289362

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