Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,169,196 T→A 16.2% E57V (GAA→GTA)  yqiA ← acyl CoA esterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,169,1960TA16.2% 559.7 / 85.9 208E57V (GAA→GTA) yqiAacyl CoA esterase
Reads supporting (aligned to +/- strand):  ref base T (88/86);  new base A (21/13);  total (109/99)
Fisher's exact test for biased strand distribution p-value = 2.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 212 total aligned reads displayed.

GTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCTTTTGCAGAGCGCGGCGAGCTGTTGAAACCGT  >  CP009273/3169048‑3169343
                                                                                                                                                    |                                                                                                                                                   
gTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGtt                                                                                                                                                     <  2:272685/149‑1 (MQ=255)
 ttCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTc                                                                                                                                                    >  1:1943811/1‑149 (MQ=255)
    tGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACAGCCATGTTCCAg                                                                                                                                                 <  2:891240/149‑1 (MQ=255)
     gACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTACCAgg                                                                                                                                                >  2:1402821/1‑149 (MQ=255)
      aCCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGa                                                                                                                                               >  1:2408488/1‑149 (MQ=255)
        cGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACa                                                                                                                                             >  1:2590325/1‑149 (MQ=255)
          agaTAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTACCAGGACAAt                                                                                                                                           >  1:2110005/1‑149 (MQ=255)
          agaTAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACAGCCATGTTCCAGGACAAt                                                                                                                                           >  2:3522518/1‑149 (MQ=255)
          agaTAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACAGCCATGTTCCAGGACAAt                                                                                                                                           >  2:3208631/1‑149 (MQ=255)
          agaTAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACAGCCATGTTCCAGGACAAt                                                                                                                                           >  2:2080952/1‑149 (MQ=255)
           gaTAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATg                                                                                                                                          >  2:3170371/1‑149 (MQ=255)
           gaTAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTACCAGGACAATg                                                                                                                                          >  1:1243085/1‑149 (MQ=255)
           gaTAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTACCAGGACAATg                                                                                                                                          >  1:3721092/1‑149 (MQ=255)
           gaTAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACAGCCATGTTCCAGGACAATg                                                                                                                                          >  2:3318155/1‑149 (MQ=255)
             tAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACAGCCATGTTCCAGGACAATGGa                                                                                                                                        <  1:1921629/149‑1 (MQ=255)
                 ccGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTcc                                                                                                                                    <  2:1078729/149‑1 (MQ=255)
                  cGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACAGCCATGTTCCAGGACAATGGATTCca                                                                                                                                   <  2:3078101/149‑1 (MQ=255)
                    tCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCcagc                                                                                                                                 >  1:2271828/1‑149 (MQ=255)
                            tCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTACCAGGACAATGGATTCCAGCAGCTCTGc                                                                                                                         <  1:3794576/149‑1 (MQ=255)
                            tCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACAGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGc                                                                                                                         <  2:4007460/149‑1 (MQ=255)
                                  gggCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACCATGGATTCCAGCAGCTCTGCCGCGTc                                                                                                                   >  2:757183/1‑149 (MQ=255)
                                      gcACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACAGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGaa                                                                                                               >  2:1436180/1‑149 (MQ=255)
                                        aCCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACAGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAgg                                                                                                             <  2:761479/149‑1 (MQ=255)
                                        aCCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACAGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAgg                                                                                                             <  2:2392168/149‑1 (MQ=255)
                                           gccgGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATa                                                                                                          >  1:3488332/1‑149 (MQ=255)
                                             cgGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACAGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATAcg                                                                                                        <  2:3928437/149‑1 (MQ=255)
                                              gggTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACAGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATAcgg                                                                                                       >  1:2062608/1‑149 (MQ=255)
                                              gggTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACAGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATAcgg                                                                                                       >  1:320836/1‑149 (MQ=255)
                                               ggTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTACCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATActgc                                                                                                      >  1:2872080/1‑149 (MQ=255)
                                                gTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACAGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATAcggcg                                                                                                     >  1:2653982/1‑149 (MQ=255)
                                                  ttACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACggcggc                                                                                                   <  1:1332183/149‑1 (MQ=255)
                                                     ccaccaCTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAAc                                                                                                >  1:3968617/1‑149 (MQ=255)
                                                     ccaccaCTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACAGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAAc                                                                                                <  2:3152690/149‑1 (MQ=255)
                                                     ccaccaCTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACAGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAAc                                                                                                <  2:3336436/149‑1 (MQ=255)
                                                       accacTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATCTCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACAGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTg                                                                                              >  1:1440099/1‑149 (MQ=255)
                                                         cacTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTACCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCg                                                                                            >  1:3123815/1‑149 (MQ=255)
                                                         cacTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTACCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCg                                                                                            >  1:237459/1‑149 (MQ=255)
                                                            tGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACAGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGaa                                                                                         >  1:3059205/1‑149 (MQ=255)
                                                               gggCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACAGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGACTGa                                                                                      >  1:39177/1‑149 (MQ=255)
                                                               gggCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACAGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGa                                                                                      >  2:1687509/1‑149 (MQ=255)
                                                                ggCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGat                                                                                     <  1:1897583/149‑1 (MQ=255)
                                                                   agcaTAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACAGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGatcat                                                                                  <  1:3610780/149‑1 (MQ=255)
                                                                        aaaaCATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACAGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCaa                                                                             <  1:1971248/149‑1 (MQ=255)
                                                                           aCATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGt                                                                          >  1:1596017/1‑149 (MQ=255)
                                                                           aCATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGt                                                                          >  1:3712008/1‑149 (MQ=255)
                                                                           aCATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGt                                                                          >  1:1596081/1‑149 (MQ=255)
                                                                            cATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACAGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTc                                                                         <  1:423568/149‑1 (MQ=255)
                                                                              ttGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAg                                                                       >  2:1580648/1‑149 (MQ=255)
                                                                               tgtgACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTACCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAgg                                                                      <  1:854645/149‑1 (MQ=255)
                                                                                  gACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACAGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTg                                                                   >  1:2201746/1‑149 (MQ=255)
                                                                                    cAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGAt                                                                 <  2:3841431/149‑1 (MQ=255)
                                                                                     aaCCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACAGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGACTGATCATCTCAACGTCAGGGTGATg                                                                >  2:2115318/1‑149 (MQ=255)
                                                                                      aCCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACAGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGt                                                               >  2:3783775/1‑149 (MQ=255)
                                                                                      aCCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACAGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGt                                                               <  1:3967692/149‑1 (MQ=255)
                                                                                         aGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACAGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTcc                                                            <  1:4033813/149‑1 (MQ=255)
                                                                                          ggTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACAGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGACTGATCATCTCAACGTCAGGGTGATGTTCCg                                                           >  2:4061364/1‑149 (MQ=255)
                                                                                           gTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACAGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCgc                                                          >  1:3250041/1‑149 (MQ=255)
                                                                                            tGGCGTAATATCCCCCCCGTGACGAACCGACAATACCCAGCGAATCACAGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCgcc                                                         <  2:652281/149‑1 (MQ=255)
                                                                                             ggCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCgcca                                                        <  2:1254942/149‑1 (MQ=255)
                                                                                                gTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACAGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCgccagcc                                                     <  2:3800870/149‑1 (MQ=255)
                                                                                                 tAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCgccagcca                                                    >  2:113934/1‑149 (MQ=255)
                                                                                                   atatCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACAGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGt                                                  <  1:412944/149‑1 (MQ=255)
                                                                                                    tatCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGtt                                                 >  2:4079608/1‑149 (MQ=255)
                                                                                                     atCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACAGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGttt                                                >  1:3409249/1‑149 (MQ=255)
                                                                                                     atCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACAGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGttt                                                <  1:1910923/149‑1 (MQ=255)
                                                                                                     atCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACAGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGttt                                                <  1:2573240/149‑1 (MQ=255)
                                                                                                       ccccccAGTGACGAACCGACAATACCCAGCGAATCACAGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGttttt                                              <  1:1000974/149‑1 (MQ=255)
                                                                                                         ccccAGTGACGAACCGACAATACCCAGCGAATCACAGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTa                                            <  1:977601/149‑1 (MQ=255)
                                                                                                              gTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAg                                       >  1:911408/1‑149 (MQ=255)
                                                                                                               tGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGc                                      <  2:2948716/149‑1 (MQ=255)
                                                                                                               tGACGAACCGACAATACCCAGCGAATCACCGCCATGTACCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGc                                      <  2:470695/149‑1 (MQ=255)
                                                                                                               tGACGAACCGACAATACCCAGCGAATCACAGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGc                                      >  1:3255001/1‑149 (MQ=255)
                                                                                                                    aaCCGACAATACCCAGCGACTCACAGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCt                                 >  2:3657608/1‑149 (MQ=255)
                                                                                                                     aCCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCtt                                <  2:3883756/149‑1 (MQ=255)
                                                                                                                          cAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCTTTTGCa                           <  1:365679/149‑1 (MQ=255)
                                                                                                                           aaTACCCAGCGACTCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCTTTTGCag                          <  2:3648806/149‑1 (MQ=255)
                                                                                                                           aaTACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCTTTTGCag                          >  2:1441357/1‑149 (MQ=255)
                                                                                                                            aTACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCTTTTGCaga                         >  1:2728927/1‑149 (MQ=255)
                                                                                                                             tACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCTTTTGCagag                        >  1:2854780/1‑149 (MQ=255)
                                                                                                                                ccAGCGAATCACAGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCTTTTGCAGAgcgc                     >  2:2905834/1‑149 (MQ=255)
                                                                                                                                ccAGCGAATCACAGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCTTTTGCAGAgcgc                     <  2:3066224/149‑1 (MQ=255)
                                                                                                                                 cAGCGAATCACAGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCTTTTGCAGAGcgcg                    <  2:1498369/149‑1 (MQ=255)
                                                                                                                                 cAGCGAATCACAGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCGTTTGCAGAGcgcg                    >  2:2883909/1‑149 (MQ=255)
                                                                                                                                     taaTCACAGCCATGTTCCAGGACAATGGAGTCCAGCAGCTCTGCGGCGTCGGAAGGATACGGAGGCAACTGCGGAATTATCATCTCAACGTCAGGGTGATGTTCCGCCTGAAAGTTTTTTAACAAGCTCGCTTTTGCAGAGCGCTGCGa                <  2:3636918/148‑1 (MQ=255)
                                                                                                                                      aaTCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCTTTTGCAGAGCGCGGCGAg               >  1:3007400/1‑149 (MQ=255)
                                                                                                                                      aaTCACCGCCATGTACCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCTTTTGCAGAGCGCGGCGAg               >  1:1383533/1‑149 (MQ=255)
                                                                                                                                      aaTCACAGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCTTTTGCAGAGCGCGGCGAg               >  2:1731192/1‑149 (MQ=255)
                                                                                                                                       aTCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCTTTTGCAGAGCGCGGCGAGc              >  2:3378407/1‑149 (MQ=255)
                                                                                                                                       aTCACAGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCTTTTGCAGAGCGCGGCGAGc              <  1:1307231/149‑1 (MQ=255)
                                                                                                                                           ccGCCATGTACCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCTTTTGCAGAGCGCGGCGAGCTGtt          >  1:4050212/1‑149 (MQ=255)
                                                                                                                                            cGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCTTTTGCAGAGCGCGGCGAGCTGTTg         >  2:2984984/1‑149 (MQ=255)
                                                                                                                                             gCCATGTACCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCTTTTGCAGAGCGCGGCGAGCTGTTGa        >  1:879242/1‑149 (MQ=255)
                                                                                                                                               cATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCTTTTGCAGAGCGCGGCGAGCTGTTGaaa      <  2:1829770/149‑1 (MQ=255)
                                                                                                                                               cATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCTTTTGCAGAGCGCGGCGAGCTGTTGaaa      <  1:1436530/149‑1 (MQ=255)
                                                                                                                                               cATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCTTTTGCAGAGCGCGGCGAGCTGTTGaaa      <  2:1568091/149‑1 (MQ=255)
                                                                                                                                               cATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCTTTTGCAGAGCGCGGCGAGCTGTTGaaa      <  1:573343/149‑1 (MQ=255)
                                                                                                                                               cATGTACCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCTTTTGCAGAGCGCGGCGAGCTGTTGaaa      <  1:2101076/149‑1 (MQ=255)
                                                                                                                                                 tGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCTTTTGCAGAGCGCGGCGAGCTGTTGAAAcc    >  2:185108/1‑149 (MQ=255)
                                                                                                                                                 tGTACCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCTTTTGCAGAGCGCGGCGAGCTGTTGAAAcc    <  2:1962085/149‑1 (MQ=255)
                                                                                                                                                   ttCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCTTTTGCAGAGCGCGGCGAGCTGTTGAAACCgt  >  1:4044399/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GTTCTGACCGAGATAGTCCGTCAGCAGTTCAAACGGGCGCACCGCCGGGTTTACCACCACTGCGGGCAGCATAAAACATTGTGACAACCAGGTGGCGTAATATCCCCCCAGTGACGAACCGACAATACCCAGCGAATCACCGCCATGTTCCAGGACAATGGATTCCAGCAGCTCTGCCGCGTCGGAAGGATACGGCGGCAACTGCGGAATGATCATCTCAACGTCAGGGTGATGTTCCGCCAGCCAGTTTTTTAACAAGCTCGCTTTTGCAGAGCGCGGCGAGCTGTTGAAACCGT  >  CP009273/3169048‑3169343

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: