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

TGATCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGATGC  >  AE015451_phaM1/5222924‑5223124
                                                                      |                                                                                                                                  
tggtctctcATCAAACAAAGTGCAACCCATTTTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTCCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAg                                                                  <  2:16838656/134‑1 (MQ=0)
tggtctctcATAAAAAAAACTGCACCCTGTTGTTTTCGAAGCACCCAATCGGCCACCGTTGGATACCCCTTTTTCTGGGCTTTGCGAAAAAAACCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAg                                                                  <  1:15002088/134‑1 (MQ=0)
tggtctctcATAAAAAAAAATGAAACCCGTTGTTTTCTAAGAACCCAGGCGGCCACCAGTTGATACCCCAGTTTCTGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAg                                                                  <  1:3554150/134‑1 (MQ=1)
     tcccatcaAACAAAATGCAACCCGTTGTTTTCGAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCtt                                                             <  2:11787830/134‑1 (MQ=0)
     tcccatcaAACAAAATGCAACCCATTTTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCtt                                                             <  2:9821634/134‑1 (MQ=1)
     tcccatcaAACAAAATGCAACCCATTGTTTTCGAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCtt                                                             <  1:4786867/134‑1 (MQ=1)
     tcccatcaAACAAAATGCAACCCATTGTTTTCGAAGCACCAAAGCGGCCACCAGTTGATACCCCGGTTTCCGGGCTTTGCGAAGAAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCtt                                                             <  1:19456466/134‑1 (MQ=1)
     tcccatcaAAAAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCGTTTGATACCCCGGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCtt                                                             <  1:11794486/134‑1 (MQ=1)
     tcccatcaAAAAAAATGAAACCCGTTGTTTTCGAAGCACCCAAGCGGCCACCAGTGGAGACCCCATTTTCCGGGCTTGGCAAAGAAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCtt                                                             <  1:12055438/134‑1 (MQ=1)
     tcccatcaAAAAAAATAAAACCCGTTGTTTTCGAAGAACTCAAGCGGCCATCAGTGGATACCCCGGTTTCTGGTCTTTGCGAAACAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCtt                                                             <  1:17433781/134‑1 (MQ=1)
     tcccataaAAAAAAATGAAACCCATTGTTTTTGAAACACTCAGGCGGCCACCAGTTAATACCCCGGTTCCCGGGTTTTGAGAAAAAAACCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCtt                                                             <  1:13421097/134‑1 (MQ=1)
     tcccataaAAAAAAAGACAACCCGTTGTTTTCGAAACACCCAAGTGGCCACCGTTGGATACCCCATTTTCCGGGTTTGGCGAAGAAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCtt                                                             <  1:6165164/134‑1 (MQ=1)
     tcccaccaaaaaaaATGCACCCCGTTTTTTTCTAAGCACCAAAGCGACCACCGGTTGATACCCCGTTTTCTGGGTTTTGCGAAGAAAACCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCtt                                                             <  1:6723449/131‑1 (MQ=1)
     ccccaccaaacaaaATGAAACCCATTTTTTTCTAAGCACTCAAGCGGCCCCCATTTGATACCCCGGTTTCTGGGCTTTGCGAAGCAAAACCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCtt                                                             <  1:8022576/131‑1 (MQ=1)
     ccccaccaaacaaaAGACAACCCGTTTTTTTCCAAGAACTCAACTGGCCACCAGTTGATACCCTGGTTTCGGTGCTTGGCGAAGAAAATCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCtt                                                             <  1:18125378/131‑1 (MQ=1)
                          ccgttGTTTTCTAAGCACTCAAGCGGCCACCGGTGGATACCCCCTTTTCTGGGTTTTGCGAAGAAAATCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGcttct                                        <  1:4217212/134‑1 (MQ=1)
                          ccgttGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCATTTTCTGGGTTTTGCGAAGAAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGcttct                                        <  1:9694319/134‑1 (MQ=1)
                             tttttttCTAAGCACTCAATCGGCCACCATTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTgg                                     <  1:18635921/134‑1 (MQ=0)
                             tttttttCTAAGCACTCAAGCGGCCACTAGTGGATACCCCAGTTCCCGGGCTGTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGTTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTgg                                     <  2:6238745/134‑1 (MQ=1)
                             tttttttCGAAGCACGAAAGCGACCACCGGTGGATACCCCAGTTTCTGGGTTTTGCAAAGAAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTgg                                     <  1:11355024/134‑1 (MQ=1)
                             tttttttCGAAACACCAAAGCGGCCACCGGTTGACACCCCAGTTTCTGGGCTTGGAGAAGAAAATCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTgg                                     <  1:904714/134‑1 (MQ=1)
                               gtttcttAAGCACTAAAACGGCCACCAGTTAATACCCCAGTTTCCGGGCTTTGCGAAGAAAAACCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGtt                                   <  1:15979210/131‑1 (MQ=1)
                                tttttgaaGCACTCAAGTGGCCACCGGTTGATACCCCAGTTTCTGGGTTTTGCGAAGAAAACCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTa                                  <  1:4800880/131‑1 (MQ=1)
                                ttttagaaGCACTCAAGCGGCCACAGGTGGATACCCCAGTTTCCGGGCTTTGAGAAAAAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTa                                  <  1:14505639/131‑1 (MQ=1)
                                tttccgaaGCACTCAAGCGGCCACCAGTTAATCCCCCAGTTTCAGGGTTTTGCGAAAAAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTa                                  <  1:663283/131‑1 (MQ=1)
                                  ttttaggCACTCAAGTGGCCACCGTTGGATACCCCATTTTCCGGGCTTGGCGAAACAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAAc                                <  1:16676694/131‑1 (MQ=1)
                                  ttttAAGCCCTCAGGCGGCCACCGGTTGATCCCCCGTTTTCCGGGTTTGGAGAGGCAGATCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAAc                                <  1:7550969/134‑1 (MQ=1)
                                   tcgaaGCCCCCAAGCGACCACCAGTTGATACCCCATTTTCTGGGCTTGGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACa                               <  1:5763514/134‑1 (MQ=1)
                                   tcgaaGCACTCAAGCGGCCACCATTTGATACCCCATTTTCTGGGCTTTGAGAAAAAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACa                               <  1:4022063/134‑1 (MQ=1)
                                   tcgaaGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCTGGGCTTTGCGAAGAAAACCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACa                               <  1:13523886/134‑1 (MQ=0)
                                   tcgaaGCACTCAAGCGGCCACCAGTGGATACCCCGGTTTCCGGGCTTGGAGAAGAAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACa                               <  1:18954358/134‑1 (MQ=1)
                                    ccaaacACTCAACCGGCCACCAGTTGCTACCCCATTTTCCGGGCTTTTCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAg                              <  1:8525504/132‑1 (MQ=1)
                                     gaaacACTCAAGCGGCCACCAGTTGATACCCCATTTTCTGGGCTGTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAgg                             <  1:14672181/133‑1 (MQ=0)
                                     gaaacACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGAAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAgg                             <  1:18374331/133‑1 (MQ=0)
                                     caaacACTCAAGCGGCCACCAGTTGATACCCCATTTTCTGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAgg                             <  1:10977871/133‑1 (MQ=0)
                                      aaacACCCAAGCGGCCACCGGTTGATACCCCAGTTTCTGGGTTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGt                            <  1:14538005/134‑1 (MQ=1)
                                       agaaCTCAAGAGGCCACCAGTTGATACCCCAGTTTCTGGGTTTTGCGAAGAAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGtt                           <  1:2762638/134‑1 (MQ=1)
                                        gcccTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGTGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTg                          >  2:17504001/4‑137 (MQ=1)
                                              aatcGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCAAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACAcc                    <  1:19114539/134‑1 (MQ=1)
                                              aaacGGCCACCAGTTGGTACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACAcc                    <  1:11986056/134‑1 (MQ=1)
                                                    aagccggTTGATACCCCATTTTCTGGTTTTTGCAAAAAAAACCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGt              <  2:1062196/131‑1 (MQ=1)
                                                           tggacaCCCCGTTTTCCGGGCTTGACAAAACAAACCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACt       <  1:5175615/132‑1 (MQ=1)
                                                            tggtACCCCGGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTg      <  2:11939835/134‑1 (MQ=1)
                                                              atcccccATTTTCCGGGCTTGGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGAt    <  2:9905676/134‑1 (MQ=0)
                                                              atcccccAGCTCCAGGGCTTGGCAAAGAAAAACCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGAt    <  1:6790285/134‑1 (MQ=1)
                                                                acccaggTTTCCGGGTTTTGAGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGATgc  <  1:765926/131‑1 (MQ=0)
                                                                      |                                                                                                                                  
TGATCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGATGC  >  AE015451_phaM1/5222924‑5223124

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