Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451_phaM15,223,0130CA42.9% ‑4.9 / 24.1 42L428F (TTG→TTT) PP_4603 CDS (transposase)
Reads supporting (aligned to +/- strand):  ref base C (1/23);  new base A (0/18);  total (1/41)
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.

GGAAGGGCTGGGTGATCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCG‑AAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGATGCGCAGTGGCCTTGATGG  >  AE015451_phaM1/5222912‑5223140
                                                                                                      |                                                                                                                               
gggaGGTCTGGTTGATCTCCCCTCAAACAAAATACAACCCATTGTTTTCCAAACATTAAATCGGCCACCAGTTGATACCCCGTTTTCCGGGTTTTGCG‑AAGCAAAACCCGGACCTGCCGGTCTTGGCTTCTGGTTTc                                                                                              <  1:18539453/134‑1 (MQ=1)
gggaGGGTTGGGTGATCTCTCATAAAAAAAAGTGCACCCCGTTTTTTTCGAAGCACTCAAGCGACCACCAGTGGATACCCCATTTTCCGGGCTTTGCG‑AAGCAAATCCCGGACCTGCCGGTCTTGGCTTCTGGTTTc                                                                                              <  1:11441581/134‑1 (MQ=1)
gggaGGGCGGGGTGACCTCCCACCAAAAAAAATGCAACCCGTTGTTTCCTAGACACTCAAGCGGCCACCATTTAATACCCCGTTTTCTGTTCTTTGCG‑AAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTc                                                                                              <  2:9978169/134‑1 (MQ=1)
ggcaGGGCTGGGTGATCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAACCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCG‑AAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTc                                                                                              <  2:14459199/134‑1 (MQ=1)
 gaggggTTGGGGGATCTCTCACAAAAAAAAAGGCAACCCATTGTTTTCGAAGCACTAAAGCGGCCACAATTTGATACCCCAGTTTCTGGGTTTGGAG‑AAGAAAATCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCg                                                                                             <  1:9167201/134‑1 (MQ=1)
            tggtctctcATAAAACAAAATGAAACCCGTTGTTTTCTAAGCACTAAAGCGGCCACCATTTGATACCCCAGTTTCCGGGTTTTGCG‑AAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAg                                                                                  <  1:3074429/134‑1 (MQ=1)
            tgctctctcGTAAAAAAAAATGCAACCCGTTTTTTTCTAAGCACTCAGGCGGCCACCGGTGGATACCCCGTTTTCTGGGTTTGGTG‑AAACAAATCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAg                                                                                  <  1:2946447/134‑1 (MQ=1)
                 tcgcatcaAAAAAAATGCAACCCATTGTTTCCTAAACACTAAAGCGGCCACCAGTTGAGACCCCGTTTTCGGGGTGTTGCG‑AGAAAAAGCCCGGACCTGCCGGTCGTGGCTTCTGGTTTCGAACTCCTCAAGACCtt                                                                             <  2:14945937/134‑1 (MQ=1)
                 tcccatcaAAAAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCATTTTCCGGGCTTTGCG‑AAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCtt                                                                             <  2:10855552/134‑1 (MQ=0)
                 tcccataaAAAAAAAGGAAACCCGTTGTTTTCGAAGCACCCAAGGTGCCACCAGTTGAGACCCCAGTTTCCGGGCTTTGCA‑AAGAAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCtt                                                                             <  1:18422540/134‑1 (MQ=1)
                     ataaaaaaaaaTGCAACCCGTTTTTTTCGAAGCACTCAAGCGGCCACCAGTGGATACCCCAGTTTCCGGGCTTTGCG‑AAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTc                                                                         <  1:3850853/134‑1 (MQ=1)
                                  caccccATTGTTTTCTAAGCACTCAAACGGCCACCAGTTGATACCCCATTTTCCGGGCTTTGTG‑AAGCAAAGCCCGGACCGGCCGTTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGc                                                            <  2:1937576/134‑1 (MQ=1)
                                   aacctgttGTTTTCGAAACACTAAAGCGGCCACAGGTGAACACCCTGTTTTTTGGGTTTGGAG‑AAACAAACCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGct                                                           <  1:9889480/131‑1 (MQ=0)
                                         tttttttCTAAGCACTCAAGCGGCCACCGGTTGATACCCCGTTTTCTGGGCTTTGCG‑AAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTgg                                                     <  2:985133/134‑1 (MQ=1)
                                         tttttttCTAAGCACTCAAGCGGCCACCGGTTGATACCCCGTTTTCCGGGCTTTGCG‑AAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTgg                                                     <  1:3065132/134‑1 (MQ=0)
                                         tttttttCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCG‑AAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTgg                                                     <  2:15694161/134‑1 (MQ=1)
                                         tttttttCGAAGCACTAAAGAGACCACCGGTTGACACCCCATTTTCCGGGCTTTGCG‑AAAAAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTgg                                                     <  1:9556001/134‑1 (MQ=1)
                                         tttttttCGAAGCACCCAATCGACCACCAGTTGATACCCCGTTTTCTGGGCTTTGCG‑AAAAAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTgg                                                     <  1:6093066/134‑1 (MQ=1)
                                         tttttttCGAAGCACCCAAGTGGCCACCGGTTGATACCCCAGTTTCTGGGCTTTGCG‑AAAAAAACCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTgg                                                     <  1:5039206/134‑1 (MQ=1)
                                         tttttttCGAAGAACCAAAGCGGCCACCGTTTAATCCCCCAGTTTCTGGGCTTGGAG‑AAGAAAAACCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTgg                                                     <  1:5965079/134‑1 (MQ=1)
                                         tttttttCCAAGCACTCAAATGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCA‑AAGAAAATCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTgg                                                     <  1:14364190/134‑1 (MQ=1)
                                         ctgctttCTAAACACTCAAGCGGCCACCAGTTGATACCCCATTTTCCGGGCGTTGCG‑AAAAAAAGCCCGGACCTGCCGGTCTTGACTCCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCCTCTTgg                                                     <  2:10958760/133‑1 (MQ=1)
                                           gtttgatAAGCACTAAAGAGGCAAAAAGTTGCCACCCCAGTTTCCGGGTTGTGCG‑AGGCAAGGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGtt                                                   <  1:3508037/131‑1 (MQ=1)
                                              ttttAAGCACCAAAGTGGCCACCGGTTGATACCCCGTTTTCTGGGTTTTGCG‑AAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAAc                                                <  1:5193038/134‑1 (MQ=1)
                                               ttgaaacACTCAACCGGCCACAAGTTGATACCCCGGTTTCTGGGCTTTGCG‑AAGCAAGGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACa                                               <  1:6148739/131‑1 (MQ=1)
                                               tcgaaGCAATCAAGGGGCCATTGTTGGATACCCTGGTTT‑CGGGGTTTGAGAAAAAAAACCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACa                                               <  1:134525/134‑1 (MQ=0)
                                                 gaaacACTCAAGCGGCCACCAGTTAATACCCCAGTTTCCGGGCTTTGAG‑AAGAAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAgg                                             <  1:7081521/133‑1 (MQ=0)
                                                 gaaacACCCAAGCGGCCACCGGTGGATACCCCATTTTCTGGGTTTGGAG‑AAGAAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAgg                                             <  1:3748212/133‑1 (MQ=1)
                                                 caaccACTCAACCGGCCACCATTTGATACCCCATTTTCCGGGCTTTGCG‑AAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAgg                                             <  1:990888/133‑1 (MQ=1)
                                                  aaacACTCAAGCGGCCACCAGTTAATCCCCCAGTTTCCGAGCCTTGCG‑AAGCAAAGCCCGGACCTGCAGGTCTTGGCTTCTGATCTCGAACCCCTCAAGACCTCGTTCAACATCAGCATGCTTCTTAGTTAACAGGt                                            >  2:1219729/4‑137 (MQ=1)
                                                      actaaggCGGCCACCGGTGGATACCCCATTTTCTGGGTGTGGAG‑AAGAAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTa                                        <  1:5649472/131‑1 (MQ=1)
                                                      acccAAGCGACCACCAGTTAACACCCCATTTTCCGGGCTTTGCG‑AAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTa                                        <  1:942799/134‑1 (MQ=1)
                                                          aatcGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCG‑AAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACAcc                                    <  2:3166111/134‑1 (MQ=1)
                                                                ccacaggTTGATACCCCGGTTTCTGGTTTTGGAG‑AAAAAAACCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGt                              <  1:4786852/131‑1 (MQ=1)
                                                                ccaacggTTGATACCCCGTTTTCTGGTCTTGGAG‑AAAAAAACCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGt                              <  1:4786889/131‑1 (MQ=1)
                                                                       ttaaTACCCCATTTTCTGGGCTGGGAA‑AAAAAAATCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACt                       <  1:10196117/134‑1 (MQ=1)
                                                                            acccaggTTTCCGGGCTTTGCG‑AAGAAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGATgc                  <  1:592760/131‑1 (MQ=0)
                                                                                           ctgtGCG‑AAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGGCGTGCCAACTGATGCGCAGTGGCCTTGATg   <  1:11930253/134‑1 (MQ=1)
                                                                                           ctctGCG‑AAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGATGCGCAGTGGCCTTGATg   <  1:6172155/134‑1 (MQ=1)
                                                                                            ttggCG‑AAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGATGCGCAGTGGCCTTGATgg  <  2:12689604/134‑1 (MQ=1)
                                                                                            ttggCG‑AAGAAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGATGCGCAGTGGCCTTGATgg  <  1:6001308/134‑1 (MQ=1)
                                                                                            ttggCG‑AAACAAATCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGATGCGCAGTGGCCTTGATgg  <  1:5915326/134‑1 (MQ=0)
                                                                                                      |                                                                                                                               
GGAAGGGCTGGGTGATCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCG‑AAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGATGCGCAGTGGCCTTGATGG  >  AE015451_phaM1/5222912‑5223140

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