Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,601,871 C→A 32.0% P90Q (CCG→CAG)  tam → trans‑aconitate methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,601,8710CA32.0% 375.2 / 160.6 208P90Q (CCG→CAG) tamtrans‑aconitate methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (66/75);  new base A (26/41);  total (92/116)
Fisher's exact test for biased strand distribution p-value = 2.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 210 total aligned reads displayed.

CCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGTACAGATGCCAGATAACTGGCTGGAGCCGACCC  >  CP009273/1601725‑1602017
                                                                                                                                                  |                                                                                                                                                  
cccAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCg                                                                                                                                                   <  1:1438062/148‑1 (MQ=255)
  cAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGt                                                                                                                                                 >  2:2079060/1‑148 (MQ=255)
     gTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACa                                                                                                                                              >  2:2163128/1‑148 (MQ=255)
     gTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACa                                                                                                                                              >  1:2537657/1‑148 (MQ=255)
      tAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTTCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGGACAg                                                                                                                                             >  1:721100/1‑148 (MQ=255)
      tAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAg                                                                                                                                             <  2:2563798/148‑1 (MQ=255)
       aaCAGCACCGCCCTTCTACAACACCGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACAGGTACAgg                                                                                                                                            <  1:1817327/148‑1 (MQ=255)
         cAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCa                                                                                                                                          <  2:3126752/148‑1 (MQ=255)
           gCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACt                                                                                                                                        >  1:1065036/1‑148 (MQ=255)
           gCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACAGGTACAGGCACt                                                                                                                                        <  2:2940533/148‑1 (MQ=255)
           gCACCGCCCTTCTACAACAACGTTGGCCTGCGGACAGGATAACAGGCATAGACTCGTCTCCGGAGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCAGCAACTGGCAACCGGTACAGGCACt                                                                                                                                        >  1:2571641/1‑148 (MQ=255)
            cACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACAGGTACAGGCACTc                                                                                                                                       >  2:775456/1‑148 (MQ=255)
               cGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGAt                                                                                                                                    <  1:2358603/148‑1 (MQ=255)
               cGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACAGGTACAGGCACTCGAt                                                                                                                                    <  1:966900/148‑1 (MQ=255)
                gCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATc                                                                                                                                   <  2:303252/148‑1 (MQ=255)
                gCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATc                                                                                                                                   >  1:2415289/1‑148 (MQ=255)
                gCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATc                                                                                                                                   >  1:1301537/1‑148 (MQ=255)
                  ccTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTg                                                                                                                                 <  1:3276178/148‑1 (MQ=255)
                   cTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACAGGTACAGGCACTCGATCTGa                                                                                                                                >  2:752117/1‑148 (MQ=255)
                     tCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGAtt                                                                                                                              <  2:2635489/148‑1 (MQ=255)
                     tCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGAGCGGAtg                                                                                                                              >  1:2099342/1‑147 (MQ=255)
                      cTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACAGGTACAGGCACTCGATCTGAttt                                                                                                                             <  2:1237407/148‑1 (MQ=255)
                       tACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGAtttt                                                                                                                            >  2:2505257/1‑148 (MQ=255)
                       tACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGAtttt                                                                                                                            <  2:1055182/148‑1 (MQ=255)
                       tACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACAGGTACAGGCACTCGATCTGAtttt                                                                                                                            >  1:2492365/1‑148 (MQ=255)
                        acaacaACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGAttttt                                                                                                                           >  1:2612519/1‑148 (MQ=255)
                         caacaaCGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTg                                                                                                                          >  1:2678351/1‑148 (MQ=255)
                         caacaaCGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTg                                                                                                                          >  1:404612/1‑148 (MQ=255)
                         caacaaCGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACAGGTACAGGCACTCGATCTGATTTTTg                                                                                                                          >  2:1721706/1‑148 (MQ=255)
                          aacaacGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGc                                                                                                                         <  1:1933095/148‑1 (MQ=255)
                          aacaacGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGc                                                                                                                         <  1:3081584/148‑1 (MQ=255)
                          aacaacGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACAGGTACAGGCACTCGATCTGATTTTTGc                                                                                                                         <  1:776957/148‑1 (MQ=255)
                           acaacGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCt                                                                                                                        >  2:2008680/1‑148 (MQ=255)
                            caacGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTa                                                                                                                       <  2:1840572/148‑1 (MQ=255)
                               cGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATg                                                                                                                    <  2:3164782/148‑1 (MQ=255)
                               cGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCAGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTTGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTGTGCTAATg                                                                                                                    <  2:1968860/148‑1 (MQ=255)
                                 ttGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGcc                                                                                                                  <  2:1172466/148‑1 (MQ=255)
                                 ttGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGcc                                                                                                                  <  2:2738623/148‑1 (MQ=255)
                                 ttGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGcc                                                                                                                  >  1:2278437/1‑148 (MQ=255)
                                       tGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTg                                                                                                            <  1:812456/148‑1 (MQ=255)
                                          ggCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCaa                                                                                                         <  1:833251/148‑1 (MQ=255)
                                          ggCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCaa                                                                                                         <  2:2472769/148‑1 (MQ=255)
                                          ggCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACAGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCaa                                                                                                         <  1:2912096/148‑1 (MQ=255)
                                             cAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATgg                                                                                                      <  2:2821030/148‑1 (MQ=255)
                                               ggATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACAGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCt                                                                                                    >  2:2860940/1‑148 (MQ=255)
                                                    aCAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCg                                                                                               >  2:1423827/1‑148 (MQ=255)
                                                    aCAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACAGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCg                                                                                               <  1:775456/148‑1 (MQ=255)
                                                      aGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGAc                                                                                             <  1:561603/148‑1 (MQ=255)
                                                      aGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGAc                                                                                             >  2:3091276/1‑148 (MQ=255)
                                                       ggCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACAGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGAcc                                                                                            >  1:746291/1‑148 (MQ=255)
                                                        gCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCa                                                                                           >  2:2412598/1‑148 (MQ=255)
                                                              aCTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGa                                                                                     <  2:1602046/148‑1 (MQ=255)
                                                               cTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGaa                                                                                    <  1:2969598/148‑1 (MQ=255)
                                                               cTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACAGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGaa                                                                                    <  1:157113/148‑1 (MQ=255)
                                                                tCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAAt                                                                                   <  1:2591207/148‑1 (MQ=255)
                                                                 cGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAAtt                                                                                  <  2:1353610/148‑1 (MQ=255)
                                                                    ctcCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGtt                                                                               >  1:1076093/1‑148 (MQ=255)
                                                                    ctcCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGtt                                                                               >  1:2635788/1‑148 (MQ=255)
                                                                      ccGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTc                                                                             <  1:1749543/148‑1 (MQ=255)
                                                                        ggCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCt                                                                           <  2:1309175/148‑1 (MQ=255)
                                                                         gCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTc                                                                          <  2:919888/148‑1 (MQ=255)
                                                                          cGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCa                                                                         >  1:2674808/1‑148 (MQ=255)
                                                                           gatgatTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACAGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCAt                                                                        <  2:1915296/148‑1 (MQ=255)
                                                                            atgatTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATc                                                                       <  1:1004239/148‑1 (MQ=255)
                                                                              gatTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGAc                                                                                             >  1:467239/1‑124 (MQ=255)
                                                                               atTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTgg                                                                    <  1:745496/148‑1 (MQ=255)
                                                                                ttGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACAGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGt                                                                   <  2:1557038/148‑1 (MQ=255)
                                                                                ttGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACAGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGt                                                                   <  1:1355905/148‑1 (MQ=255)
                                                                                ttGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACAGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATCGTTTCCTCATCTGGt                                                                   <  2:1558244/148‑1 (MQ=255)
                                                                                 tGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACAGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGtt                                                                  <  2:3219126/148‑1 (MQ=255)
                                                                                    tGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACAGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCt                                                               <  2:2950078/148‑1 (MQ=255)
                                                                                      aaGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCttt                                                             <  2:1425162/148‑1 (MQ=255)
                                                                                      aaGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACAGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCttt                                                             >  2:1311617/1‑148 (MQ=255)
                                                                                        gcgcgcAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTAc                                                           >  2:654932/1‑148 (MQ=255)
                                                                                        gcgcgcAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTAc                                                           >  1:1583062/1‑148 (MQ=255)
                                                                                         cgcgcAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACAGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACt                                                          >  1:1429716/1‑148 (MQ=255)
                                                                                           cgcAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACAGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTa                                                        <  2:1181978/148‑1 (MQ=255)
                                                                                             cAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAAt                                                      <  2:3106044/148‑1 (MQ=255)
                                                                                             cAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAAt                                                      <  2:1529062/148‑1 (MQ=255)
                                                                                               gTGCTTTGCCAGACTGCCAGTTCGTGGACGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATcc                                                    >  2:2447718/1‑148 (MQ=255)
                                                                                               gTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACAGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATcc                                                    >  1:760700/1‑148 (MQ=255)
                                                                                                    ttGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAgg                                               >  2:3087577/1‑148 (MQ=255)
                                                                                                    ttGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAgg                                               >  2:1285588/1‑148 (MQ=255)
                                                                                                     tGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACAGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAggg                                              >  1:339635/1‑148 (MQ=255)
                                                                                                     tGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACAGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAggg                                              >  1:339538/1‑148 (MQ=255)
                                                                                                          taCTGCCAGTTGGTGTATGCATAAATCCGAAAGTGGAAACAGTTACAGGAAATCTATGTTTTTTTTGATATTGCCTAACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCTTCTGGTTTATTTATTTAATCCGCTGGGCGTGc                                         <  2:392082/147‑1 (MQ=255)
                                                                                                           aCTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACAGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCt                                        <  1:909948/148‑1 (MQ=255)
                                                                                                              gCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGc                                     >  2:2795045/1‑148 (MQ=255)
                                                                                                                  gTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGTa                                 >  1:1176131/1‑148 (MQ=255)
                                                                                                                  gTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGTa                                 >  1:2147133/1‑148 (MQ=255)
                                                                                                                   ttCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGTAc                                >  1:2085752/1‑148 (MQ=255)
                                                                                                                    tCGTGGAAGCAGACATCCGCAACTGGCAACAGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGTACa                               >  2:154597/1‑148 (MQ=255)
                                                                                                                            gCAGACATCCGCAACTGGCAACAGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGTACAGATGCCAg                       >  2:347785/1‑148 (MQ=255)
                                                                                                                             cAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGTACAGATGCCAGa                      <  1:2435297/148‑1 (MQ=255)
                                                                                                                               gACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGTACAGATGCCAGATa                    >  1:2833352/1‑148 (MQ=255)
                                                                                                                                aCATCCGCAACTGGCAACAGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGTACAGATGCCAGATaa                   <  1:1705893/148‑1 (MQ=255)
                                                                                                                                    ccGCAACTGGCAACAGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGTACAGATGCCAGATAActgg               >  1:16270/1‑148 (MQ=255)
                                                                                                                                         aCTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGTACAGATGCCAGATAACTGGCTGGa          >  2:50152/1‑148 (MQ=255)
                                                                                                                                         aCTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGTACAGATGCCAGATAACTGGCTGGa          >  2:246892/1‑148 (MQ=255)
                                                                                                                                          cTGGCAACAGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGTACAGATGCCAGATAACTGGCTGGAg         >  2:537039/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGTACAGATGCCAGATAACTGGCTGGAGCCGACCC  >  CP009273/1601725‑1602017

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: ‑