Predicted mutation
evidence seq id position mutation annotation gene description
RA AE015451_phaM1 4,740,816 T→G intergenic (+53/+700) gltA CDS (citrate synthase) → / ← yeiW CDS (putative oxidoreductase) –/–

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451_phaM14,740,8160TG95.2% 1064.2 / ‑5.8 311intergenic (+53/+700)gltA CDS (citrate synthase)/yeiW CDS (putative oxidoreductase)–/–
Reads supporting (aligned to +/- strand):  ref base T (0/15);  new base G (166/130);  total (166/145)
Fisher's exact test for biased strand distribution p-value = 7.12e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 311 total aligned reads displayed.

AGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTT  >  AE015451_phaM1/4740690‑4740943
                                                                                                                              |                                                                                                                               
aGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGTAGAAGCGCGACATCGTTGCCCTGAAGGACCGCAAGTAATCTGCAGTAGCAAACTTCAAGAGGGCTGCCCTCGGGGAGTC‑TTTTTTATTTGTCCTGTGCCgg                                                                                                                      >  2:3343955/1‑137 (MQ=255)
aGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAG‑CATTTTTTATTTGTCCTGTGCCgg                                                                                                                      >  2:6387680/1‑137 (MQ=255)
    cGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCctc                                                                                                                  <  1:8454184/137‑1 (MQ=255)
    cGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCctc                                                                                                                  <  2:2063050/137‑1 (MQ=255)
         gTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCACGAAGGCTGCCCTCGAGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTcgc                                                                                                             >  2:5681077/1‑137 (MQ=255)
         gTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTtgc                                                                                                             >  2:8243587/1‑137 (MQ=255)
         gTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTcgc                                                                                                             >  2:9765183/1‑137 (MQ=255)
         gTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTcgc                                                                                                             >  1:6936271/1‑137 (MQ=255)
         gTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTcgc                                                                                                             >  2:6559401/1‑137 (MQ=255)
         atCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGACAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTcga                                                                                                             >  1:15523252/2‑136 (MQ=255)
           ccGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCgg                                                                                                           <  1:10639966/137‑1 (MQ=255)
           ccGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCgg                                                                                                           <  2:9597675/137‑1 (MQ=255)
           ccGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCgg                                                                                                           <  1:14775187/137‑1 (MQ=255)
              cgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCt                                                                                                        >  2:4641452/1‑137 (MQ=255)
               gcCAGCTGTACACCGGCGAGCAGAAGCTCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCtt                                                                                                       <  1:15184623/137‑1 (MQ=255)
               gcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCtt                                                                                                       >  1:13237675/1‑137 (MQ=255)
                         cacCGGCGGGCGGAAGTGCGACATCGTTGCCCGGAAAGACCGCAAGGAGGCTGCATTAGCCGAATGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTTTCCTGTGCTGGCCTCTTCGCGGCCTTGCCCGCGAag                                                                                             <  2:3029576/137‑1 (MQ=255)
                               cGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCg                                                                                       >  1:10862537/1‑137 (MQ=255)
                                gAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCCCTTCGCGGCCTTGCCCGCGAAGAGGCCgg                                                                                      >  1:8121404/1‑137 (MQ=255)
                                 aGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTTTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGc                                                                                     <  1:6203557/137‑1 (MQ=255)
                                 aGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTTCCCGCAAAAAAGCCAGc                                                                                     >  2:4503540/1‑137 (MQ=255)
                                 aGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGGGGACGGc                                                                                     >  1:9936896/1‑137 (MQ=255)
                                 aGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGc                                                                                     >  2:1025342/1‑137 (MQ=39)
                                 aGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGc                                                                                     >  2:14148229/1‑137 (MQ=39)
                                 aGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGc                                                                                     >  2:1126147/1‑137 (MQ=255)
                                 aGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGc                                                                                     >  1:7738160/1‑137 (MQ=39)
                                 aGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGc                                                                                     >  1:2137468/1‑137 (MQ=39)
                                 aGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCCCTTCGCGGCCTTGCCCCCGAAGAGGCCgga                                                                                     >  1:10315970/1‑136 (MQ=255)
                                 aGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGACTCTTCGCGGCCTTGCCCGCGAAGAGGCCgga                                                                                     >  1:14064711/1‑136 (MQ=255)
                                 aGCAGAAGCGCGAAATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCCTGCCCACAAAGAGGcatga                                                                                     >  1:7309984/1‑133 (MQ=255)
                                  gCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGca                                                                                    >  2:11462960/1‑137 (MQ=39)
                                  gCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGca                                                                                    >  2:437718/1‑137 (MQ=255)
                                   cAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAATGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGGAGAGGCCGGcac                                                                                   <  2:13053499/137‑1 (MQ=255)
                                   aaGAAGCGCGACATTGTTGCCCTGAAAGACCGAAAGGAAGCTGCAGTAGCCGAACGAAAGAAGGCTGCCTTTGGGCA‑TCTTTTTTTGTTTTTCTTGTGCCGGCCTCTTTGCGGCCTTTCCCGCGAAGAGGCCGGcac                                                                                   <  1:8243587/136‑1 (MQ=255)
                                    aGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTTTCCGCGAAGAGGCCAGcaca                                                                                  >  1:8747075/1‑137 (MQ=255)
                                    aGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTTCCCGCGAAAAAGCCAGcaca                                                                                  >  1:7401171/1‑137 (MQ=255)
                                    aGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGTCCCCTAAGAGGCCGGaaca                                                                                  >  1:10992167/1‑137 (MQ=255)
                                    aGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGcaca                                                                                  >  2:2522722/1‑137 (MQ=255)
                                    aGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGcaca                                                                                  >  2:8756530/1‑137 (MQ=39)
                                           cgACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGAcacca                                                                           >  2:4199905/1‑137 (MQ=255)
                                              cATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACaa                                                                        >  2:14692871/1‑137 (MQ=255)
                                              cATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACaa                                                                        >  2:2977101/1‑137 (MQ=255)
                                                tCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTTCCCACTAAGAGGCCAGCACAGACAACCCaaaa                                                                      >  1:5874332/1‑137 (MQ=255)
                                                           aaaGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGAAACCACAAAAACTCAGttttt                                                           >  2:7193061/1‑137 (MQ=255)
                                                            aaGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTtc                                                          >  1:3685250/1‑137 (MQ=39)
                                                              gACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTCTTtctc                                                        >  2:537317/1‑137 (MQ=255)
                                                              gACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCCTGCCCACAAAAAGGCCGGCACCGACAACCCCAAAACTCAATTTTtctc                                                        >  1:11241028/1‑137 (MQ=255)
                                                              gACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCCTCGCGGCCCTTCCCGCGAAGAGGCCGGCACAGACACCCCAAAAACCCAGGTTTtctc                                                        >  1:3249051/1‑137 (MQ=255)
                                                                 tgCAAGGAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAG‑CTTTTTTTGTTTGTCTTGTGCTGGCCTTTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTccgc                                                     <  1:9403876/136‑1 (MQ=255)
                                                                 cGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTccgc                                                     <  2:5459772/137‑1 (MQ=255)
                                                                  gCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTccgcc                                                    <  2:7866094/137‑1 (MQ=255)
                                                                  gCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTccgcc                                                    <  1:11708384/137‑1 (MQ=255)
                                                                   cAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGTCTTGCCCGCGAAGAGGCCGGCACCGACACCACAAAAATTCAGTTTTTCTccgccc                                                   >  2:1855445/1‑137 (MQ=38)
                                                                    aaGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCtt                                                                                                                <  2:16004807/75‑1 (MQ=255)
                                                                    aaGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTAGCGGCCTTGCCCGCTAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTccgcccg                                                  >  2:14895792/1‑137 (MQ=255)
                                                                      gTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCt                                                <  2:11781573/137‑1 (MQ=255)
                                                                          gCTGCATTAGCTGAACGCAAGAGGGCTGCCCTCGGGCAGCC‑TTTTTTATTTTTCCTGTGCCGGCCTTNTTGTGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACg                                            <  1:14547435/137‑1 (MQ=255)
                                                                          gCTGCATTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGGGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACg                                            <  2:13369070/137‑1 (MQ=255)
                                                                          gCTGCAGTAGCCGAACGCAAGAAGGCTGCCTTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACg                                            <  2:2410261/137‑1 (MQ=255)
                                                                          gCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTTGGGCAGCC‑TTTTTTATTTTTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACg                                            <  1:7456885/137‑1 (MQ=255)
                                                                          gCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTTTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACg                                            <  1:5785901/137‑1 (MQ=255)
                                                                           ttGCATTAGCCGAACGCAAGAAGGCTGCCCTCGGGAAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCAAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGc                                           <  1:4793636/136‑1 (MQ=255)
                                                                           cTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGc                                           <  2:6228692/137‑1 (MQ=255)
                                                                           cTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGc                                           <  2:13807882/137‑1 (MQ=255)
                                                                           cTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGc                                           <  2:13182356/137‑1 (MQ=255)
                                                                             gCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAg                                         <  2:6749009/137‑1 (MQ=255)
                                                                              cAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑GTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCCTGCCCGCGAAGAGGACGGCACAGACACCACAAAAACTCAGTTTTTCTCCGCCCGCTCACGCAGc                                        >  1:16296485/1‑137 (MQ=255)
                                                                               aGTAGCCAAACGCAAGAAGGCTGCCCTCGGGAAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGTGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGcc                                       <  1:16644370/137‑1 (MQ=255)
                                                                                 tAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACCGACACCACAAAACCTAAGTTTTTCTCCGCCCGCTCACGCAGCCCt                                     >  2:7365392/1‑137 (MQ=255)
                                                                                  aaccGAACGCAAGAAGTCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTTCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCtt                                    <  1:8888558/135‑1 (MQ=255)
                                                                                  aaccGAACGAAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCtt                                    <  1:16461429/135‑1 (MQ=255)
                                                                                  aGCCGAACGCAAGAATGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTTTGCCGGCCCCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCtt                                    <  2:1694341/137‑1 (MQ=255)
                                                                                  aGCCGAACGCAAGAAGGCTGCCCTTGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAGGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCtt                                    <  1:11062844/137‑1 (MQ=255)
                                                                                  aGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCtt                                    <  2:12760455/137‑1 (MQ=255)
                                                                                  aGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGGGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCtt                                    <  2:16546534/137‑1 (MQ=255)
                                                                                  aGCCGAACGCAAGAAGGCTACCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCtt                                    <  1:16048123/137‑1 (MQ=255)
                                                                                  aGCCGAACGAAAGAAGGCTGCCCTCGGGCAG‑CTTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCtt                                    <  1:2991277/137‑1 (MQ=255)
                                                                                  aGCAGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTGTTTGTCTTGTGCCGGCCTCTTCGCGGCTTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCtt                                    <  1:1202944/137‑1 (MQ=255)
                                                                                    ccGAACGCAAGAAGGCTTCCCTCGTGCAGCC‑TTTTTTATTTGCCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTc                                  <  1:663507/137‑1 (MQ=255)
                                                                                    ccGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTTCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTc                                  <  2:16352255/137‑1 (MQ=255)
                                                                                      gAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGACTTGCCCGCGAAGAGGAAGGCACAGACACCCCAAAAACTCAGTTTTTCTCCGCCCGCTCCCACAACCCTTTCAg                                >  1:15255092/1‑137 (MQ=255)
                                                                                       aaCGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCCCTTCTCGGCCTTTCCCACGAAAAGGCCGGCACCAAACCCACAAAAACCCAGTTTTTCTCCGCCCGCTCACACAACCCTTTCAgg                               >  1:13131244/1‑137 (MQ=255)
                                                                                           cAAGAATGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTAt                           <  2:15479449/137‑1 (MQ=255)
                                                                                           cAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGTGGCCTTGCCCGCGAAGAGGCCAGCACAAAAACCACAAAAACTCAGTTTTTCTCCGCCCGACCACGCAGCCCTTTCAGGGTAt                           >  1:11821220/1‑137 (MQ=255)
                                                                                           cAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGTCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGATTACGCAGCCCTTTCAGGGtct                           >  2:12169777/1‑135 (MQ=255)
                                                                                           cAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGTGAAGAGGCCGGCACAGACACCACAAAAACTCAGTTTTTCTCCAACCGCTCACGCAGGACCTTCAGGGTAt                           >  1:12130089/1‑137 (MQ=255)
                                                                                           cAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTAt                           >  2:8081480/1‑137 (MQ=255)
                                                                                           cAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTAt                           <  1:2216864/137‑1 (MQ=255)
                                                                                           cAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAAACCCCGTTTTTCTCCACCCCCCCACGCAACCCTTTCAGGGTAt                           >  2:2241132/1‑137 (MQ=255)
                                                                                            aagaagGCTGCCCTCGGGCATCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTAtt                          <  1:12916374/137‑1 (MQ=255)
                                                                                            aaaaagGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCAGGCTTGCCCGCGAAGAGGAAAGCACAGAAAACACAAAAAACCAGTTTTTCTCCGCCCACCCACACAACCCTTTCAAGGTAtt                          >  1:1349128/4‑137 (MQ=255)
                                                                                             agaagGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTg                         <  2:373376/137‑1 (MQ=255)
                                                                                              gaagGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTTCCCGCGAAGAGGCCGGCACAGACACCCCAAAAACTCAGTTTTTCTCCGCCCGCTCCCGCAGCCCTTTTAAGGTGTTGa                        >  1:10098658/1‑137 (MQ=255)
                                                                                              gaagGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTTCCCGCGAAAAGGCCGGCACAGACACCCCAAAAACTCAGTTTTTCTCCGCCCACTCACGCAGCCCTTTCCGGGtctgtc                        >  2:7583983/1‑132 (MQ=255)
                                                                                                 ggCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCCCTTCGCGGCCTTGCCCGCGCAGAGGCCGGAAAAGACACCAAAAAACCTCAGTTTTTCTCCACCCGCCCACACAGCCCCTTCGAGGTGTTGAACg                     >  2:1641997/1‑137 (MQ=255)
                                                                                                   cTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGgc                   <  2:6072572/137‑1 (MQ=255)
                                                                                                         tCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGAAACCACAAAAACTCAGTTTTTCTCCGCCCGCTCCCACAGCACTTTCAGGGTATTGAACGGCGCATcc             >  2:10558765/1‑137 (MQ=37)
                                                                                                          cGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGTAGCCCTTTCAGGGTATTGAACGGCGCATcca            >  1:14398180/1‑137 (MQ=39)
                                                                                                                agccttttttATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACg       <  2:3639706/133‑1 (MQ=255)
                                                                                                                    ttttttATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAAGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCCTCCCCCACGAACt   >  2:5165859/1‑137 (MQ=255)
                                                                                                                              |                                                                                                                               
AGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTT  >  AE015451_phaM1/4740690‑4740943

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