Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451_phaM15,222,9940GT54.5% ‑4.8 / 15.5 33L435M (CTG→ATG) PP_4603 CDS (transposase)
Reads supporting (aligned to +/- strand):  ref base G (0/15);  new base T (0/18);  total (0/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.35e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GAAGGGCTGGGTGATCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTG  >  AE015451_phaM1/5222913‑5223120
                                                                                 |                                                                                                                              
gaggggTTGGGGGATTTCTCATAAAACAAAGTGCAACCCGTTGTTTTCTAGGCACTCAGGTGGCCACTGGTTGATACCCCGTTTTCCGGGCTTGGCGAAGAAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCg                                                                         <  1:8588501/134‑1 (MQ=1)
gaggggCTGGGTGATCTTTCAGCAAACAAACTGCAACCCATGTTTTTCAAATCACTCAAATGTCCACCGGTTGACACCCCGGTTTCTGGGCTTTGTGAAGCAAAGCCCGGACCTGCCGGTCGTGGTTTCTGGTTTCg                                                                         <  1:4049259/134‑1 (MQ=0)
           tggtctctcATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCATTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAGAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAg                                                              <  2:10269649/134‑1 (MQ=0)
           tggtctctcACCAAACAAAATGAAACCCATTTTTTTCTAAGCATTCAAATGGCCACCAGTTGAGACCCCGGTTTCCGGGCTTTGCGAAGCAAAGCCCGTACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAg                                                              <  2:14972047/134‑1 (MQ=1)
                tttcctcaAACAAACTACAACCCCTTTTTTTCTAGCCACTCAGGCGGCCACCGGTTGACCCCCCATTTCCTGGTCTTTACGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCtt                                                         <  1:15703697/132‑1 (MQ=1)
                tctaagcaaaaaaaATACAACCCGTTGTTTTCGAAGCACTCAAGCGGCCACCAGTTGATACCCCGGTTTCCGGGTTTGGCGAAGAAAACCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCtt                                                         <  1:3001336/131‑1 (MQ=1)
                tcccatcaAACAAAATTCAACCCATTGTCTTCTAAGCACTCAAGCGGCCACCATTTGATACCCCATTTTCCGGGCTTTGCGAATCAAATCCCGGACCTGCCGGTCTTTGCTTCTGGTTTCGAACTCCTCAAGACCtt                                                         <  2:6164418/134‑1 (MQ=1)
                tcccatcaAACAAAATGCAACCTATTGTCTTCCAAGCACTCAAGCTGCCACCAGTTGATACCTCAGTTTCCGGGCTTGGCGAGGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCtt                                                         <  2:14374766/134‑1 (MQ=1)
                tcccataaAACAAACTACAACCCATTGTTTTCGAAGCACTCAGGCGGCCACCAGTTGATACCCCGTTTTCTGGGTTTTGCGAAACAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCtt                                                         <  1:1936528/134‑1 (MQ=1)
                tcccataaAAAAAAAGACAACCCGTTGTTTCCGAAACACCCAGGTGGCCACCGGTTGATACCCCGTTTTCGGGGTTTGGAGAAAAAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCtt                                                         <  1:4618515/134‑1 (MQ=0)
                ccccaccaaacaaaGTGCAACCCGTTTTTTTCGAAACACTCAAGCGGCCACAGGTTGACACCCAATTTTCTGGTCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCtt                                                         <  1:6192510/131‑1 (MQ=1)
                               tgaaaCCCGTTGTTTTCGAAGCACTCAAGCGGCCACCAGTGGATCCCCCAGTTTCCGGGCTTGGAGAAAAAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCt                                          <  1:15366452/134‑1 (MQ=1)
                                        tttttttttAAGCACTCAAGCGGCCACCATGTGATACCCCATTTTCCGGGCTTTGCGAAACAAAGCCCGGACCTGCCGGTCGTGGCTTCTGGTTTCGAACTCCTCAAGACTGTGGTCAATATAAGCCTGCTTCTTgg                                 <  2:10670147/134‑1 (MQ=1)
                                        tttttttCTAAGCACTCAAGCGGCCACCGGTTGATCCCCCATTTTCCGGGCTTTGCGAAGCAAACCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTgg                                 <  1:4094990/134‑1 (MQ=1)
                                        tttttttCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCTGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTgg                                 <  1:2932326/134‑1 (MQ=1)
                                        tttttttCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGTGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTgg                                 <  2:5618242/134‑1 (MQ=1)
                                        tttttttCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTgg                                 <  1:11222354/134‑1 (MQ=1)
                                        tttttttCTAAGCACCCAAGCGGCCACCATTTGAGCCCCCGTTTTCCGGGCTTGGAGAAAAAAAACCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTgg                                 <  1:14935535/134‑1 (MQ=1)
                                        tttttttCTAAACACTCAAGCGTCCACCAGTTAATCCCCCATTTTCCGGGCTTTGTGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTgg                                 <  1:9248904/134‑1 (MQ=1)
                                        tttttttCGAAACACCCAAGCGGCCACCAGTTAACACCCCATTTTCTGGGCTTGGAGAAGCAAATCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTgg                                 <  1:4939854/134‑1 (MQ=1)
                                        tttttttCAAAGCACTCAAGAGGCAACCGTTTAACACCCCGTTTTCTGGGTTTGGAGAAAAAAATCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTgg                                 <  1:729991/134‑1 (MQ=1)
                                          gtcttCTAAGCACTCAAGCGGCCATCGGTTGCTACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGtt                               <  2:742422/134‑1 (MQ=0)
                                           tttttgaaGCACTCAAGTGGCCACCGGTTGATACCCCATTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTa                              <  1:3217515/131‑1 (MQ=1)
                                              tcgaaGCACTAAAGCGGCCACCGGTTGATACCCCGTTTTCTGGGTTTTGCGAAGAAAATCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACa                           <  1:7050610/134‑1 (MQ=1)
                                              tcgaaGCACCAAAGCGGCCACCGGTTGATACCCTGTTTTTGGGGTTGTGAGAAGAAAATCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACa                           <  1:3600664/134‑1 (MQ=1)
                                              tcgaaGAACTCAAGCGGCAACCAGTTGATACCCAAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACa                           <  1:8359289/134‑1 (MQ=0)
                                              tcaaaGCACCCAAGCGGCCACCAGTGGATACCCCATTTTCCGGGCTTTGCGAAACAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACa                           <  1:4751451/134‑1 (MQ=1)
                                               tgaaacACTCAAGCGGCCACCGGTGGATACCCAGGTTTCTGGGTTTGGAGAAAAAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAg                          <  1:8297731/132‑1 (MQ=1)
                                               cgaaacACTCAAGCGGCCACCAGTTGATACCCCGGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAg                          <  1:4214526/132‑1 (MQ=0)
                                               cgaaacACTAAAGTGGCCACCGGTGGATACCCCGGTTTCCGGGTTTTGAGAAGAAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAg                          <  1:1268624/132‑1 (MQ=1)
                                                gaatcACTCAAGCGGCCACCATTGGATACCCCATTTTCTGGTCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAgg                         <  1:13613767/133‑1 (MQ=1)
                                                         aagtgaccaccgGTTGATACCCCAGTTTCCGGGTTTTGCAAAGAAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACAcc                <  1:14471072/131‑1 (MQ=1)
                                                                       tggtACCCCATTTTCCGGGCTTTGAGAAAAAAACCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTg  <  1:9268641/134‑1 (MQ=1)
                                                                                 |                                                                                                                              
GAAGGGCTGGGTGATCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTG  >  AE015451_phaM1/5222913‑5223120

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑