Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451_phaM15,223,0130CA47.1% ‑6.2 / 28.8 51L428F (TTG→TTT) PP_4603 CDS (transposase)
Reads supporting (aligned to +/- strand):  ref base C (1/26);  new base A (0/24);  total (1/50)
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.

TGATCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGATGCGCAGTGGCCTTGATGG  >  AE015451_phaM1/5222924‑5223140
                                                                                         |                                                                                                                               
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)
                                                                             ggttttGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGATGCGCAGTGGCCTTGa     <  2:2510010/134‑1 (MQ=1)
                                                                             ggtttggCGAAGAAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGATGCGCAGTGGCCTTGa     <  2:257078/131‑1 (MQ=0)
                                                                               ctttaagAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGATGCGCAGTGGCCTTGATg   <  1:17206109/131‑1 (MQ=1)
                                                                                tttgaaaaGAAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGATGCGCAGTGGCCTTGATgg  <  1:15179382/131‑1 (MQ=0)
                                                                                ttggCGAAGCAAATCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGATGCGCAGTGGCCTTGATgg  <  1:7776658/134‑1 (MQ=1)
                                                                                         |                                                                                                                               
TGATCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAACTGATGCGCAGTGGCCTTGATGG  >  AE015451_phaM1/5222924‑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: ‑