Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,664,115 C→A 6.1% R123R (CGG→AGG)  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,1150CA 6.1% 595.3 / 11.2 157R123R (CGG→AGG) hcaF3‑phenylpropionate dioxygenase, small (beta) subunit
Reads supporting (aligned to +/- strand):  ref base C (72/75);  new base A (5/5);  total (77/80)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 158 total aligned reads displayed.

CTGGAGCGGCGAATCGCCCGTCTGGAAACGGGCATGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTG  >  CP009273/2663968‑2664261
                                                                                                                                                   |                                                                                                                                                  
cTGGAGCGGCGAATCGCCCGTCTGGAAACGGGCATGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATc                                                                                                                                                    <  1:1806089/148‑1 (MQ=255)
    aGCGGCGAATCGCCCGTCTGGAAACGGGCATGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGc                                                                                                                                                >  1:2180513/1‑148 (MQ=255)
      cggcgAATCGCCCGTCTGGAAACGGGCATGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGcac                                                                                                                                              >  2:1204888/1‑148 (MQ=255)
        gcgAATCGCCCGTCTGGAAACGGGCATGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACaa                                                                                                                                            >  2:2826265/1‑148 (MQ=255)
          gAATCGCCCGTCTGGAAACGGGCATGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACaaaa                                                                                                                                          <  1:153592/148‑1 (MQ=255)
          gAATCGCCCGTCTGGAAACGGGCATGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACaaaa                                                                                                                                          <  1:153628/148‑1 (MQ=255)
            aTCGCCCGTCTGGAAACGGGCATGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAag                                                                                                                                        >  2:595044/1‑148 (MQ=255)
            aTCGCCCGTCTGGAAACGGGCATGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAag                                                                                                                                        >  2:915355/1‑148 (MQ=255)
                cccGTCTGGAAACGGGCATGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAgcg                                                                                                                                    >  1:1217305/1‑148 (MQ=255)
                  cGTCTGGAAACGGGCATGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGcgcg                                                                                                                                  <  1:2278290/148‑1 (MQ=255)
                   gTCTGGAAACGGGCATGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGa                                                                                                                                 <  2:2107150/148‑1 (MQ=255)
                   gTCTGGAAACGGGCATGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGa                                                                                                                                 >  1:2664382/1‑148 (MQ=255)
                    tCTGGACACGGGCAGGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCACACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGAt                                                                                                                                <  2:2527817/148‑1 (MQ=255)
                        gAAACGGGCATGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGaaa                                                                                                                            <  1:2876566/148‑1 (MQ=255)
                         aaaCGGTCATGGCCTGGGCAGAAGAGCCGCCGTCACGCACACGTCAATTAATCAGAAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGAGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGacac                                                                                                                           <  1:45694/148‑1 (MQ=255)
                         aaaCGGGCATGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAAc                                                                                                                           <  2:2203629/148‑1 (MQ=255)
                          aaCGGGCATGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACa                                                                                                                          <  2:1637908/148‑1 (MQ=255)
                          aaCGGGCATGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACa                                                                                                                          >  2:3217254/1‑148 (MQ=255)
                           aCGGGCATGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACAt                                                                                                                         >  1:2931145/1‑148 (MQ=255)
                             gggCATGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTc                                                                                                                       <  2:2581265/148‑1 (MQ=255)
                                   ggCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGtt                                                                                                                 <  2:266773/148‑1 (MQ=255)
                                      cTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAACAGAGCGCGATGACACATTCTCTGTTGGa                                                                                                              >  2:1518136/1‑148 (MQ=255)
                                      cTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGa                                                                                                              >  2:1518135/1‑148 (MQ=255)
                                      cTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGa                                                                                                              <  1:1411570/148‑1 (MQ=255)
                                      cTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGa                                                                                                              <  1:2427040/148‑1 (MQ=255)
                                        gggCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAAc                                                                                                            >  1:2106095/1‑148 (MQ=255)
                                         ggCAGAAGAGCCGCCGTCACGCACTCGTCACTTAATCAGCAACTGCCAGGTAAGCGAAACCGACATCCCAAATGTATTTGCCGTGAGGGTGAATTATCTACTTTATCGGGCACAAAAAGAGCGCGATGAAACATTATATGTTGGAAcg                                                                                                           >  2:2940852/1‑148 (MQ=255)
                                           cAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAAcgcg                                                                                                         <  1:304114/148‑1 (MQ=255)
                                             gaagaGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTCTGTTGGAACGCGtt                                                                                                       >  2:655483/1‑148 (MQ=255)
                                              aagaGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGttt                                                                                                      >  1:2218116/1‑148 (MQ=255)
                                              aagaGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGttt                                                                                                      <  2:705540/148‑1 (MQ=255)
                                              aagaGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGttt                                                                                                      <  1:2476396/148‑1 (MQ=255)
                                                       cgTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGt                                                                                             <  2:1263427/148‑1 (MQ=255)
                                                       cgTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGt                                                                                             <  2:355605/148‑1 (MQ=255)
                                                          cACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCg                                                                                          >  1:2354562/1‑148 (MQ=255)
                                                          cACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCg                                                                                          >  1:3118415/1‑148 (MQ=255)
                                                          cACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCg                                                                                          <  2:769558/148‑1 (MQ=255)
                                                            cGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGcc                                                                                        <  2:2136413/148‑1 (MQ=255)
                                                            cGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGcc                                                                                        <  2:2136566/148‑1 (MQ=255)
                                                            cGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGcc                                                                                        >  1:1963986/1‑148 (MQ=255)
                                                               aCCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTc                                                                                     <  1:1036108/148‑1 (MQ=255)
                                                                cccGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCt                                                                                    <  2:2729395/148‑1 (MQ=255)
                                                                cccGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCt                                                                                    <  2:2933717/148‑1 (MQ=255)
                                                                 ccGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTg                                                                                   >  2:390041/1‑148 (MQ=255)
                                                                 ccGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTg                                                                                   >  1:1381708/1‑148 (MQ=255)
                                                                  cGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTgg                                                                                  <  2:2574335/148‑1 (MQ=255)
                                                                     cACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAg                                                                               <  2:3170431/148‑1 (MQ=255)
                                                                        ttAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATg                                                                            <  1:2671001/148‑1 (MQ=255)
                                                                         tAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGa                                                                           <  1:2145883/148‑1 (MQ=255)
                                                                         tAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGa                                                                           >  2:2632290/1‑148 (MQ=255)
                                                                         tAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGa                                                                           <  2:920414/148‑1 (MQ=255)
                                                                             cAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAAc                                                                       >  1:1426265/1‑148 (MQ=255)
                                                                             cAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAAc                                                                       <  1:759884/148‑1 (MQ=255)
                                                                                cAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTgg                                                                    >  1:2557316/1‑148 (MQ=255)
                                                                                cAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGATTCGACAAAGTTCGCCGTCAGGAAGATGACAACAgg                                                                    >  2:2211972/1‑148 (MQ=255)
                                                                                   cTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAGCGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGgcgc                                                                 >  2:8920/1‑148 (MQ=255)
                                                                                    tGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCt                                                                >  2:2471700/1‑148 (MQ=255)
                                                                                      ccAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACACCAGGCGCTTg                                                              >  2:3162868/1‑148 (MQ=255)
                                                                                       cAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGc                                                             <  2:186496/148‑1 (MQ=255)
                                                                                          aTAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTgg                                                          >  1:2876455/1‑148 (MQ=255)
                                                                                            aaGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATAGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGaa                                                        >  2:714804/1‑148 (MQ=255)
                                                                                              gCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACg                                                      >  2:1763109/1‑148 (MQ=255)
                                                                                               cGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACgg                                                     >  1:1032717/1‑148 (MQ=255)
                                                                                               cGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACgg                                                     <  1:862256/148‑1 (MQ=255)
                                                                                                 aaaCCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGa                                                   >  2:7949/1‑148 (MQ=255)
                                                                                                  aaCCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGat                                                  <  2:1217305/148‑1 (MQ=255)
                                                                                                  aaCCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGat                                                  <  2:1217393/148‑1 (MQ=255)
                                                                                                   aCCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGata                                                 >  1:1662148/1‑148 (MQ=255)
                                                                                                       aCATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGt                                             <  1:3255960/148‑1 (MQ=255)
                                                                                                        cATCCCACCCGTATTTGCTGTGCGGGTAACTTATCTGCGTTATCGGGCACCAAAACAGCGCGATGACTCATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTc                                            <  2:15746/148‑1 (MQ=255)
                                                                                                        cATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACCAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTc                                            <  2:1692858/148‑1 (MQ=255)
                                                                                                        cATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTc                                            <  2:15759/148‑1 (MQ=255)
                                                                                                        cATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATAGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTc                                            <  2:2254950/148‑1 (MQ=255)
                                                                                                            ccAAATGTATTTGCCGTGAGGGTGAATTATCTACTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTgg                                        <  2:2064619/148‑1 (MQ=255)
                                                                                                             cAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGa                                       <  1:3240649/148‑1 (MQ=255)
                                                                                                              aaaCGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGAt                                      <  1:2471700/148‑1 (MQ=255)
                                                                                                               aaCGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATc                                     >  2:3313123/1‑148 (MQ=255)
                                                                                                                    aTTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCg                                >  1:2633487/1‑148 (MQ=255)
                                                                                                                    aTTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCg                                <  1:66923/148‑1 (MQ=255)
                                                                                                                      ttGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGt                              >  2:2865136/1‑148 (MQ=255)
                                                                                                                           gtgCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCa                         >  2:483754/1‑148 (MQ=255)
                                                                                                                           gtgCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCa                         <  2:654278/148‑1 (MQ=255)
                                                                                                                             gCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACt                       <  1:2826280/148‑1 (MQ=255)
                                                                                                                             gCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACt                       <  1:2826265/148‑1 (MQ=255)
                                                                                                                                ggTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTcc                    >  2:1501209/1‑148 (MQ=255)
                                                                                                                                 gTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTccc                   >  2:3197290/1‑148 (MQ=255)
                                                                                                                                 gTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTccc                   <  2:76977/148‑1 (MQ=255)
                                                                                                                                  tAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCa                  >  1:2517910/1‑148 (MQ=255)
                                                                                                                                     aTTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATaa               >  2:1390989/1‑148 (MQ=255)
                                                                                                                                     aTTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATaa               <  2:1426265/148‑1 (MQ=255)
                                                                                                                                     aTTATCTGCTTTATAGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATaa               >  2:486625/1‑148 (MQ=255)
                                                                                                                                      ttATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAAc              >  2:1284670/1‑148 (MQ=255)
                                                                                                                                      ttATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAAc              <  1:472567/148‑1 (MQ=255)
                                                                                                                                      ttATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAAc              >  2:847387/1‑148 (MQ=255)
                                                                                                                                      ttATCTGCTTTATAGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAAc              <  2:935825/148‑1 (MQ=255)
                                                                                                                                          cTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGa          >  1:876628/1‑148 (MQ=255)
                                                                                                                                             cTTTATAGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAgtg       >  1:1702673/1‑148 (MQ=255)
                                                                                                                                                 aTCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACt   >  1:1856474/1‑148 (MQ=255)
                                                                                                                                                 aTCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACt   >  1:1046828/1‑148 (MQ=255)
                                                                                                                                                  tCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTg  <  1:390041/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CTGGAGCGGCGAATCGCCCGTCTGGAAACGGGCATGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTG  >  CP009273/2663968‑2664261

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: