Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451_phaM15,485,1530AG32.1% 70.0 / 12.5 55intergenic (+59/‑88)purH CDS (phosphoribosylaminoimidazolecarboxamide formyltransferase_IMP cyclohydrolase)/purD CDS (Phosphoribosylamine‑‑glycine ligase)–/–
Reads supporting (aligned to +/- strand):  ref base A (7/29);  new base G (0/17);  total (7/48)
Fisher's exact test for biased strand distribution p-value = 8.25e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTG‑CCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTG  >  AE015451_phaM1/5485017‑5485285
                                                                                                                                         |                                                                                                                                    
gATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTG‑CCCGCGAATACGa                                                                                                                                      <  2:6373647/137‑1 (MQ=255)
 aTGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTG‑CCCGCGAATACGAt                                                                                                                                     <  2:1726810/137‑1 (MQ=255)
 aTGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTG‑CCCGCGAATACGAt                                                                                                                                     <  2:13312235/137‑1 (MQ=255)
 aTGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTG‑CCCGCGAATACGAt                                                                                                                                     <  2:3016852/137‑1 (MQ=255)
    cTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTG‑CCCGCGAATACGAtggt                                                                                                                                  <  1:4358318/137‑1 (MQ=255)
    cTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTG‑CCCGCGAATACGAtggt                                                                                                                                  <  1:4358357/137‑1 (MQ=255)
         gTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTG‑CCCGCGAATACGATGGTGGAtt                                                                                                                             <  1:5550065/137‑1 (MQ=255)
             gcGCTGCGGACGACGAGGCCGGCATCGCGATGTTCTTCACTGGCAGGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTG‑CCCGCGAATACGATGGTGGATTCAcc                                                                                                                         <  1:15861517/136‑1 (MQ=255)
               gctgctGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGGTTG‑TCCGCGAACAAAAAGGTGGAGTccccac                                                                                                                       >  2:16414584/1‑135 (MQ=255)
                   ctgcCGGCGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTG‑CCCGCGAATACGATGGTGGATTCACCGCCGCa                                                                                                                   <  2:6361515/137‑1 (MQ=255)
                   ctgcCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGNTCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTG‑CCCGCGAATACGATGGTGGATTCACCGCCGCa                                                                                                                   <  1:14552591/137‑1 (MQ=255)
                   ctgcCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGATGCAAGTCGAGATCCTGTGGGAGCGGGCTTG‑CCCGCGAATACGATGGTGGATTCACCGCCGCa                                                                                                                   <  2:16713316/137‑1 (MQ=255)
                   ctgcCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTG‑CCCGCGAATACGATGGTGGATTCACCGCCGCa                                                                                                                   <  1:10772264/137‑1 (MQ=255)
                   ctgcCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTG‑CCCGCGAATACGATGGTGGATTCACCGCCGCa                                                                                                                   <  2:1495675/137‑1 (MQ=255)
                   ctgcCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTG‑CCCGCGAATACGATGGTGGATTCACCGCCGCa                                                                                                                   <  1:11924020/137‑1 (MQ=255)
                   ctgcCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTG‑CCCGCGAATACGATGGTGGATTCACCGCCGCa                                                                                                                   <  1:3335484/137‑1 (MQ=255)
                   ctgcCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTG‑CCCGCGAATACGATGGTGGATTCACCGCCGCa                                                                                                                   <  2:6361547/137‑1 (MQ=255)
                    tgcCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTG‑CCCGCGAATACGATGGTGGATTCACCGCCGCAt                                                                                                                  <  2:15566965/137‑1 (MQ=255)
                              gCCGGCATCGCGATGGTCTTCATGGGCATGCACCACTTCCGCCACTAATTACGCGGATCCTGAGGCAAGTCGAGATCTTGTGGGAGCGGGCTTG‑CCCGCGAATACGAGGGTGGATTCACCGCCGCATTCGCGGACaa                                                                                                        <  2:14662539/137‑1 (MQ=255)
                                             tccTTCACTGTCATCCCCCATTTCCCCCACCAACTACGCGGATCCCGAAGAAATTCGAAGCCCGGGGGAAGTGGGCTTT‑CCCGAGACCCAGGGGGTGAATCCACCGCCCCGTCAGAGGGCAAGCCCGCTCCCACAgt                                                                                         <  2:1939790/135‑1 (MQ=255)
                                                    ttGGCGTGCGCCATTTCCGCCACCAATTACACGGACCCCGAGACAAGTTGAGATCCGGTGGGGGTGGGTTTC‑CCCGAGAGTAAGGGGGTGGATCCACCGCCGCATTCGCGGGCAACCCCGCTCCCACAGTGTTCGgc                                                                                  <  2:3612244/136‑1 (MQ=255)
                                                    ttGGCATGCACCCTTTCCGCCACTAATTACGTGAATCCCGGGGAAGGTTGAGATCTTGTGGGAGTGTGTTTG‑CCCGAAAATATGGTGGGGGATTCCCCCCCGCATTTGCGGGAAACCCCGCTCCCACAGTGTTCGgc                                                                                  <  2:3612245/136‑1 (MQ=255)
                                                    tgggCGTGCGCCATTTCCGCCACGAATTACGCGGACCCCGGGGAAAGTCGAGATCTTGGGGGAG‑GGTATTGACCCGCGAATATGAGGGGGGATTCACCCCCCCATTCGAGGGAAATCCCGCTCCCACAGTGTTCGgc                                                                                  <  2:14560158/135‑1 (MQ=255)
                                                        cATGCACCACTTCCACCACTAATGCCGCGGACCCCGAGGAAGGTGGAGGTCGGGGGGGGGGGGGTTTG‑CCCGAGAGGAGGGGGGTGGATTCCCCGCCGCATTCGGGGGAAACCCCGCTCCCACAGTGTTCGGCGCaa                                                                              <  2:3645490/137‑1 (MQ=255)
                                                           gcgcCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTG‑CCCGAGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGcc                                                                           >  2:10854240/1‑137 (MQ=255)
                                                             gcCATTTCCCCCACAGCTTACGTGGATCCCGAGGAAATTAGGGATCTTGGGGGAGGGGGCTTC‑CCCGAGAATACGGTGTGGGTTCCACCCCCACGTTCGGGGGAAACCCCGCTCCCACAGTGTTCGGCGCAAGCCAc                                                                         <  2:2864431/137‑1 (MQ=255)
                                                               cactcacGCCATTAATCAAGCGAGTCACGGGACAATTCGAGATCTTGGGGGAGCGGGTTTC‑CCCACGAATACGTTGGTGGCTCCACCTCATCCTCCGCGAGCAATCCCGCTCCCACAGTGTTCGGCGCAAGCCACgg                                                                       <  2:2937338/131‑1 (MQ=255)
                                                                aCTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTG‑CCCGAGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACggg                                                                      >  1:15412930/1‑137 (MQ=255)
                                                                  ttCCCCCACTACTTACTCGGCCCCCGGGACAGGTCGAGATCTTGGGGGAGTGGTTTTA‑CCCGGAAGGAGGGGGGGGGGTTCACCCCCGCGTTCGCGGGCAACCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGAt                                                                    <  2:11456919/137‑1 (MQ=255)
                                                                                              acAAGTCGAGCTCTTGTGGGAGTGGTCTTC‑CCCGCAAAGACGGTGGGGGATTCACCGCCGCATTCGCGGGAAACCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAg                                        <  2:16479430/136‑1 (MQ=255)
                                                                                                   tCGAGCTCTTGTGGGAGTGGGCTTG‑CCCGCGAATACGAGGGGGGATTCCCCCCCGCGTTTACGGACAACCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTg                                   <  2:3721033/137‑1 (MQ=255)
                                                                                                         tCCTGTGGGAGCGGGCTTG‑CCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGATCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGatcatc                             <  2:10401014/137‑1 (MQ=255)
                                                                                                             gtgGGAGCGGTCTTA‑CCCACGAATACTGGGTTTACTTCCCCTCCCCATTCACGGGCAAACCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCa                         <  2:13489004/137‑1 (MQ=255)
                                                                                                             gggggAGCGGGCTTA‑CCCGCAAATACGATGGGGGATCCACCGCCGCATCCGTGGACAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCa                         <  2:3587470/135‑1 (MQ=255)
                                                                                                             gggggAGCGGGATGC‑CCCGAGAAGAAGGTGGGGGATCCACCGCCACGTTCGTGGGAAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCa                         <  2:4947191/135‑1 (MQ=255)
                                                                                                              ggggAGCGGTCTTG‑CCCGAGAATATGATGGTGGATCCCCCGCCGCATTCGCGGGCAATCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAg                        <  2:15728689/136‑1 (MQ=255)
                                                                                                                   gggttcTTC‑CCCGCGAAGAGGGGGGTGGCTCCCCCGCCGTGTTCGGGGACAATCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTg                   <  1:14389049/132‑1 (MQ=255)
                                                                                                                    cgtgCTCG‑CCCGAAAAGAAGATGGGGGATCCACCGCCACGTCCGCGGAAAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTgg                  <  2:11356584/134‑1 (MQ=255)
                                                                                                                        tttG‑CCAGAGAAGAGGGTGGGGGCTTCACCTCTGCATTCTCGGGCAACCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGt              <  2:1303957/136‑1 (MQ=255)
                                                                                                                        cTTG‑CCCGAGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGt              <  2:15412930/137‑1 (MQ=255)
                                                                                                                         tttacccGAGAATAAGAGGGGGGATTCACCGCCGCGTTTGGGGGAAACCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGt              <  2:6190008/133‑1 (MQ=255)
                                                                                                                         ttG‑CCCGAAAATAAGGTGGGGGGTCCCCCGCCACGTTCGCGGAAAACCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTg             <  2:6129218/137‑1 (MQ=255)
                                                                                                                          ttccccGCGCATACGGTGTGGGGTCCCCCCCCACTTTCGCGGGCACCCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTg             <  2:4333952/134‑1 (MQ=255)
                                                                                                                          ttccccACGAATACGGGGTGGGATCCACCCCCACCTTCGCGGGAAATCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTg             <  2:14954987/134‑1 (MQ=255)
                                                                                                                               cgcgGATACGTTGGTGGATTCACCGCCGCATTCGCGGACAACCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCAc        <  2:3953269/137‑1 (MQ=255)
                                                                                                                               cgcgAATACGAGGGTGGGTCCACCGCCATATTCGTGAGAAATCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCAc        <  2:15711647/137‑1 (MQ=255)
                                                                                                                               cgcaAATATGATGTTGACTTCCCCTCCCCCTCCGCGAACAATCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCAc        <  2:8667084/137‑1 (MQ=255)
                                                                                                                               cgagAATACGGTGGTGGGTTCACCGCCGCGTTCGGGGGCAAACCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCAc        <  2:12491233/134‑1 (MQ=255)
                                                                                                                               cgaaaaTACGATGTGGGATCCACCCCCACTTTCGCGGAAAATCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCAc        <  2:7843783/133‑1 (MQ=255)
                                                                                                                               cagaaaTACGGTGTGGGATTCACCGCCGCATTAGTGGGCAACCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCAc        <  2:15711648/133‑1 (MQ=255)
                                                                                                                               cacgAATACGGTGGTGGCTCCCCCCCCGACTTCTCGGAAAACCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCAc        <  1:11555648/135‑1 (MQ=255)
                                                                                                                               cacgAATAAGAGGTTGGATTCCCTTCCGCATTCGCGGAAAACCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCAc        <  2:11471970/135‑1 (MQ=255)
                                                                                                                                   aaTACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGccc    >  2:13486942/1‑137 (MQ=255)
                                                                                                                                     tACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctg  >  1:9674224/1‑137 (MQ=255)
                                                                                                                                     tACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctg  >  1:2587236/1‑137 (MQ=255)
                                                                                                                                     tACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctg  >  1:13185579/1‑137 (MQ=255)
                                                                                                                                         |                                                                                                                                    
GATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTG‑CCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTG  >  AE015451_phaM1/5485017‑5485285

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