Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,429,3450TG52.6% ‑4.3 / 14.7 19R99R (CGT→CGG) stfRputative prophage side tail fiber protein StfR
Reads supporting (aligned to +/- strand):  ref base T (3/6);  new base G (6/4);  total (9/10)
Fisher's exact test for biased strand distribution p-value = 3.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.77e-02
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

AGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGT‑ACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGTGGCACAGAACACGGCAGCCGCGAAGAAGTCAGCCAGTGATGCCAGCACATCAGCCCGTGAGGCGGCAACCCATGCGGCTGATGCTGCGGACTCAGCACGCGCAGCCAGCACGTCAGCCGGACAGGCCGCGTCGTCGGCTCAGTCAGCGTCTTCCAGCGCAGGAACGGCATCAACAAAGG  >  NC_000913/1429101‑1429586
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                 
aGAACTGCACAAACCGGCGGAAAGCAAAAGGTAGCAGCACCACGGGGGTGGTGAACGCGGTGCCCTCAGAAAATCCGGATGAATCCGTCTGTTACAGCATGGACGTTG‑GTGAGGGTAATTATGGCGTCACTGGGTTTGTGGAAGGATTCCCGCCGTCACATGCCGGGGCTATTGCCGTGTATGAAGATTCTCAGCCCGGGACGCTGAATGATTTTCTCGGTCCCAGGACGGAGGATGATGCCCGTCCGGAg                                                                                                                                                                                                                                             <  1:198554/251‑1 (MQ=255)
         caATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGT‑ACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTgcgc                                                                                                                                                                                                                                    >  1:169723/1‑251 (MQ=255)
                                                       cacGTTGGCCTAAGAAAATCCGGATGAAGCCGGGCGTTACAGCAGGGAGGTGGAGT‑ACGGTCAGTACAGCGTTATTCTGTGGGTGGAAGGCTTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACGGATGGTGGAAGAGGTGGTGCGTAACGCGTCCGCgg                                                                                                                                                                                      <  1:172378/251‑1 (MQ=255)
                                                               gcTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGT‑ACGGTCAGTAAAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGCTTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCAATGACGGAGGATGATGCCCGGCCGGAGGTGCTGCGTCGTATTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGt                                                                                                                                                                                            >  1:150114/2‑237 (MQ=255)
                                                               cctaaGAAAATCCGGATGAAGCCGGGCATTCCAGCATGGACGTTGAGT‑ACGGTCGGTACAGCGGTATTCTGTTGGTGGATGGTTTCCCGCCGTCACATGCCGTGACCACTACCGTGTATGAAGAGTCCCAACCCGGTACGCTGAATGATTGTCTCGGTGACATGACGGAGGATGATGCCCGTCCGGAGGCCCTGCGCCGTTTTGAACTGATGGTGGACGAGGTGGGGCTTAACGCGTGCGCGGTGGCACAg                                                                                                                                                                              <  1:16949/247‑1 (MQ=255)
                                                                      gaaTCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGT‑ACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGAGGTGCTGCGTCGTCTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGGGGTGGCACAGagtacgg                                                                                                                                                                       >  1:93871/2‑245 (MQ=255)
                                                                            ggATGAAGCCGGGCGTTACAGCATGGATGTGGAGT‑ACGGCCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCCCCATCGCACGCCGGGACCATCACCGTGTATGCAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGAGGTGCTGCGTCGTCTTGAACTGATGGTGGCAGAGGTGGAGCGTAACGCGTCCGTGGTGGCACAGAGTACGGCAgacg                                                                                                                                                                 >  1:29840/1‑248 (MQ=255)
                                                                                      gggCGTTACAGCATGGACGTTGAGT‑ACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCAGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGTg                                                                                                                                                                                    <  1:17771/222‑1 (MQ=255)
                                                                                          gTTACAGCATGGCTGTGGAGT‑ACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGCTGACGCCCGGCCGGAGGTGCTGCGTCGTCTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGTGGTGGCACAGAGTACGGCAGACGCGAAGAAATCAGcc                                                                                                                                                   >  1:218268/1‑251 (MQ=255)
                                                                                          gTTACAGCATGGATGTGGAGT‑ACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGATGTGCTGCGTCGTCGTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGTGGTGGCACAGAGTACGGCAGACGCGAAGAAATCAGcc                                                                                                                                                   <  1:73144/251‑1 (MQ=255)
                                                                                          gTTACAGCATGGACGTTGAGT‑ACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGTGGCACAGAACACGGCAGCCGCGAAGAAGTCAGcc                                                                                                                                                   >  1:99726/1‑251 (MQ=255)
                                                                                                                                                        gccgTTACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGTGACACAGAACACGGCAGcc                                                                                                                                                                  >  1:16872/1‑175 (MQ=255)
                                                                                                                                                         ccatcgcacgCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGAGGTGCTGCGTCGTCTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGTGGTGGCCCAGa                                                                                                                                                                             <  1:4876/154‑1 (MQ=255)
                                                                                                                                                               caTGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGGTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGTGGCACAGAACACGGCAGCCGCGAAGAAGTCAGCCAGTGATGCCAGCACATCAGCCCGTGAGGCGGCAACCCATGCGGCTGATGCTGCGGACTCAGCACgcgc                                                                               <  1:148192/251‑1 (MQ=255)
                                                                                                                                                               caTGCCGGGACCATGACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGGCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCTCGTAACGCGTCCGCGGTGGCACAGAACACGGCAGCCGCGAAGAAGTCAGCCAGTGATGCCAGCACATCAGCCCGTGAGGCGGCAACCCATGCGGCTGATGCTGCGGCCTCAGCACgcgc                                                                               <  1:148199/251‑1 (MQ=255)
                                                                                                                                                               caTGCCGGGAACATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTGTCTCGGTTCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCACCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGTGGCACAGAACACGGCAGCCGCGAAGAAGTCAGCCATTGATGCCCGCACATCAGCCCGTGAGGCTGCAACCCATGCGGCTGATGCTGCGGACTCAGCACgcgc                                                                               <  1:148193/251‑1 (MQ=255)
                                                                                                                                                                       gACCATCACCGTGTATGCAGATTCACAACTGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGAGGTGCTGCGTCGTCTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCgtggtg                                                                                                                                                                                    >  1:66169/1‑142 (MQ=255)
                                                                                                                                                                                                                                            tgatgCCCGGCCGGAGGTGCTGCGTCGTCTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGAGGTGGCACAGAGTACGGCAGACGCGAAGAAATCAGCCGGCGATGCCAGTGCATCAGCTGCTCAGGTCGCGGCCCTTGTGACTGATGCAACTGACTCAGCACGCGCCGCCAGCACGTCCGCCGGACAGGCTGCATCGTCAGCTCAGGAAGCGTCCTCCGGCGCAGAAGCGGCATCAGCAAAgg  >  1:97650/1‑251 (MQ=255)
                                                                                                                                                                                                                                                gCCCGGCCGGAGGTGCTGCGTCGTCTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGTGGTGGCACAGAGTACGGCAGACGCGAAGAAATCAGCCGGCGATGCCAGTGCATCAGCTGCTCAGGTCGCGGCCCTTGTGACTGATGCAACTGACTCAGCACGCGCCGCCAGCACGTCCGCCGGACAgg                                                        <  1:197851/193‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                 
AGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGT‑ACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGTGGCACAGAACACGGCAGCCGCGAAGAAGTCAGCCAGTGATGCCAGCACATCAGCCCGTGAGGCGGCAACCCATGCGGCTGATGCTGCGGACTCAGCACGCGCAGCCAGCACGTCAGCCGGACAGGCCGCGTCGTCGGCTCAGTCAGCGTCTTCCAGCGCAGGAACGGCATCAACAAAGG  >  NC_000913/1429101‑1429586

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑