Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451_phaM15,223,0130CA42.1% ‑5.3 / 20.1 38L428F (TTG→TTT) PP_4603 CDS (transposase)
Reads supporting (aligned to +/- strand):  ref base C (0/22);  new base A (0/16);  total (0/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
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.

GAAGGGCTGGGTGATCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGATGCGCAGTGGCCTTGATGGCGC  >  AE015451_phaM1/5222913‑5223143
                                                                                                    |                                                                                                                                  
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)
                                                                                          tttggCGAAGAAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGATGCGCAGTGGCCTTGATg      <  1:5419622/133‑1 (MQ=0)
                                                                                          tttggCGAAAAAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGATGCGCAGTGGCCTTGATg      <  1:1811265/133‑1 (MQ=0)
                                                                                           ttggCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTGCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTANACCAGACGTGCCAACTGATGCGCAGTGGCCTTGATgg     <  1:17006710/134‑1 (MQ=0)
                                                                                           ttggCGAAAAAAACCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGATGCGCAGTGGCCTTGATgg     <  1:4662036/134‑1 (MQ=0)
                                                                                              gcaaaGAAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGATGCGCAGTGGCCTTGATGgcgc  <  1:7344941/134‑1 (MQ=1)
                                                                                                    |                                                                                                                                  
GAAGGGCTGGGTGATCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGATGCGCAGTGGCCTTGATGGCGC  >  AE015451_phaM1/5222913‑5223143

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: ‑