Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,451,422 A→T 43.6% 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,4220AT43.6% 153.0 / inf 314T246S (ACG→TCG) paaEring 1,2‑phenylacetyl‑CoA epoxidase, NAD(P)H oxidoreductase component
Reads supporting (aligned to +/- strand):  ref base A (95/82);  new base T (76/61);  total (171/143)
Fisher's exact test for biased strand distribution p-value = 8.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

Only 100 of 315 total aligned reads displayed.

AATTACAGTCACTTGGGGCCTCGCTCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGCGGCATTGCGCCAGGGGGCGGATCT  >  CP009273/1451273‑1451570
                                                                                                                                                     |                                                                                                                                                    
  ttACAGTCACTTGGGGCCTCGCTCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCg                                                                                                                                                    >  1:2666297/1‑150 (MQ=255)
  ttACAGTCACTTGGGGCCTCGCTCATTAATTTTCGTCTGTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCg                                                                                                                                                    <  2:3944252/150‑1 (MQ=255)
    aCAGTCACTTGGGGCCTCGCTCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGcc                                                                                                                                                  <  2:2095538/150‑1 (MQ=255)
    aCAGTCACTTGGGGCCTCGCTCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGcc                                                                                                                                                  <  2:3017074/150‑1 (MQ=255)
      aGTCACTTGGGGCCTCGCTCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTg                                                                                                                                                <  2:2510734/150‑1 (MQ=255)
        tCACTTGGGGCCTCGCTCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGc                                                                                                                                              <  2:1328398/150‑1 (MQ=255)
         cACTTGGGGCCTCGCTCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCa                                                                                                                                             <  1:2164028/150‑1 (MQ=255)
          aCTTGGGGCCTCGCTCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCAc                                                                                                                                            <  1:4067493/150‑1 (MQ=255)
              ggggCCTCGCTCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCAcgcgc                                                                                                                                        <  2:356159/150‑1 (MQ=255)
                ggCCTCGCTCATTAATTTTCGTCTTTATGATGAGTCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGt                                                                                                                                      <  1:755715/150‑1 (MQ=255)
                ggCCTCGCTCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGt                                                                                                                                      >  1:2280672/1‑150 (MQ=255)
                ggCCTCGCTCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGt                                                                                                                                      >  1:2566642/1‑150 (MQ=255)
                ggCCTCGCTCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGt                                                                                                                                      <  2:1524609/150‑1 (MQ=255)
                 gCCTCGCTCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTc                                                                                                                                     >  2:3690876/1‑150 (MQ=255)
                 gCCTCGCTCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTc                                                                                                                                     <  2:4743675/150‑1 (MQ=255)
                 gCCTCGCTCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTc                                                                                                                                     <  1:1203221/150‑1 (MQ=255)
                  ccTCGCTCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCa                                                                                                                                    >  2:3231571/1‑150 (MQ=255)
                  ccTCGCTCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCa                                                                                                                                    <  1:2867913/150‑1 (MQ=255)
                   cTCGCTCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCaa                                                                                                                                   <  2:1152071/150‑1 (MQ=255)
                    tCGCTCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCaaa                                                                                                                                  <  1:109105/150‑1 (MQ=255)
                    tCGCTCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCaaa                                                                                                                                  >  1:906193/1‑150 (MQ=255)
                    tCGCTCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCaaa                                                                                                                                  >  1:1689126/1‑150 (MQ=255)
                        tCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTa                                                                                                                              <  1:480078/150‑1 (MQ=255)
                         cATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAg                                                                                                                             >  2:1488139/1‑150 (MQ=255)
                          aTTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGc                                                                                                                            >  2:34421/1‑150 (MQ=255)
                            tAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGt                                                                                                                          >  1:3245417/1‑150 (MQ=255)
                            gtaTTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGt                                                                                                                          <  2:2916966/148‑1 (MQ=255)
                               ttttCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTaa                                                                                                                       <  2:2323941/150‑1 (MQ=255)
                                tttCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAAc                                                                                                                      >  2:1624415/1‑150 (MQ=255)
                                tttCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAAc                                                                                                                      >  2:1569837/1‑150 (MQ=255)
                                tttCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAAc                                                                                                                      >  1:3668473/1‑150 (MQ=255)
                                  tCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGt                                                                                                                    <  1:378748/150‑1 (MQ=255)
                                       tttATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCaaa                                                                                                               <  1:2905198/150‑1 (MQ=255)
                                          atgatgAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTg                                                                                                            >  2:794164/1‑150 (MQ=255)
                                            gatgaGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGAc                                                                                                          <  2:241859/150‑1 (MQ=255)
                                             atgaGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACg                                                                                                         <  2:120272/150‑1 (MQ=255)
                                             atgaGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACg                                                                                                         >  2:1275209/1‑150 (MQ=255)
                                               gaGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGa                                                                                                       <  2:5053971/150‑1 (MQ=255)
                                                aGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGAc                                                                                                      <  2:531556/150‑1 (MQ=255)
                                                aGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGAc                                                                                                      >  1:93504/1‑150 (MQ=255)
                                                aGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGAc                                                                                                      <  1:4733088/150‑1 (MQ=255)
                                                 ggCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGAcc                                                                                                     >  2:3316637/1‑149 (MQ=255)
                                                   cATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACaaa                                                                                                   <  2:398644/150‑1 (MQ=255)
                                                   cATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACaaa                                                                                                   <  2:719455/150‑1 (MQ=255)
                                                     tttatttGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACaaaaa                                                                                                 >  2:4818826/1‑150 (MQ=255)
                                                       tatttGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGt                                                                                               >  2:1911020/1‑150 (MQ=255)
                                                       tatttGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGt                                                                                               >  1:628771/1‑150 (MQ=255)
                                                       tatttGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGt                                                                                               >  2:1040760/1‑150 (MQ=255)
                                                         tttGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGa                                                                                             >  1:2423198/1‑150 (MQ=255)
                                                          ttGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGAc                                                                                            <  1:424026/150‑1 (MQ=255)
                                                           tgtgGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACt                                                                                           >  1:4465431/1‑150 (MQ=255)
                                                           tgtgGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACt                                                                                           <  1:4638082/150‑1 (MQ=255)
                                                             tgGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGt                                                                                         <  1:1389935/150‑1 (MQ=255)
                                                                 ccGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAAGGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGt                                                                                     <  2:4851427/150‑1 (MQ=255)
                                                                   ggcggcGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCa                                                                                   <  2:772128/150‑1 (MQ=255)
                                                                     cggcgATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAgg                                                                                 >  1:2055483/1‑150 (MQ=255)
                                                                      ggcgATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAATGTGACTGTACGTCAGGa                                                                                >  2:3829461/1‑150 (MQ=255)
                                                                        cgATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATg                                                                              <  1:3836464/150‑1 (MQ=255)
                                                                         gatgatGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATgg                                                                             >  2:4610666/1‑150 (MQ=255)
                                                                           tgatgGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGc                                                                           <  2:1959819/150‑1 (MQ=255)
                                                                             atggatgACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCgg                                                                         >  2:2041655/1‑150 (MQ=255)
                                                                             atggatgACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCgg                                                                         >  1:3114745/1‑150 (MQ=255)
                                                                                   gACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTCGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCgg                                                                   >  1:4623582/1‑150 (MQ=255)
                                                                                   gACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCgg                                                                   <  2:4209825/150‑1 (MQ=255)
                                                                                      gcgGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGaa                                                                >  1:1537906/1‑150 (MQ=255)
                                                                                        ggAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAAt                                                              >  1:1015243/1‑150 (MQ=255)
                                                                                         gAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATc                                                             >  1:4887116/1‑150 (MQ=255)
                                                                                               gCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCtt                                                       <  1:5116371/150‑1 (MQ=255)
                                                                                                ccTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTa                                                      <  2:4077651/150‑1 (MQ=255)
                                                                                                 cTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTaa                                                     >  2:992180/1‑150 (MQ=255)
                                                                                                  ttAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAAt                                                    >  1:4308779/1‑150 (MQ=255)
                                                                                                   tAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATg                                                   >  1:4545334/1‑150 (MQ=255)
                                                                                                   tAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATg                                                   <  1:440700/150‑1 (MQ=255)
                                                                                                     aaaGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGcc                                                 >  2:2599596/1‑150 (MQ=255)
                                                                                                      aaGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCcg                                                >  2:2982518/1‑150 (MQ=255)
                                                                                                        gCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCcgac                                              <  2:3361589/150‑1 (MQ=255)
                                                                                                          aCTCGGAAACCCAGATAAAACAATTCATCAGGAGAGGTATAATACGCCACGCACGCGCGACAAACGTAGAGTAAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGCGAAATCGTGCTTAATGCagacga                                            >  2:2628364/1‑150 (MQ=255)
                                                                                                             gggAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGa                                         <  2:93118/150‑1 (MQ=255)
                                                                                                             gggAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAAGTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGa                                         >  2:2673550/1‑150 (MQ=255)
                                                                                                              ggATTGCCTGATAAAACCATTCATCTGGAGGGGTTTAATACGCCTGGCACGCGCGTCAAACGTTGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGaa                                        <  2:2856598/150‑1 (MQ=255)
                                                                                                                aaTGCCAGATAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAg                                      >  1:2155763/1‑150 (MQ=255)
                                                                                                                aaTGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAg                                      <  1:3334192/150‑1 (MQ=255)
                                                                                                                    ccAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCAtt                                  <  2:4383256/150‑1 (MQ=255)
                                                                                                                        aTAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCg                              >  1:1216542/1‑150 (MQ=255)
                                                                                                                         tAAAACCATTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGa                             <  2:2624359/150‑1 (MQ=255)
                                                                                                                           aaaCCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATg                           <  1:4126050/150‑1 (MQ=255)
                                                                                                                           aaaCCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATg                           <  1:3461995/150‑1 (MQ=255)
                                                                                                                                aTTCATCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGCGGCa                      >  2:3981283/1‑150 (MQ=255)
                                                                                                                                aTTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGCGGCa                      <  1:1857375/150‑1 (MQ=255)
                                                                                                                                aTTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGCGGCa                      <  1:2931677/150‑1 (MQ=255)
                                                                                                                                   cATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGCGGCATTg                   >  1:645393/1‑150 (MQ=255)
                                                                                                                                    aTCTGGAGCGGTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGCGGCATTgc                  >  1:3966250/1‑150 (MQ=255)
                                                                                                                                    aTCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGCGGCATTgc                  >  2:4935249/1‑150 (MQ=255)
                                                                                                                                    aTCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGCGGCATTgc                  <  2:3792139/150‑1 (MQ=255)
                                                                                                                                       tGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGCGGCATTGCGcc               >  2:3300159/1‑150 (MQ=255)
                                                                                                                                             ggTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGCGGCATTGCGCCAggggg         >  2:2820501/1‑150 (MQ=255)
                                                                                                                                              gTTTAATTCGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGCGGCATTGCGCCAGGGGGc        <  1:4824890/150‑1 (MQ=255)
                                                                                                                                               tttAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGCGGCATTGCGCCAGGGGGCg       <  1:4999516/150‑1 (MQ=255)
                                                                                                                                                ttAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGCGGCATTGCGCCAGGGGGCgg      >  1:1831493/1‑150 (MQ=255)
                                                                                                                                                   aTACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGCGGCATTGCGCCAGGGGGCGGATc   >  2:3130/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AATTACAGTCACTTGGGGCCTCGCTCATTAATTTTCGTCTTTATGATGAGGCATTTATTTGTGGTCCGGCGGCGATGATGGATGACGCGGAAACCGCCTTAAAAGCACTGGGAATGCCAGATAAAACCATTCATCTGGAGCGGTTTAATACGCCTGGCACGCGCGTCAAACGTAGCGTTAACGTGCAAAGTGACGGACAAAAAGTGACTGTACGTCAGGATGGGCGGGATCGGGAAATCGTGCTTAATGCCGACGATGAAAGCATTCTCGATGCGGCATTGCGCCAGGGGGCGGATCT  >  CP009273/1451273‑1451570

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: