Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,664,138 A→C 58.7% E130D (GAA→GAC hcaF → 3‑phenylpropionate dioxygenase, small (beta) subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,664,1380AC58.7% 212.2 / 225.1 206E130D (GAA→GAChcaF3‑phenylpropionate dioxygenase, small (beta) subunit
Reads supporting (aligned to +/- strand):  ref base A (39/46);  new base C (48/73);  total (87/119)
Fisher's exact test for biased strand distribution p-value = 3.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01

Only 100 of 211 total aligned reads displayed.

GAAACGGGCATGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGCGT  >  CP009273/2663992‑2664285
                                                                                                                                                  |                                                                                                                                                   
gAAACGGGCATGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGacac                                                                                                                                                   >  1:3783284/1‑149 (MQ=255)
  aaCGGGCATGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACAt                                                                                                                                                 <  2:3913327/149‑1 (MQ=255)
      ggCATGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTa                                                                                                                                             >  2:2961877/1‑149 (MQ=255)
         aTGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAGATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGt                                                                                                                                          <  2:3184426/149‑1 (MQ=255)
         aTGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGt                                                                                                                                          >  2:1111139/1‑149 (MQ=255)
         aTGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGt                                                                                                                                          >  2:3163758/1‑149 (MQ=255)
          tGGTCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGtt                                                                                                                                         <  1:492624/149‑1 (MQ=255)
          tGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCACACGGATGTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGtt                                                                                                                                         <  2:3215631/149‑1 (MQ=255)
          tGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGtt                                                                                                                                         <  2:2107944/149‑1 (MQ=255)
          tGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGtt                                                                                                                                         <  1:945556/149‑1 (MQ=255)
          tGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATAGGGCACAAAAAGAGCGCGATGAAACATTCTATGtt                                                                                                                                         >  2:186559/1‑149 (MQ=255)
              cTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGaa                                                                                                                                     >  1:3062326/1‑149 (MQ=255)
              cTGGGCAGAAGAGCCGCCGTCACGCAACAGTCACTTAATCAGCAACTGCCAGATAAGCGACACCGACATAACAAACGTATTTGATGTGCGGGAAAATAAACTGCTTTATCGGGCCCAAAAAGAGCGCGATGACAGCGTCTCTGTTGGaa                                                                                                                                     >  2:261659/1‑149 (MQ=255)
               tGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAAc                                                                                                                                    >  2:2409099/1‑149 (MQ=255)
                gggCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAAcg                                                                                                                                   >  1:1447055/1‑149 (MQ=255)
                gggCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAAcg                                                                                                                                   >  2:156555/1‑149 (MQ=255)
                    agaagaGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGtt                                                                                                                               >  1:1335287/1‑149 (MQ=255)
                      aagaGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTc                                                                                                                             <  1:3571787/149‑1 (MQ=255)
                          gccgccGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACa                                                                                                                         >  1:1712027/1‑149 (MQ=255)
                            cgccgTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACaaa                                                                                                                       <  2:553977/149‑1 (MQ=255)
                             gccgTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAg                                                                                                                      >  1:1477648/1‑149 (MQ=255)
                             gccgTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAg                                                                                                                      >  1:1044997/1‑149 (MQ=255)
                               cgTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGtt                                                                                                                    <  1:786951/149‑1 (MQ=255)
                                gTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTc                                                                                                                   <  2:2168652/149‑1 (MQ=255)
                                 tCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCg                                                                                                                  >  1:560099/1‑149 (MQ=255)
                                  cACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGc                                                                                                                 <  1:2940456/149‑1 (MQ=255)
                                  cACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGc                                                                                                                 <  2:1244594/149‑1 (MQ=255)
                                       aCCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCt                                                                                                            >  1:2564980/1‑149 (MQ=255)
                                             cACTTAATCAGCAACTGCCAGATAATCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGa                                                                                                      >  1:3219987/1‑149 (MQ=255)
                                             cACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGa                                                                                                      >  1:1164420/1‑149 (MQ=255)
                                                       gCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTgg                                                                                            >  1:1764334/1‑149 (MQ=255)
                                                        cAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGgc                                                                                           <  2:45980/149‑1 (MQ=255)
                                                           cTGCCAGATAAGCGAAACCGACATCCCACACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCt                                                                                        <  2:2822310/149‑1 (MQ=255)
                                                           cTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCt                                                                                        <  1:2266399/149‑1 (MQ=255)
                                                           cTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCt                                                                                        <  1:1978925/149‑1 (MQ=255)
                                                            tGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCtt                                                                                       <  2:2868413/149‑1 (MQ=255)
                                                             gCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTg                                                                                      >  2:1367987/1‑149 (MQ=255)
                                                              ccAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGc                                                                                     <  2:1638003/149‑1 (MQ=255)
                                                                  aTAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGa                                                                                 >  2:3147293/1‑149 (MQ=255)
                                                                      gCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACgg                                                                             <  2:807991/149‑1 (MQ=255)
                                                                      gCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACgg                                                                             >  1:2046350/1‑149 (MQ=255)
                                                                      gCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACgg                                                                             <  1:606395/149‑1 (MQ=255)
                                                                        gAAACCGACATCCCACACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGa                                                                           <  1:3144449/149‑1 (MQ=255)
                                                                        gAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGa                                                                           <  1:1857764/149‑1 (MQ=255)
                                                                          aaCCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGata                                                                         <  1:53130/149‑1 (MQ=255)
                                                                           aCCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGatat                                                                        <  2:3770991/149‑1 (MQ=255)
                                                                              gACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGt                                                                     <  2:2639440/149‑1 (MQ=255)
                                                                               aCATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTc                                                                    >  1:1652124/1‑149 (MQ=255)
                                                                               aCATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTc                                                                    <  2:2371845/149‑1 (MQ=255)
                                                                               aCATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTc                                                                    <  1:605208/149‑1 (MQ=255)
                                                                               aCATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTc                                                                    >  2:3066653/1‑149 (MQ=255)
                                                                                cATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTcc                                                                   <  1:4010237/149‑1 (MQ=255)
                                                                                cATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTcc                                                                   <  1:3906595/149‑1 (MQ=255)
                                                                                  tCCCAAACTTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTg                                                                 <  2:3296361/149‑1 (MQ=255)
                                                                                     cAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGAt                                                              <  1:2389458/149‑1 (MQ=255)
                                                                                      aaaCGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATc                                                             <  2:360793/149‑1 (MQ=255)
                                                                                          gTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGc                                                         <  1:683745/149‑1 (MQ=255)
                                                                                          gTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGc                                                         >  2:2990616/1‑149 (MQ=255)
                                                                                          gTATTTGCTGTGCGGGTAAATTAACTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGAACAAGc                                                         <  2:3696903/149‑1 (MQ=255)
                                                                                            aTTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCgg                                                       <  2:4108382/149‑1 (MQ=255)
                                                                                              ttGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAGAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTa                                                     >  2:3359344/1‑149 (MQ=255)
                                                                                               tGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTaa                                                    <  1:1928177/149‑1 (MQ=255)
                                                                                               tGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTaa                                                    <  2:3410862/149‑1 (MQ=255)
                                                                                                 cTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATc                                                  >  1:1421662/1‑149 (MQ=255)
                                                                                                  tgtgCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCa                                                 <  2:3486168/149‑1 (MQ=255)
                                                                                                   gtgCGGTTAATTTATCTTCTTTATCGGTCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGAGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCAc                                                <  2:2609333/149‑1 (MQ=255)
                                                                                                   gtgCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCAc                                                <  2:2609379/149‑1 (MQ=255)
                                                                                                         gTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCa                                          >  1:2151731/1‑149 (MQ=255)
                                                                                                          tAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCAt                                         <  2:1193014/149‑1 (MQ=255)
                                                                                                            aaTTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATaa                                       <  1:2680523/149‑1 (MQ=255)
                                                                                                            aaTTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATaa                                       <  2:154212/149‑1 (MQ=255)
                                                                                                                aTCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTg                                   >  2:2015765/1‑149 (MQ=255)
                                                                                                                  cTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAg                                 <  2:560099/149‑1 (MQ=255)
                                                                                                                    gCTTTATCGGGCACAAAACGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAgtg                               <  2:3973127/149‑1 (MQ=255)
                                                                                                                    gCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAgtg                               <  2:2112934/149‑1 (MQ=255)
                                                                                                                     cTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAgtgt                              <  1:2628981/149‑1 (MQ=255)
                                                                                                                      tttATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTa                             <  2:1907064/149‑1 (MQ=255)
                                                                                                                        tataGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACt                           >  2:3439797/1‑149 (MQ=255)
                                                                                                                         atagggCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTg                          <  1:2711982/146‑1 (MQ=255)
                                                                                                                         aTCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTg                          >  2:2378667/1‑149 (MQ=255)
                                                                                                                          tCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGt                         >  1:705618/1‑149 (MQ=255)
                                                                                                                             ggCACAAACAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCt                      <  2:1973138/149‑1 (MQ=255)
                                                                                                                                 caAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCtgatg                  >  2:2623496/1‑149 (MQ=255)
                                                                                                                                 caAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCtgatg                  >  1:4052756/1‑149 (MQ=255)
                                                                                                                                  aaaaaGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCtgatga                 <  1:2961877/149‑1 (MQ=255)
                                                                                                                                    aaaGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAAt               <  2:1426191/149‑1 (MQ=255)
                                                                                                                                     aaGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATc              >  1:3046421/1‑149 (MQ=255)
                                                                                                                                      agagCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCg             >  1:138356/1‑149 (MQ=255)
                                                                                                                                      agagCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCg             >  1:2212956/1‑149 (MQ=255)
                                                                                                                                       gagCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGa            >  1:319074/1‑149 (MQ=255)
                                                                                                                                         gcgcGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAAt          <  2:2390761/149‑1 (MQ=255)
                                                                                                                                          cgcgATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAAtt         <  2:1501452/149‑1 (MQ=255)
                                                                                                                                          cgcgATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAAtt         <  2:4065421/149‑1 (MQ=255)
                                                                                                                                           gcgATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAAttt        >  2:2230095/1‑149 (MQ=255)
                                                                                                                                             gATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAGTTTAt      <  1:3359344/149‑1 (MQ=255)
                                                                                                                                             gATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTAt      <  1:1367987/149‑1 (MQ=255)
                                                                                                                                              aTGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATg     >  1:581024/1‑149 (MQ=255)
                                                                                                                                               tGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGc    <  2:3635965/149‑1 (MQ=255)
                                                                                                                                               tGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGc    >  2:1692778/1‑149 (MQ=255)
                                                                                                                                                 aaaCATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGCgt  <  2:4051624/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GAAACGGGCATGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGCGT  >  CP009273/2663992‑2664285

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: