Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913331,8790GA58.6% 1.7 / 35.6 29intergenic (+3/‑492)yahV/pdeLprotein YahV/DNA‑binding transcriptional regulator/c‑di‑GMP phosphodiesterase PdeL
Reads supporting (aligned to +/- strand):  ref base G (3/9);  new base A (8/9);  total (11/18)
Fisher's exact test for biased strand distribution p-value = 2.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.45e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

ACCTACATTTGGGCTGTTGCGAATGTTCAATAAGTTTAGTCTTATTTAATGTAAATATTGCTGATCATTTGAAATGACGCATTATTCATGAGAAATGTGTATCGTAAATCAACTGAAATTAACGCAACCATTTGTTATTTAAGGTTTAATTATCTGTGTGTGATATTTTATTGAATGTTTTAAATATTGTTTTTATTGGCATTGCTATAATATTGGTTATCATTTGCTGAATGGATTCAGTCTTAATGAGTGGGTTTTTAAGGGACAGGCATAGAGTAATGATACGTATGCATAACCAACATCTTTACTCATTATGTCATTGAATGTTGACGCTATGTGTTTATGAGGGAGAGGTATTTTCAGTTGATCTGGATTGTTAAATTCATATAATGCGCCTTTGCTCATGAATGGATGCCAGTATGTAGTGGGAAATTATAAATATTGAAATAGTCCAACTACTTCTTTATTACCAAAAATGA  >  NC_000913/331647‑332125
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      
aCCTACATTTGGGCTGTTGCGAATGTTCAATAAGTTTAGTCTTATTTAATGTAAATATTGCTGATCATTTGAAATGACGCCTTATTCATGAGAAATGTGTATCGTAAATCAACTGAAATTAACGCAACCATTTGTTATTTAAGGTTTAATTATCTGTGTGTGATATTTTATTGAAGGTTTTAAATATTGTTTTTATTGGCATTGCTATAATATTGGTTATCATTTGCTGAATGGATTCAGTCTTAATTAGt                                                                                                                                                                                                                                      >  2:83584/1‑251 (MQ=255)
     cATTTGGGCTGTTGCGAATGTTCAATAAGTTTAGTCTTATTTAATGTAAATATTGCTGATCATTTGAAATGACGCATTATTCATGAGAAATGTGTATCGTAAATCAACTGAAATTAACGCAACCATTTGTTATTTAAGGTTTAATTATCTGTGTGTGATATTTTATTGAATGTTTTAAATATTGTTTTTATTGGCATTGCTATAATATTGGTTATCATTTGCTGAATGGATTCAGTCTTAATGAGTGGGtt                                                                                                                                                                                                                                 <  1:241108/251‑1 (MQ=255)
         tGGGCTGTTGCGAATGTTCAATAAGTTTAGTCTTATTTAATGTAAATATTGCTGATCATTTGAAATGACGCATTATTCATGAGAAATGTGTATCGTAAATCAACTGAAATTAACGCAACCATTTGTTATTTAAGGTTTAATTATCTGTGTGTGATATTTTATTGAATGTTTTAAATATTGTTTTTATTGGCATTGCTATAATATTGGTTATCATTTGCTGAATAGATTCAGTCTTAATGAGTGGGttttt                                                                                                                                                                                                                              >  1:2615/1‑250 (MQ=255)
                        gTTCAATAAGTTTAGTCTTATTTAATGTAAATATTGCTGATCATTTGAAATGACGCATTATTCATGAGAAATGTGTATCGTAAATCAACTGAAATTAACGCAACCATTTGTTATTTAAGGTTTAATTATCTGTGTGTGATATTTTATTGAATGTTTTAAATATTGTTTTTATTGGCATTGCTATAATATTGGTTATCATTTGCTGAATGGATTCAGTCGTAATGAGTGGGTTTTTAAGGGACAGGCATaga                                                                                                                                                                                                              <  1:115230/251‑1 (MQ=255)
                               aaGTTTAGTCTTATTTAATGTAAATATTGCTGATCATTTGAAATGACGCATTATTCATGAGAAATGTGTATCGTAAATCAACTGAAATTAACGCAACCATTTGTTATTTAAGGTTTAATTATCTGTGTGTGATATTTTATTGAATGTTTTAAATATTGTTTTTATTGGCATTGCTATAATATTGGTTATCATTTGCTGAATGGATTCAGTCTTAATGAGTGGGTTTTTAAGGGACAGGCATAGAGTAATGa                                                                                                                                                                                                       <  2:139407/251‑1 (MQ=255)
                                       tCTTATTTAATGTAAATATTGCTGATCATTTGAAATGACGCATTATTCATGAGAAATGTGTATCGTAAATCAACTGAAATTAACGCAACCATTTGTTATTTAAGGTTTAATTATCTGTGTGTGATATTTTATTGAATGTTTTAAATATTGTTTTTATTGGCATTGCTATAATATTGGTTATCATTTGCTGAATAGATTCAGTCTTAATGAGTGGGTTTTTAAGGGACAGGCATAGAGTAATGATACGTATg                                                                                                                                                                                               >  1:230673/1‑251 (MQ=255)
                                            tttAATGTAAATATTGCTGATCATTTGAAATGACGCATTATTCATGAGAAATGTGTATCGTAAATCAACTGAAATTAACGCAACCATTTGTTATTTAAGGTTTAATTATCTGTGTGTGATATTTTATTGAATGTTTTAAATATTGTTTTTATTGGCATTGCTATAATATTGGTTATCATTTGCTGAATAGATTCAGTCTTAATGAGTGGGTTTTTAAGGGACAGGCATAGAGTAATGATACGTATGCATa                                                                                                                                                                                           >  2:172946/1‑250 (MQ=255)
                                              tAATGTAAATATTGCTGATCATTTGAAATGACGCATTATTCATGAGAAATGTGTATCGTAAATCAACTGAAATTAACGCAACCATTTGTTATTTAAGGTTTAATTATCTGTGTGTGATATTTTATTGAATGTTTTAAATATTGTTTTTATTGGCATTGCTATAATATTGGTTATCATTTGCTGAATGGATTCAGTCTTAATGAGTGGGTTTTTAAGGGACAGGCATAGAGTAATGATACGTATGCATAAc                                                                                                                                                                                         <  1:229918/250‑1 (MQ=255)
                                                    aaaTATTGCTGATCATTTGAAATGACGCATTATTCATGAGAAATGTGTATCGTAAATCAACTGAAATTAACGCAACCATTTGTTATTTAAGGTTTAATTATCTGTGTGTGATATTTTATTGAATGTTTTAAATATTGTTTTTATTGGCATTGCTATAATATTGGTTATCATTTGCTGAATGGATTCAGTCTTAATGAGTGGGTTTTTAAGGGACAGGCATAGAGTAATGATACGTATGCATAACCAACAt                                                                                                                                                                                   <  1:157747/250‑1 (MQ=255)
                                                     aaTATTGCTGATCATTTGAAATGACGCATTATTCATGAGAAATGTGTATCGTAAATCAACTGAAATTAACGCAACCATTTGTTATTTAAGGTTTAATTATCTGTGTGTGATATTTTATTGAATGTTTTAAATATTGTTTTTATTGGCATTGCTATAATATTGGTTATCATTTGCTGAATAGATTCAGTCTTAATGAGTGGGTTTTTAAGGGACAGGCATAGAGTAATGATACGTATGCATAACCAACATCt                                                                                                                                                                                 >  2:191050/1‑251 (MQ=255)
                                                             tGATCATTTGAACTGACGCATTATTCATGAGAAATGTGTATCGTAAATCAACTGAAATTAACGCAACCATTTGTTATTTAAGGTTTAATTATCTGTGTGTGATATTTTATTGAATGTTTTAAATATTGTTTTTATTGGCATTGCTATAATATTGGTTATCATTTGCTGAATAGATTCAGTCTTAATGAGTGGGTTTTTAAGGGACAGGCATAGAGTAATGATACGTATGCATAACCAACATCTTTACTCAt                                                                                                                                                                         >  2:112995/1‑251 (MQ=255)
                                                             tGATCATTTGAAATGACGCATTATTCATGAGAAATGTGTATCGTAAATCAACTGAAATTAACGCAACCATTTGTTATTTAAGGTTTAATTATCTGTGTGTGATATTTTATTGAATGTTTTAAATATTGTTTTTATTGGCATTGCTATAATATTGGTTATCATTTGCTGAATAGATTCAGTCTTAATGAGTGGGTTTTTAAGGGACAGGCATAGAGTAATGATACGTATGCATAACCAACATCTTTACTCAt                                                                                                                                                                         >  2:65400/1‑251 (MQ=255)
                                                                        aaTGACGCATTATTCATGAGAAATGTGTATCGTAAATCAACTGAAATTAACGCAACCATTTGTTATTTAAGGTTTAATTATCTGTGTGTGATATTTTATTGAATGTTTTAAATATTGTTTTTATTGGCATTGCTATAATATTGGTTATCATTTGCTGAATGGATTCAGTCTTAATGAGTGGGTTTTTAAGGGACAGGCATAGAGTAATGATTCGTATGCATAACCAACATCTTTACTCATTATGTCATTGa                                                                                                                                                              <  1:305674/251‑1 (MQ=255)
                                                                             cGCATTATTCATGAGAAATGTGTATCGTAAATCAACTGAAATTATCGCCACCATTTGTTATTTAATGTTTAATTATCTGTTTGTGATATTTTATTGAATGTTTTAAATATTGTTTTTATTGGCATTGCTATAATATTGGTTATCATTTGCTTAATAGATTCAGTCTTAATGAGTGGGTTTTTAAGGGACAGGCATATAGTAATGATACGTATGCATATCCAACATCTTTACTCATTATGTCATTGAATGtt                                                                                                                                                         <  1:294961/251‑1 (MQ=255)
                                                                             cGCATTATTCATGAGAAATGTGTATCGTAAATCAACTGAAATTAACGCAACCATTTGTTATTTAAGGTTTAATTATCTGTGTGTGATATTTTATTGAATGTTTTAAATATTGTTTTTATTGGTATTGCTATAATATTGGTTATCATTTGCTGAATGGATTCAGTCTTAATGAGTGGGTTTTTAAGGGACAGGCATAGAGTAATGATACGTATGCATAACCAACATCTTTACTCATTATGTCATTGAATGtt                                                                                                                                                         >  1:60762/1‑251 (MQ=255)
                                                                               cATTATTCATGAGAAATGTGTATCGTAAATCAACTGAAATTAACGCAACCATTTGTTATTTAAGGTTTAATTATCTGTGTGTGATATTTTATTGAATGTTTTAAATATTGTTTTTATTGGCATTGCTATAATATTGGTTATCATTTGCTGAATGGATTCAGTCTTAATGAGTGGGTTTTTAAGGGACAGGCATAGAGTAATGATACGTATGCATAACCAACATCTTTACTCATTATGTCATTGAATGTTGa                                                                                                                                                       <  1:39206/251‑1 (MQ=255)
                                                                                                      cGTAAATCAACTGAAATTAACGCAACCATTTGTTATTTAAGGTTTAATTATCTGTGTGTGATATTTTATTGAATGTTTTAAATATTGTTTTTATTGGCATTGCTATAATATTGGTTATCATTTGCTGAATGGATTCAGTCTTAATGAGTGGGTTTTTAAGGGACAGGCATAGAGTAATGATACGTATGCATAACCAACATCTTTACTCATTATGTCATTGAATGTTGACGCTAtgtgctgatgaatcccct                                                                                                                                <  1:197141/251‑15 (MQ=255)
                                                                                                               aCTGAAATTAACGCAACCATTTGTTATTTAAGGTTTAATTATCTGTGTGTGATATTTTATTGAATGTTTTAAATATTGTTTTTATTGGCATTGCTATAATATTGGTTATCATTTGCTGAATGGATTCAGTCTTAATGAGTGGGTTTTTAAGGGACAGGCATAGAGTAATGATACGTATGCATAACCAACATCTTTACTCATTATGTCATTGAATGTTGACGCTAtgtgctgatgaatcccctaatgatttt                                                                                                                       >  1:275777/1‑228 (MQ=255)
                                                                                                                             aaCCATTTGTTATTTAAGGTTTAATTATCTGTGTGTGATATTTTATTGAATGTTTTAAATATTGTTTTTATTGGCATTGCTATAATATTGGTTATCATTTGCTGAATAGATTCAGTCGTAATGAGTGGGTTTTTAAGGGACAGGCATAGAGTAATGATACGTATGCATAACCAACATCTTTACTCATTATGTCATTTAATGTTGTCGCTATGTGTTTATGAGGGAGAGGTATTTTCAGTTGATCTGGATTg                                                                                                         <  1:43004/251‑1 (MQ=255)
                                                                                                                                                           gtgtgtgATATTTTATTGAATGTTTTAAATATTGTTTTTATTGGCATTGCTATAATATTGGTTATCATTTGCTGAATAGATTCAGTCTTAATGAGTGGGTTTTTAAGGGACAGGCATAGAGTAATGATACGTATGCATAACCAACATCTTTACTCATTATGTCATTGAATGTTGACGCTATGTGTTTATGAGGGAGAGGTATTTTCAGTTGATCTGGATTGTTAAATTCATATAATGCGCCTTTGCTc                                                                              <  1:24935/248‑1 (MQ=255)
                                                                                                                                                                                                   ttGGCATTGCTATAATATTGGTTATCATTTGCTGAATAGATTCAGTCTTAATGAGTGGGTTTTTAAGGGACAGGCATAGAGTAATGATACGTATGCATAACCAACATCTTTACTCATTATGTCATTGAATGTTGACGCTATGTGTTTATGAGGGAGAGGTATTTTCAGTTGATCTGGATTGTTAAATTCata                                                                                              <  1:234030/192‑1 (MQ=255)
                                                                                                                                                                                                   ttGGCATTGCTATAATATTGGTTATCATTTGCTGAATAGATTCAGTCTTAATGAGTGGGTTTTTAAGGGACAGGCATAGAGTAATGATACGTATGCATAACCAACATCTTTACTCATTATGTCATTGAATGTTGACGCTATGTGTTTATGAGGGAGAGGTATTTTCAGTTGATCTGGATTGTTAAATTCata                                                                                              >  2:234030/1‑192 (MQ=255)
                                                                                                                                                                                                                      ggTTATCATTTGCTGAATAGATTCAGTCTTAATGAGTGGGTTTTTAAGGGACAGGCATAGAGTAATGATACGTATGCATAACCAACATCTTTACTCATTATGTCATTGAATGTTGACGCTATGTGTTTATGAGGGAGAGGTATTTTCAGTTGATCTGGATTGTTAAATTCATATAATGCGCCTTTGCTCATTAATGGATGCCAGTATGTAGTGGGAAATTATAAATATTGAAATAGTCCAACTACTTCttt                <  1:177701/251‑1 (MQ=255)
                                                                                                                                                                                                                          aTCATTTGCTGAATAGATTCAGTCTTAATGAGTGGGTTTTTAAGGGACAGGCATAGAGTAATGATACGTATGCATAACCAACATCTTTACTCATTATGTCATTGAATGTTGACGCTATGTGTTTATGAGGGAGAGGTATTTTCAGTTGATCTGGATTGTTAAATTCATATAATGCGCCTTTGCTCATGAATGGATGCCAGTATGTAGTGGGAAATTATAAATATTGAAATAGTCCAact                        <  2:148955/239‑1 (MQ=255)
                                                                                                                                                                                                                          aTCATTTGCTGAATAGATTCAGTCTTAATGAGTGGGTTTTTAAGGGACAGGCATAGAGTAATGATACGTATGCATAACCAACATCTTTACTCATTATGTCATTGAATGTTGACGCTATGTGTTTATGAGGGAGAGGTATTTTCAGTTGATCTGGATTGTTAAATTCATATAATGCGCCTTTGCTCATGAATGGATGCCAGTATGTAGTGGGAAATTATAAATATTGAAATAGTCCAact                        >  1:148955/1‑239 (MQ=255)
                                                                                                                                                                                                                          aTCATTTGCTGAATAGATTCAGTCTTAATGAGTGGGTTTTTAAGGGACAGGCATAGAGTAATGATACGTATGCATAACCAACATCTTTACTCATTATGTCATTGAATGTTGACGCTATGTGTTTATGAGGGAGAGGTATTTTCAGTTGATCTGGATTGTTAAATTCATATAATGCGCCTTTGCTCATGAATGGATGCCAGTATGTAGTGGGAAATTATAAATATTGAAATAGTCCAACTACTTCTttatt             <  1:292191/250‑1 (MQ=255)
                                                                                                                                                                                                                                tGCTGAATGTATTCAGTCTTAATGAGTGGGTTTTTAAGGGACAGGCATAGAGTAATGATACGTATGCATAACCAACATCTTGACTCATTATGTCATTGAATGTTGACGCTATGTGTTGATGAGGGAGAGGTATTGTCAGTTGATCTGGATTGTTAAAGTCATATAATGCGCCTTTGGTCATGAATGGATGGCAGTATGTATTTGGAAATTATAAATATTGAAATAGTCCAACTACTTCTTTATTACCaaa       <  1:178955/250‑1 (MQ=255)
                                                                                                                                                                                                                                    gAATAGATTCAGTCTTAATGAGTGGGTTTTTAAGGGACAGGCATAGAGTAATGATACGTATGCATAACCAACATCTTTACTCATTATGTCATTGAATGTTGACGCTATGTGTTTATGAGGGAGAGGTATTTTCAGTTGATCTGGATTGTTAAATTCATATAATGCTCCTTTGCTCATGAATGGATGCCAGTATGTAGTGGGAAATTATAAATATTGAAATAGTCCAACTACTTCTTTATTACCAAAAAt    <  1:145408/249‑1 (MQ=255)
                                                                                                                                                                                                                                    gAATAGATTCAGTCTTAATGAGTGGGTTTTTAAGGGACAGGCATAGAGTAATGATACGTATGCATAACCAACATCTTTACTCATTATGTCATTGAATGTTGACGCTATGTGTTTATGAGGGAGAGGTATTTTCAGTTGATCTGGATTGTTAAATTCATATAATGCGCCTTTGCTCATGAATGGATGCCAGTATGTAGTGGGAAATTATAAATATTGAAATAGTCCAACTACTTCTTTATTACCAAAAATGa  <  2:177983/251‑1 (MQ=255)
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      
ACCTACATTTGGGCTGTTGCGAATGTTCAATAAGTTTAGTCTTATTTAATGTAAATATTGCTGATCATTTGAAATGACGCATTATTCATGAGAAATGTGTATCGTAAATCAACTGAAATTAACGCAACCATTTGTTATTTAAGGTTTAATTATCTGTGTGTGATATTTTATTGAATGTTTTAAATATTGTTTTTATTGGCATTGCTATAATATTGGTTATCATTTGCTGAATGGATTCAGTCTTAATGAGTGGGTTTTTAAGGGACAGGCATAGAGTAATGATACGTATGCATAACCAACATCTTTACTCATTATGTCATTGAATGTTGACGCTATGTGTTTATGAGGGAGAGGTATTTTCAGTTGATCTGGATTGTTAAATTCATATAATGCGCCTTTGCTCATGAATGGATGCCAGTATGTAGTGGGAAATTATAAATATTGAAATAGTCCAACTACTTCTTTATTACCAAAAATGA  >  NC_000913/331647‑332125

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