Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,601,871 C→A 21.8% 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,8710CA21.8% 538.4 / 100.4 210P90Q (CCG→CAG) tamtrans‑aconitate methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (76/88);  new base A (16/30);  total (92/118)
Fisher's exact test for biased strand distribution p-value = 1.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 217 total aligned reads displayed.

GCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGTACAGATGCCAGATAACTGGCTGGAGCCGACCCA  >  CP009273/1601724‑1602018
                                                                                                                                                   |                                                                                                                                                   
gCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGACGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAAcc                                                                                                                                                     >  2:3243716/1‑148 (MQ=255)
   cAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGt                                                                                                                                                  <  2:536372/148‑1 (MQ=255)
   cAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACAGGt                                                                                                                                                  <  2:2606329/148‑1 (MQ=255)
     ggTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTAc                                                                                                                                                >  2:472586/1‑148 (MQ=255)
     ggTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACAGGTAc                                                                                                                                                <  1:2276695/148‑1 (MQ=255)
      gTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACa                                                                                                                                               <  1:2270396/148‑1 (MQ=255)
        aaCAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAgg                                                                                                                                             >  2:3577513/1‑148 (MQ=255)
          cAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCa                                                                                                                                           <  2:3746086/148‑1 (MQ=255)
          cAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCa                                                                                                                                           <  1:2710372/148‑1 (MQ=255)
            gCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACt                                                                                                                                         >  1:1623614/1‑148 (MQ=255)
             cACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTc                                                                                                                                        >  2:665770/1‑148 (MQ=255)
               ccGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGa                                                                                                                                      <  2:722222/148‑1 (MQ=255)
               ccGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACAGGTACAGGCACTCGa                                                                                                                                      <  2:1486503/148‑1 (MQ=255)
                cGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGAt                                                                                                                                     >  2:2962978/1‑148 (MQ=255)
                cGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGAt                                                                                                                                     <  1:371495/148‑1 (MQ=255)
                cGCCCTTCTACAAAAACGTTGTCCTGCGGCAATGATAACAGTCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGGAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGAt                                                                                                                                     <  1:2569588/148‑1 (MQ=255)
                 gCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATc                                                                                                                                    <  1:3212662/148‑1 (MQ=255)
                  cccTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCt                                                                                                                                   >  1:1357214/1‑148 (MQ=255)
                    cTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGa                                                                                                                                 <  1:2623595/148‑1 (MQ=255)
                    cTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACAGGTACAGGCACTCGATCTGa                                                                                                                                 <  1:401169/148‑1 (MQ=255)
                        tACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGAtttt                                                                                                                             <  1:717655/148‑1 (MQ=255)
                        tACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGAtttt                                                                                                                             <  1:718137/148‑1 (MQ=255)
                        tACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGAtttt                                                                                                                             <  1:1611848/148‑1 (MQ=255)
                           aacaacGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGc                                                                                                                          <  1:1149772/148‑1 (MQ=255)
                             caacGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTa                                                                                                                        >  2:2723899/1‑148 (MQ=255)
                              aacGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACAGGTACAGGCACTCGATCTGATTTTTGCTaa                                                                                                                       <  2:3103347/148‑1 (MQ=255)
                                  ttGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACAGGTACAGGCACTCGATCTGATTTTTGCTAATGcc                                                                                                                   >  1:989810/1‑148 (MQ=255)
                                     gCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCa                                                                                                                <  2:1115545/148‑1 (MQ=255)
                                     gCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCa                                                                                                                <  2:1418255/148‑1 (MQ=255)
                                        tGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACAGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTg                                                                                                             <  1:1413792/148‑1 (MQ=255)
                                         gCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGc                                                                                                            >  2:1583299/1‑148 (MQ=255)
                                              ctggATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGTAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATgg                                                                                                       <  2:1355470/146‑1 (MQ=255)
                                              cAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATgg                                                                                                       >  2:3066739/1‑148 (MQ=255)
                                              cAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATgg                                                                                                       <  1:2755777/148‑1 (MQ=255)
                                               aGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGc                                                                                                      <  1:726862/148‑1 (MQ=255)
                                                 gATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTg                                                                                                    >  2:2303105/1‑148 (MQ=255)
                                                 gATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTg                                                                                                    >  1:2209525/1‑148 (MQ=255)
                                                      cAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGa                                                                                               >  1:2672358/1‑148 (MQ=255)
                                                           aTCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACAGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACt                                                                                          <  1:3763593/148‑1 (MQ=255)
                                                               aCTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGa                                                                                      <  1:2519960/148‑1 (MQ=255)
                                                                cTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACAGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGaa                                                                                     <  2:2977510/148‑1 (MQ=255)
                                                                cTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGAAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGaa                                                                                     <  1:269435/148‑1 (MQ=255)
                                                                 tCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACAGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAAt                                                                                    >  1:1272461/1‑148 (MQ=255)
                                                                   gTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGAGATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTg                                                                                  <  2:1848595/148‑1 (MQ=255)
                                                                     ctcCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACAGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGtt                                                                                >  1:3358589/1‑148 (MQ=255)
                                                                        cGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTcc                                                                             <  2:2221687/148‑1 (MQ=255)
                                                                        cGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTcc                                                                             <  2:1520870/148‑1 (MQ=255)
                                                                            gatgatTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCAt                                                                         <  2:1782925/148‑1 (MQ=255)
                                                                            gatgatTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACAGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCAt                                                                         <  2:2406751/148‑1 (MQ=255)
                                                                              tgatTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCt                                                                       >  2:496605/1‑148 (MQ=255)
                                                                              tgatTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCt                                                                       >  2:2885961/1‑148 (MQ=255)
                                                                              tgatTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCt                                                                       <  1:3765427/148‑1 (MQ=255)
                                                                                atTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTgg                                                                     >  1:3402887/1‑148 (MQ=255)
                                                                                atTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTgg                                                                     <  2:917759/148‑1 (MQ=255)
                                                                                atTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACAGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTgg                                                                     <  1:1889725/148‑1 (MQ=255)
                                                                                   gCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGttt                                                                  <  2:2509464/148‑1 (MQ=255)
                                                                                   gCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGttt                                                                  >  1:355095/1‑148 (MQ=255)
                                                                                     tGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGTCCACTACGAATGTTTTCCTCATCTGTTTTCt                                                                <  1:381579/148‑1 (MQ=255)
                                                                                        aGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTa                                                             <  1:1991830/148‑1 (MQ=255)
                                                                                        aGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTa                                                             >  1:1789918/1‑148 (MQ=255)
                                                                                            cgcAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCTCTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTa                                                         <  1:676344/148‑1 (MQ=255)
                                                                                            cgcAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTa                                                         <  2:1019337/148‑1 (MQ=255)
                                                                                            cgcAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTa                                                         <  1:2592973/148‑1 (MQ=255)
                                                                                             gcAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTaa                                                        <  1:2962136/148‑1 (MQ=255)
                                                                                                  gCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGc                                                   >  1:2974906/1‑148 (MQ=255)
                                                                                                   cTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACAGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCa                                                  <  1:1638748/148‑1 (MQ=255)
                                                                                                    tttGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACAGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAg                                                 >  1:591258/1‑148 (MQ=255)
                                                                                                          aGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTg                                           >  1:124360/1‑148 (MQ=255)
                                                                                                          aGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTg                                           <  2:1357214/148‑1 (MQ=255)
                                                                                                              tGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTgg                                       <  2:1668883/148‑1 (MQ=255)
                                                                                                              tGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTgg                                       <  2:3562094/148‑1 (MQ=255)
                                                                                                               gCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGc                                      >  1:3312731/1‑148 (MQ=255)
                                                                                                                ccAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCa                                     >  2:1161730/1‑148 (MQ=255)
                                                                                                                ccAGTTCGTGGAAGCAGACATCCGCAACTGGCAACAGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCa                                     >  2:1889022/1‑148 (MQ=255)
                                                                                                                  aGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGt                                   <  1:1146777/148‑1 (MQ=255)
                                                                                                                  aGTTCGTGGAAGCAGACATCCGCAACTGGCAACAGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGt                                   >  2:2781030/1‑148 (MQ=255)
                                                                                                                   gTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGTa                                  <  1:2303105/148‑1 (MQ=255)
                                                                                                                   gTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGTa                                  <  2:3046591/148‑1 (MQ=255)
                                                                                                                       gTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGTACAGa                              <  1:226658/148‑1 (MQ=255)
                                                                                                                       gTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGTACAGa                              <  1:496605/148‑1 (MQ=255)
                                                                                                                         ggAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGTACAGATg                            <  2:2736613/148‑1 (MQ=255)
                                                                                                                            aGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGTACAGATGCCa                         >  2:3033053/1‑148 (MQ=255)
                                                                                                                            aGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGTACAGATGCCa                         >  2:3687771/1‑148 (MQ=255)
                                                                                                                             gCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGTACAGATGCCAg                        >  2:2608215/1‑148 (MQ=255)
                                                                                                                             gCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGTACAGATGCCAg                        >  2:1398315/1‑148 (MQ=255)
                                                                                                                              cAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGTACAGATGCCAGa                       >  1:2122829/1‑148 (MQ=255)
                                                                                                                                gACATCCGCAACTGGCAACAGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGTACAGATGCCAGATa                     >  1:2507575/1‑148 (MQ=255)
                                                                                                                                 aCATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGTACAGATGCCAGATaa                    >  1:658155/1‑148 (MQ=255)
                                                                                                                                 aCATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGTACAGATGCCAGATaa                    >  2:3194513/1‑148 (MQ=255)
                                                                                                                                   aTCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGTACAGATGCCAGATAAct                  <  1:1671952/148‑1 (MQ=255)
                                                                                                                                            tGGCAACAGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGTACAGATGCCAGATAACTGGCTGGAGc         <  2:2073110/148‑1 (MQ=255)
                                                                                                                                             ggCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGTACAGATGCCAGATAACTGGCTGGAGcc        >  2:3196319/1‑148 (MQ=255)
                                                                                                                                             ggCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGTACAGATGCCAGATAACTGGCTGGAGcc        >  2:3197021/1‑148 (MQ=255)
                                                                                                                                              gCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGTACAGATGCCAGATAACTGGCTGGAGCCg       >  1:904983/1‑148 (MQ=255)
                                                                                                                                               cAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGTACAGATGCCAGATAACTGGCTGGAGCCGa      >  2:950610/1‑148 (MQ=255)
                                                                                                                                               cAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGTACAGATGCCAGATAACTGGCTGGAGCCGa      >  2:1435706/1‑148 (MQ=255)
                                                                                                                                                 aCCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGTACAGATGCCAGATAACTGGCTGGAGCCGAcc    >  2:1255602/1‑148 (MQ=255)
                                                                                                                                                 aCCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGTACAGATGCCAGATAACTGGCTGGAGCCGAcc    <  2:2056242/148‑1 (MQ=255)
                                                                                                                                                  ccGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGTACAGATGCCAGATAACTGGCTGGAGCCGAccc   >  1:3449747/1‑148 (MQ=255)
                                                                                                                                                   cGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGTACAGATGCCAGATAACTGGCTGGAGCCGACCCa  <  2:3214298/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGTACAGATGCCAGATAACTGGCTGGAGCCGACCCA  >  CP009273/1601724‑1602018

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: