Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,451,422 A→T 100% T246S (ACG→TCG)  paaE → ring 1,2‑phenylacetyl‑CoA epoxidase, NAD(P)H oxidoreductase component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,451,4220AT98.6% 953.5 / 1.4 287T246S (ACG→TCG) paaEring 1,2‑phenylacetyl‑CoA epoxidase, NAD(P)H oxidoreductase component
Reads supporting (aligned to +/- strand):  ref base A (3/1);  new base T (137/146);  total (140/147)
Fisher's exact test for biased strand distribution p-value = 3.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 287 total aligned reads displayed.

TACAGTCACTTGGGGCCTCGCTCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGCGGCATTGCGCCAGGGGGCGGA  >  CP009273/1451276‑1451567
                                                                                                                                                  |                                                                                                                                                 
tACAGTCACTTGGGGCCTCGCTCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCg                                                                                                                                                 <  1:919633/149‑1 (MQ=255)
 aCAGTCACTTGGGGCCTCGCTCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGc                                                                                                                                                >  1:210227/1‑149 (MQ=255)
  cAGTCACTTGGGGCCTCGCTCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGcc                                                                                                                                               <  1:1293394/149‑1 (MQ=255)
    gTCACTTGGGGCCTCGCTCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTg                                                                                                                                             <  1:3415070/149‑1 (MQ=255)
          tGGGGCCTCGCTCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCAcgc                                                                                                                                       <  2:3267839/149‑1 (MQ=255)
          tGGGGCCTCGCTCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCAcgc                                                                                                                                       <  2:3208593/149‑1 (MQ=255)
              gCCTCGCTCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGt                                                                                                                                   <  2:3900739/149‑1 (MQ=255)
              gCCTCGCTCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGt                                                                                                                                   >  2:1506987/1‑149 (MQ=255)
              gCCTCGCTCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGt                                                                                                                                   >  1:2289405/1‑149 (MQ=255)
               ccTCGCTCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTc                                                                                                                                  <  1:1091460/149‑1 (MQ=255)
                 tCGCTCATTAATTTTCGTCTTTATGATTAGGCATTTATTTGTGGTCCGTCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAAACATTAATCTGTAGAGGTTTAATTCGCCTCGCAAGCGCGTCaa                                                                                                                                >  2:1917906/1‑149 (MQ=255)
                  cGCTCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCaaa                                                                                                                               >  2:3138607/1‑149 (MQ=255)
                     tCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGt                                                                                                                            <  1:682568/149‑1 (MQ=255)
                        ttAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGc                                                                                                                         >  2:2216053/1‑149 (MQ=255)
                           aTTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGtt                                                                                                                      >  2:870532/1‑149 (MQ=255)
                           aTTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGtt                                                                                                                      <  1:988630/149‑1 (MQ=255)
                             tttCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTaa                                                                                                                    <  1:17821/149‑1 (MQ=255)
                              ttCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAAc                                                                                                                   >  1:3443538/1‑149 (MQ=255)
                                 gTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTCGCGTTAACGTg                                                                                                                >  1:4033277/1‑149 (MQ=255)
                                    tttATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCaa                                                                                                             >  2:2722747/1‑149 (MQ=255)
                                     ttATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCaaa                                                                                                            <  1:31828/149‑1 (MQ=255)
                                       atgatgAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGt                                                                                                          <  1:1416497/149‑1 (MQ=255)
                                          atgaGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGAc                                                                                                       <  1:2300977/149‑1 (MQ=255)
                                           tgaGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACg                                                                                                      >  2:14704/1‑149 (MQ=255)
                                            gaGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACgg                                                                                                     <  1:2992710/149‑1 (MQ=255)
                                            gaGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACgg                                                                                                     <  1:2991953/149‑1 (MQ=255)
                                            gaGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACgg                                                                                                     >  1:2265882/1‑149 (MQ=255)
                                             aGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGa                                                                                                    <  2:1895103/149‑1 (MQ=255)
                                             aGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGa                                                                                                    >  2:3117672/1‑149 (MQ=255)
                                              ggCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGAc                                                                                                   >  2:3587749/1‑149 (MQ=255)
                                               gCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACa                                                                                                  <  2:2601372/149‑1 (MQ=255)
                                                cATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACaa                                                                                                 >  2:1303032/1‑149 (MQ=255)
                                                cATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACaa                                                                                                 >  2:2710161/1‑149 (MQ=255)
                                                  tttatttGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACaaaa                                                                                               >  2:1976587/1‑149 (MQ=255)
                                                  tttatttGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACaaaa                                                                                               <  2:3403717/149‑1 (MQ=255)
                                                  tttatttGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACaaaa                                                                                               >  1:204673/1‑149 (MQ=255)
                                                  tttatttGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACaaaa                                                                                               <  2:3403669/149‑1 (MQ=255)
                                                  tttatttGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACaaaa                                                                                               >  2:4018763/1‑149 (MQ=255)
                                                  tttatttGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACaaaa                                                                                               >  1:2864663/1‑149 (MQ=255)
                                                    tatttGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAg                                                                                             >  1:1313996/1‑149 (MQ=255)
                                                     atttGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGt                                                                                            >  1:3814563/1‑149 (MQ=255)
                                                     atttGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGt                                                                                            <  1:2784183/149‑1 (MQ=255)
                                                     atttGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGt                                                                                            <  2:2503692/149‑1 (MQ=255)
                                                            gTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTa                                                                                     <  1:3488477/149‑1 (MQ=255)
                                                            gTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTa                                                                                     >  2:443131/1‑149 (MQ=255)
                                                             tCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTAc                                                                                    >  2:2860043/1‑149 (MQ=255)
                                                             tCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTAc                                                                                    >  2:1542357/1‑149 (MQ=255)
                                                               cggcggCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGt                                                                                  <  2:266571/149‑1 (MQ=255)
                                                                 gcggcgATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCa                                                                                <  2:671355/149‑1 (MQ=255)
                                                                     cgATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGAt                                                                            <  2:3392432/149‑1 (MQ=255)
                                                                       atgatgGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATgg                                                                          <  1:739953/149‑1 (MQ=255)
                                                                       atgatgGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATgg                                                                          <  1:3983910/149‑1 (MQ=255)
                                                                       atgatgGATGAAGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATgg                                                                          >  2:3471553/1‑149 (MQ=255)
                                                                           tggatgACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCgg                                                                      <  2:1661552/149‑1 (MQ=255)
                                                                             gatgACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGa                                                                    <  1:335411/149‑1 (MQ=255)
                                                                             gatgACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGa                                                                    >  2:1827190/1‑149 (MQ=255)
                                                                                gACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCTAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCg                                                                 <  2:1749230/149‑1 (MQ=255)
                                                                                      gAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAAt                                                           >  2:2081248/1‑149 (MQ=255)
                                                                                      gAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAAt                                                           >  1:2286137/1‑149 (MQ=255)
                                                                                       aaaCCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATc                                                          >  1:1850691/1‑149 (MQ=255)
                                                                                       aaaCCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATc                                                          >  1:198640/1‑149 (MQ=255)
                                                                                         aCCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGt                                                        >  2:3811969/1‑149 (MQ=255)
                                                                                          ccGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTg                                                       >  1:2538554/1‑149 (MQ=255)
                                                                                               ttAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTaa                                                  >  1:2395099/1‑149 (MQ=255)
                                                                                                   aaGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACTCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACAGTACGTCAGGATTGGCGGGATCGGGAAATCGTGCTTAATGcc                                              >  1:1206597/1‑149 (MQ=255)
                                                                                                     gCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCcga                                            <  1:3582864/149‑1 (MQ=255)
                                                                                                      cACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCcgac                                           <  2:2003233/149‑1 (MQ=255)
                                                                                                         tGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGACATCGTGCTTAATGCCGACGAt                                        >  2:1214251/1‑149 (MQ=255)
                                                                                                         tGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGAt                                        >  2:3705504/1‑149 (MQ=255)
                                                                                                          gggAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATg                                       >  1:305162/1‑149 (MQ=255)
                                                                                                            gAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGaa                                     <  1:1702069/149‑1 (MQ=255)
                                                                                                              aTGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAg                                   >  2:3042564/1‑149 (MQ=255)
                                                                                                                gCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCa                                 <  2:1570570/149‑1 (MQ=255)
                                                                                                                 ccAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCAt                                <  2:392128/149‑1 (MQ=255)
                                                                                                                 ccAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCAt                                <  2:653926/149‑1 (MQ=255)
                                                                                                                  cAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCAtt                               >  2:2780243/1‑149 (MQ=255)
                                                                                                                  cAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCAtt                               >  2:3321089/1‑149 (MQ=255)
                                                                                                                    gATAACACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATtct                             <  2:208224/149‑1 (MQ=255)
                                                                                                                     aTAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATtctc                            >  2:2650351/1‑149 (MQ=255)
                                                                                                                     aTAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATtctc                            >  2:3660324/1‑149 (MQ=255)
                                                                                                                      tAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCg                           >  2:2092754/1‑149 (MQ=255)
                                                                                                                        aaaCCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGAt                         <  2:1416788/149‑1 (MQ=255)
                                                                                                                        aaaCCAATCATCTGGAGCGGTTTAAATCGCCAGGCACGCTCGTCAAACGTATCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTAAGGATGGGCGTGATCGGGAAATCGTGCTTAATGCCGACGATGAAACCATTCTCGAt                         >  1:4026685/1‑149 (MQ=255)
                                                                                                                         aaCCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATg                        <  1:2216053/149‑1 (MQ=255)
                                                                                                                          aCCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGc                       <  1:3755329/149‑1 (MQ=255)
                                                                                                                           ccATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGCg                      >  2:3855323/1‑149 (MQ=255)
                                                                                                                            cATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGCgg                     >  2:3276305/1‑149 (MQ=255)
                                                                                                                              ttCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGCGGCa                   <  2:4001014/149‑1 (MQ=255)
                                                                                                                              ttCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGCGGCa                   <  1:3887113/149‑1 (MQ=255)
                                                                                                                                 aTCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGCGGCATTg                <  2:219662/149‑1 (MQ=255)
                                                                                                                                  tCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGCGGCATTgc               <  1:3086976/149‑1 (MQ=255)
                                                                                                                                    tGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGCGGCATTgcgc             <  1:561572/149‑1 (MQ=255)
                                                                                                                                     ggAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGCGGCATTGCGcc            <  2:1314258/149‑1 (MQ=255)
                                                                                                                                     ggAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGCGGCATTGCGcc            <  1:1506987/149‑1 (MQ=255)
                                                                                                                                         cGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGCGGCATTGCGCCAggg        <  2:739827/149‑1 (MQ=255)
                                                                                                                                          ggTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGCGGCATTGCGCCAgggg       <  1:3600820/149‑1 (MQ=255)
                                                                                                                                           gTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGCGGCATTGCGCCAggggg      <  2:2864663/149‑1 (MQ=255)
                                                                                                                                            tttAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGCGGCATTGCGCCAGGGGGc     >  1:2320826/1‑149 (MQ=255)
                                                                                                                                             ttAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGCGGCATTGCGCCAGGGGGCg    <  1:1141123/149‑1 (MQ=255)
                                                                                                                                               aaTTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGCGGCATTGCGCCAGGGGGCGGa  <  2:915281/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
TACAGTCACTTGGGGCCTCGCTCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGCGGCATTGCGCCAGGGGGCGGA  >  CP009273/1451276‑1451567

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