Predicted mutation
evidence seq id position mutation annotation gene description
RA AE015451_phaM1 4,740,811 (T)7→6 intergenic (+48/+705) 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,8050T.95.2% 1362.9 / 44.2 299intergenic (+42/+711)gltA CDS (citrate synthase)/yeiW CDS (putative oxidoreductase)–/–
Reads supporting (aligned to +/- strand):  ref base T (0/14);  new base . (155/126);  total (158/141)
Fisher's exact test for biased strand distribution p-value = 2.05e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 377 total aligned reads displayed.

AATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCAC  >  AE015451_phaM1/4740670‑4740937
                                                                                                                                       |                                                                                                                                    
aaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCtt                                                                                                                                     <  1:2628207/137‑1 (MQ=255)
   gCTCTCCAGCCCTTACAAGATCGGCCGTCCGCGCCAGCCGTACACCTGCGAGCAGAAGCGCGACATCGTTGCCCTGAAATACCGCAAGTAAGCTGCAGTAGCCAAACGCAAGAAGGCTGCCCTCGGGCAGCCttttt                                                                                                                                  <  2:13369342/137‑1 (MQ=255)
   gCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttt                                                                                                                                  <  2:11085109/137‑1 (MQ=255)
   gCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttt                                                                                                                                  <  1:8010447/137‑1 (MQ=255)
   gCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttt                                                                                                                                  <  1:6886337/137‑1 (MQ=255)
   gCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttt                                                                                                                                  <  1:15753873/137‑1 (MQ=255)
   gCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttt                                                                                                                                  <  2:400611/137‑1 (MQ=255)
   gCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttt                                                                                                                                  <  1:13337297/137‑1 (MQ=255)
   gCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCGGCCCTCGGGCAGCCttttt                                                                                                                                  <  1:4311213/137‑1 (MQ=255)
   gCTCTCCAGCCCGTACAAGAGCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGTCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAATGCTGCCCTCTGGCAGCCttttt                                                                                                                                  <  2:6257935/137‑1 (MQ=255)
       tcCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttatt                                                                                                                              >  2:15773231/1‑134 (MQ=255)
       tcCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttatt                                                                                                                              >  2:5876063/1‑134 (MQ=255)
        ccAGCCCGTACAAGATCGGCCTTCCGCGCCAGCTGTACACCGGCGAGTAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGTCAGCCttttttattt                                                                                                                             <  2:16664566/137‑5 (MQ=255)
        ccAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttattt                                                                                                                             <  1:2062471/137‑5 (MQ=255)
        ccAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttattt                                                                                                                             <  2:14031237/137‑5 (MQ=255)
        ccAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttattt                                                                                                                             <  1:14613532/137‑5 (MQ=255)
        ccAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttattt                                                                                                                             <  2:9772249/137‑5 (MQ=255)
        ccAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttattt                                                                                                                             <  1:13174268/137‑5 (MQ=255)
         cAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttatttg                                                                                                                            >  2:5610489/1‑132 (MQ=255)
                        cGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCctc                                                                                                            >  2:1276862/1‑137 (MQ=255)
                        cGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCctc                                                                                                            <  2:2063050/137‑1 (MQ=255)
                        cGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCctc                                                                                                            <  1:8454184/137‑1 (MQ=255)
                         ggtcGTCCGCGCCAGCTGTACACCGTCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCctct                                                                                                           <  1:2192430/134‑1 (MQ=255)
                             gTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTcgt                                                                                                       >  1:10663872/1‑136 (MQ=255)
                             gTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTcgc                                                                                                       >  2:83044/1‑137 (MQ=255)
                             gTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTcgc                                                                                                       >  1:3592783/1‑137 (MQ=255)
                             gTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTcgc                                                                                                       >  2:3093765/1‑137 (MQ=255)
                              tCCGCGCCAGCTGTACACCGTCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCATTAGCCGAACGCAAGAAGGCTGCCCTCGTGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTcgcg                                                                                                      <  1:14702830/137‑1 (MQ=255)
                                         tGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAGGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCcgc                                                                                           <  2:13246498/137‑1 (MQ=255)
                                                cGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAgg                                                                                    >  2:2673757/1‑137 (MQ=255)
                                                  gCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGcc                                                                                  <  2:15013173/137‑1 (MQ=255)
                                                     aGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTTCCGGCCCTGCCCGCGAAGAGGCCGGc                                                                               >  1:4879223/1‑137 (MQ=255)
                                                     aGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGc                                                                               >  1:7738160/1‑137 (MQ=39)
                                                     aGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGc                                                                               >  2:7220548/1‑137 (MQ=39)
                                                     aGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGACTCTTCGCGGCCTTGCCCGCGAAGAGGCCgga                                                                               >  1:14064711/1‑136 (MQ=255)
                                                      gCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAGGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGca                                                                              >  2:11758034/1‑137 (MQ=38)
                                                      gCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGca                                                                              >  2:437718/1‑137 (MQ=255)
                                                       cAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGcac                                                                             <  1:14866564/137‑1 (MQ=255)
                                                       cAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGcac                                                                             <  2:6957632/137‑1 (MQ=255)
                                                       aaGAAGCGCGACATCTTTGCCCTGAAAGACCACAAGGAAGCTGCAGTAGCCGAAAAAAAGAAGGCTGCCCTCGGGCAGCCTTTTTTGTTTTTCCT‑GTGCCGGCCTCTTCGTGGCCTTGCCCGCGAAGAGGCCGGcac                                                                             <  1:11049591/136‑1 (MQ=255)
                                                       aaGAAGCGCGACATCGTTGCCCTGAAAGACAGCAAGGAAGCTGCAGGAGCCGAACGCAAGAAGGCTGCCCTCGGGCA‑CCCTTTTTGATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGcac                                                                             <  1:7342263/136‑1 (MQ=255)
                                                        aGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTTTCCGCGAAGAGGCCAGcaca                                                                            >  1:8747075/1‑137 (MQ=255)
                                                        aGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGcaca                                                                            >  2:9464496/1‑137 (MQ=39)
                                                        aGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGcaca                                                                            >  1:9800336/1‑137 (MQ=39)
                                                        aGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGcaca                                                                            >  1:8203979/1‑137 (MQ=255)
                                                        aGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGACTTTCCCGCGAAGAGGCCGGcaaa                                                                            >  1:5471083/1‑135 (MQ=255)
                                                                  cATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGTCGGAACAAACACCACaa                                                                  >  1:10947560/1‑137 (MQ=255)
                                                                  cATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACaa                                                                  >  2:2977101/1‑137 (MQ=255)
                                                                  cATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACaa                                                                  >  2:14692871/1‑137 (MQ=255)
                                                                    tCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGGAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACaaaa                                                                <  2:6423570/137‑1 (MQ=255)
                                                                        tGCCCTGAAAGATCGCAAGGAAGCGGCAGGAGCTGAACGCAAGAAGGCTTCCCTCGGGC‑GTCTTTTTTTGTTTGTCCTTTTCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTc                                                            <  1:11508297/137‑1 (MQ=25)
                                                                                aaGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTG‑CGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTtct                                                   <  2:10110443/137‑1 (MQ=255)
                                                                                 aGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGTCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTtct                                                   <  2:6198052/137‑1 (MQ=255)
                                                                                 aGACCGCAAGGAGGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTGTTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTtct                                                   <  1:10497378/137‑1 (MQ=255)
                                                                                  gACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCCTGCCCACAAAAAGGCCGGCACCGACAACCCCAAAACTCAATTTTtctc                                                  >  1:11241028/1‑137 (MQ=255)
                                                                                  gACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCCTCGCGGCCCTTCCCGCGAAGAGGCCGGCACAGACACCCCAAAAACCCAGGTTTtctc                                                  >  1:3249051/1‑137 (MQ=255)
                                                                                     cGCAGGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTTCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTccgc                                               <  1:896588/137‑1 (MQ=255)
                                                                                     cGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCTGGCCTTTTCGCGGCCTTGCCCGTGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTccgc                                               <  1:6634744/137‑1 (MQ=255)
                                                                                     cGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTccgc                                               <  2:13065195/137‑1 (MQ=255)
                                                                                     cGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTccgc                                               <  2:16271819/137‑1 (MQ=255)
                                                                                       cAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGTCTTGCCCGCGAAGAGGCCGGCACCGACACCACAAAAATTCAGTTTTTCTccgccc                                             >  2:1855445/1‑137 (MQ=38)
                                                                                        aaGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGGCCAGACACCACAAAAACTCAGTTTTTCTccgcccg                                            >  2:3691853/1‑137 (MQ=38)
                                                                                          ggaaGCTGAAGTAGCCGAACGAAAGAAGGTTTCCCTTGGGC‑GGCTTTTTTTGTTTTTCTTTTTCTGGCTTTTTCGCGGCCTTTCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCt                                          <  1:5876063/135‑1 (MQ=255)
                                                                                           tAAACTGCAGTAGCCAAACACAAGAAGGCTGCCCTCGTGCATCC‑TTTTTTATTTGTCCTGGGCTGGCCTCTTCGTGGCCTGGCCCGCGAAAAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTc                                         <  1:2935239/137‑1 (MQ=255)
                                                                                            aaGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCa                                        <  1:15200704/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)
                                                                                              gCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACg                                      <  2:12397054/137‑1 (MQ=255)
                                                                                               ctacAGGAGCCGAACGCAAGAAGGCTGCCCGCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGc                                     <  2:3591771/134‑1 (MQ=255)
                                                                                               cTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGc                                     <  2:4611757/137‑1 (MQ=255)
                                                                                                  cAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑GTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCCTGCCCGCGAAGAGGACGGCACAGACACCACAAAAACTCAGTTTTTCTCCGCCCGCTCACGCAGc                                  >  1:16296485/1‑137 (MQ=255)
                                                                                                    gTAGCCAAACACAAGAAGGCTGCCCTCGGGCAACC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTCCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGccc                                <  1:15412956/137‑1 (MQ=255)
                                                                                                     tAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCCTTCCCCCTAAGAGGCCAGCCCAAAAAACACAAAAACTCAATTTTTCTCCGCCCACTCCCACCGCCCt                               >  1:1292499/1‑137 (MQ=255)
                                                                                                     tAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑‑TTTTTATTTTTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTCTCTCCGCCCGCTCACGCAGCCCtt                              <  2:7942510/137‑1 (MQ=39)
                                                                                                      aGCCGAACGCAAGAAGGCTGCCTTTGGGCAG‑CTTTTTTTGTTTGTCTTGTGCCGGCCTTTTCGCGGCCTTTCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCtt                              <  1:9287698/137‑1 (MQ=255)
                                                                                                      aGCCGAACGCAAGAAGGCTACCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCtt                              <  1:16048123/137‑1 (MQ=255)
                                                                                                      aGCCAAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCTTCTTCGTGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCtt                              <  1:725606/137‑1 (MQ=255)
                                                                                                      aGCAGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTGTTTGTCTTGTGCCGGCCTCTTCGCGGCTTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCtt                              <  1:1202944/137‑1 (MQ=255)
                                                                                                        ccGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTc                            <  1:16486870/137‑1 (MQ=255)
                                                                                                        ccGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTc                            >  2:6709775/1‑137 (MQ=255)
                                                                                                          gAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAg                          >  2:5240039/1‑137 (MQ=39)
                                                                                                               cAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGTGGCCTTGCCCGCGAAGAGGCCAGCACAAAAACCACAAAAACTCAGTTTTTCTCCGCCCGACCACGCAGCCCTTTCAGGGTAt                     >  1:11821220/1‑137 (MQ=255)
                                                                                                               cAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGTGAAGAGGCCGGCACAGACACCACAAAAACTCAGTTTTTCTCCAACCGCTCACGCAGGACCTTCAGGGTAt                     >  1:12130089/1‑137 (MQ=255)
                                                                                                               cAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTAt                     >  2:8081480/1‑137 (MQ=255)
                                                                                                               cAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAAACCCCGTTTTTCTCCACCCCCCCACGCAACCCTTTCAGGGTAt                     >  2:2241132/1‑137 (MQ=255)
                                                                                                                aaaaagGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCAGGCTTGCCCGCGAAGAGGAAAGCACAGAAAACACAAAAAACCAGTTTTTCTCCGCCCACCCACACAACCCTTTCAAGGTAtt                    >  1:1349128/4‑137 (MQ=255)
                                                                                                                  gaagGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTTCCCGCGAAGAGGCCGGCACAGACACCCCAAAAACTCAGTTTTTCTCCGCCCGCTCCCGCAGCCCTTTTAAGGTGTTGa                  >  1:10098658/1‑137 (MQ=255)
                                                                                                                  gaagGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTTCCCGCGAAAAGGCCGGCACAGACACCCCAAAAACTCAGTTTTTCTCCGCCCACTCACGCAGCCCTTTCCGGGtctgtc                  >  2:7583983/1‑132 (MQ=255)
                                                                                                                   aagGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGACGAGGCCGTCCCAGACACCACAAAAACTCCGTTTTTCTCCGCACGCTCACGCAGCACTTTCAGGGTATGGaa                 >  2:13776029/1‑137 (MQ=255)
                                                                                                                   aagGCTGCCCTCGGGCAGCC‑TTTTTGATTTGTCCGGTGCTGGCCTCTTTGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGaa                 <  2:4676437/137‑1 (MQ=255)
                                                                                                                     ggCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCCCTTCGCGGCCTTGCCCGCGCAGAGGCCGGAAAAGACACCAAAAAACCTCAGTTTTTCTCCACCCGCCCACACAGCCCCTTCGAGGTGTTGAACg               >  2:1641997/1‑137 (MQ=255)
                                                                                                                       cTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGgc             <  2:16856745/137‑1 (MQ=255)
                                                                                                                           ccTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGGCAGCACAGACACCACAAAACCTCCATTTTTCTCCGCCCACTCACACAAACCCTTCAAGGTATTTAACAGCGCac         >  1:14476687/1‑136 (MQ=255)
                                                                                                                            cTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGAAACCACAAAAACCCAGTTTTTCTCCACCCACTCACGCAACCCCTTCAAGGTATTGAACGGCGCATc        >  1:10244234/1‑137 (MQ=37)
                                                                                                                             tCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAGGAGGCCGGCGCAGACACCACAAAAACCCAGTTTTTCGCGGCCCGATCACGCAGCCCTTTCCGGGTATTGAACGGAGCATcc       >  2:11288438/1‑137 (MQ=255)
                                                                                                                             tCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATcc       >  2:5457281/1‑137 (MQ=255)
                                                                                                                             tCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGAAACCACAAAAACTCAGTTTTTCTCCGCCCGCTCCCACAGCACTTTCAGGGTATTGAACGGCGCATcc       >  2:10558765/1‑137 (MQ=37)
                                                                                                                              cGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATcca      <  2:15911299/137‑1 (MQ=255)
                                                                                                                               gggCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGGGGCCTTGCCCGCGAAGAGGCAGGTACAGACACCACAAAACCTCAGTTTTTCTCCGCCCTATCAAGCAGCCCTTTCAGGGTATTGAAAGTCGCCTccac     >  2:11948754/1‑137 (MQ=255)
                                                                                                                                       |                                                                                                                                    
AATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCAC  >  AE015451_phaM1/4740670‑4740937

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