Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,303,463 C→A 9.9% T355T (ACG→ACT)  pnp ← polynucleotide phosphorylase/polyadenylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,303,4630CA 9.9% 734.9 / 37.3 212T355T (ACG→ACT) pnppolynucleotide phosphorylase/polyadenylase
Reads supporting (aligned to +/- strand):  ref base C (107/84);  new base A (11/10);  total (118/94)
Fisher's exact test for biased strand distribution p-value = 8.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01

Only 100 of 217 total aligned reads displayed.

ACCGTGACCAATTTCACGACGCTTCGGAGAACCGACCATGCCGGTTTCGCCTACGGAGTACGGAGGGAAGTTGTAGTGGAACAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACATCCAGACCACGGATCATATCTTTTTCACGACCGTCGATACGCGGTTCGCCTGCCAGTACGCGGCTA  >  CP009273/3303316‑3303610
                                                                                                                                                   |                                                                                                                                                   
    tGACCAATTTCACGACGCTTCGGAGAACCGACCATGCCGGTTTCGCCTACGGAGTACGGAGGGAAGTTGTAGTGGAACAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTg                                                                                                                                                 <  2:2538613/148‑1 (MQ=255)
     gACCAATTTCACGACGCTTCGGAGAACCGACCATGCCGGTTTCGCCTACGGAGTACGGAGGGAAGTTGTAGTGGAACAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGc                                                                                                                                                <  1:1092483/148‑1 (MQ=255)
      aCCAATTTCACGACGCTTCGGAGAACCGACCATGCCGGTTTCGCCTACGGAGTACGGAGGGAAGTTGTAGTGGAACAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCg                                                                                                                                               <  2:890708/148‑1 (MQ=255)
       ccAATTTCACGACGCTTCGGAGAACCGACCATGCCGGTTTCGCCTACGGAGTACGGAGGGAAGTTGTAGTGGAACAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCgg                                                                                                                                              <  1:1399885/148‑1 (MQ=255)
       ccAATTTCACGACGCTTCGGAGAACCGACCATGCCGGTTTCGCCTACGGAGTACGGAGGGAAGTTGTAGTGGAACAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCgg                                                                                                                                              <  2:2344643/148‑1 (MQ=255)
          aTTTCACGACGCTTCGGAGAACCGACCATGCCGGTTTCGCCTACGGAGTACGGAGGGAAGTTGTAGTGGAACAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTaa                                                                                                                                           >  2:2536473/1‑148 (MQ=255)
          aTTTCACGACGCTTCGGAGAACCGACCATGCCGGTTTCGCCTACGGAGTACGGAGGGAAGTTGTAGTGGAACAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGAGTTGCGGTaa                                                                                                                                           >  1:386229/1‑148 (MQ=255)
             tCACGACGCTTCGGAGAACCGACCATGCCGGTTTCGCCTACGGAGTACGGAGGGAAGTTGTAGTGGAACAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCa                                                                                                                                        <  2:550322/148‑1 (MQ=255)
              cACGACGCTTCGGAGAACCGACCATGCCGGTTTCGCCTACGGAGTACGGAGGGAAGTTGTAGTGGAACAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAg                                                                                                                                       <  1:3050423/148‑1 (MQ=255)
               acgacgCTTCGGAGAACCGACCATGCCGGTTTCGCCTACGGAGTACGGAGGGAAGTTGTAGTGGAACAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAgc                                                                                                                                      >  1:2686573/1‑148 (MQ=255)
                 gacgCTTCGGAGAACCGACCATGCCGGTTTCGCCTACGGAGTACGGAGGGAAGTTGTAGTGGAACAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAgcgc                                                                                                                                    <  1:2704524/148‑1 (MQ=255)
                 gacgCTTCGGAGAACCGACCATGCCGGTTTCGCCTACGGAGTACGGAGGGAAGTTGTAGTGGAACAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAgcgc                                                                                                                                    <  1:2704315/148‑1 (MQ=255)
                  acgCTTCGGAGAACCGACCATGCCGGTTTCGCCTACGGAGTACGGAGGGAAGTTGTAGTGGAACAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGcc                                                                                                                                   >  2:267402/1‑148 (MQ=255)
                       tCGGAGAACCGACCATGCCGGTTTCGCCTACGGAGTACGGAGGGAAGTTGTAGTGGAACAGGAAGGTCTCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGt                                                                                                                              >  2:1795483/1‑148 (MQ=255)
                        cGGAGAACCGACCATGCCGGTTTCGCCTACGGAGTACGGAGGGAAGTTGTAGTGGAACAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGtt                                                                                                                             <  1:1720355/148‑1 (MQ=255)
                        cGGAGAACCGACCATGCCGGTTTCGCCTACGGAGTACGGAGGGAAGTTGTAGTGGAACAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGtt                                                                                                                             <  2:2413761/148‑1 (MQ=255)
                         ggAGAACCGACCATGCCGGTTTCGCCTACGGAGTACGGAGGGAAGTTGTAGTGGAACAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGttt                                                                                                                            <  1:74019/148‑1 (MQ=255)
                          gagaACCGACCATGCCGGTTTCGCCTACGGAGTACGGAGGGAAGTTGTAGTGGAACAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTc                                                                                                                           >  1:912915/1‑148 (MQ=255)
                                   ccATGCCGGTTTCGCCTACGGAGTACGGAGGGAAGTTGTAGTGGAACAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGt                                                                                                                  <  1:2475326/148‑1 (MQ=255)
                                     aTGCCGGGTTCGCCTACGGAGTACGGAGGGAAGTTGTAGTGGACCAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGAGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGa                                                                                                                <  2:104149/148‑1 (MQ=255)
                                        ccGGTTTCGCCTACGGAGTACGGAGGGAAGTTGTAGTGGAACAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACa                                                                                                             <  1:2081777/148‑1 (MQ=255)
                                        ccGGTTTCGCCTACGGAGTACGGAGGGAAGTTGTAGTGGAACAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACa                                                                                                             <  1:2281837/148‑1 (MQ=255)
                                        ccGGTTTCGCCTACGGAGTACGGAGGGAAGTTGTAGTGGAACAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGAGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAAcc                                                                                                             >  2:2296954/1‑147 (MQ=255)
                                           gTTTCGCCTACGGAGTACGGAGGGAAGTTGTAGTGGAACAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAgcg                                                                                                          <  1:790039/148‑1 (MQ=255)
                                           gTTTCGCCTACGGAGTACGGAGGGAAGTTGTAGTGGAACAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAgcg                                                                                                          >  2:862414/1‑148 (MQ=255)
                                               cGCCTACGGAGTACGGAGGGAAGTTGTAGTGGAACAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGa                                                                                                      <  2:274864/148‑1 (MQ=255)
                                                gCCTACGGAGTACGGAGGGAAGTTGTAGTGGAACAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGaa                                                                                                     <  1:1794134/148‑1 (MQ=255)
                                                 ccTACGGAGTACGGAGGGAAGTTGTAGTGGAACAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAAc                                                                                                    <  1:1188771/148‑1 (MQ=255)
                                                  cTACGGAGTACGGAGGGAAGTTGTAGTGGAACAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAAcc                                                                                                   >  2:357100/1‑148 (MQ=255)
                                                    acggagtacggagGGAAGTTGTAGTGGAACAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGt                                                                                                 >  1:295860/1‑148 (MQ=255)
                                                              gagGGAAGTTGTAGTGGAACAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCgg                                                                                       >  2:3244677/1‑148 (MQ=255)
                                                               agGGAAGTTGTAGTGGAACAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGgc                                                                                      >  1:2430101/1‑148 (MQ=255)
                                                               agGGAAGTTGTAGTGGAACAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGgc                                                                                      >  2:365736/1‑148 (MQ=255)
                                                                 ggAAGTTGTAGTGGAACAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGgcag                                                                                    >  1:1038677/1‑148 (MQ=255)
                                                                          aGTGGAACAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGAGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGt                                                                           <  2:1614445/148‑1 (MQ=255)
                                                                            tGGAACAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTAc                                                                         >  1:515223/1‑148 (MQ=255)
                                                                            tGGAACAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTAc                                                                         >  2:1451217/1‑148 (MQ=255)
                                                                             ggAACAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACg                                                                        >  1:40923/1‑148 (MQ=255)
                                                                              gAACAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGc                                                                       >  2:2493251/1‑148 (MQ=255)
                                                                              gAACAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGc                                                                       >  1:3018867/1‑148 (MQ=255)
                                                                              gAACAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGAGAGTACGCGGCAGCACGCCAGTACGc                                                                       >  2:2826951/1‑148 (MQ=255)
                                                                                aCAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGcac                                                                     <  1:1869661/148‑1 (MQ=255)
                                                                                aCAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGcac                                                                     <  2:1418193/148‑1 (MQ=255)
                                                                                 cAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGcaca                                                                    >  1:1134804/1‑148 (MQ=255)
                                                                                 cAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGAGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGcaca                                                                    <  2:812022/148‑1 (MQ=255)
                                                                                  aGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACAt                                                                   <  1:2849019/148‑1 (MQ=255)
                                                                                   ggAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACATc                                                                  >  2:2687589/1‑148 (MQ=255)
                                                                                    gAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACATcc                                                                 >  2:2605832/1‑148 (MQ=255)
                                                                                      aGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACATCCAg                                                               >  1:994709/1‑148 (MQ=255)
                                                                                      aGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACATCCAg                                                               <  1:1921460/148‑1 (MQ=255)
                                                                                       ggTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACATCCAGa                                                              <  2:1548776/148‑1 (MQ=255)
                                                                                        gTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACATCCAGAc                                                             >  1:3229384/1‑148 (MQ=255)
                                                                                        gTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACATCCAGAc                                                             <  1:2360040/148‑1 (MQ=255)
                                                                                         gaTCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACATCCAGAcc                                                            <  2:1851659/147‑1 (MQ=255)
                                                                                            cGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACATCCAGACCACg                                                         >  2:2589358/1‑148 (MQ=255)
                                                                                            cGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACATCCAGACCACg                                                         >  2:2466965/1‑148 (MQ=255)
                                                                                               tACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACATCCAGACCACGGAt                                                      >  2:361490/1‑148 (MQ=255)
                                                                                               tACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACATCCAGACCACGGAt                                                      >  2:1770946/1‑148 (MQ=255)
                                                                                               gaCGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGGACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACATCCAGACCACGGAt                                                      <  2:619443/147‑1 (MQ=255)
                                                                                                aCGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACATCCAGACCACGGATc                                                     >  2:2147055/1‑148 (MQ=255)
                                                                                                 cGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACATCCAGACCACGGATCa                                                    >  1:2694227/1‑148 (MQ=255)
                                                                                                  gTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGAGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACATCCAGACCACGGATCat                                                   <  1:1415599/148‑1 (MQ=255)
                                                                                                   ttCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACATCCAGACCACGGATCata                                                  <  2:2908712/148‑1 (MQ=255)
                                                                                                    tCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACATCCAGACCACGGATCatat                                                 <  1:63075/148‑1 (MQ=255)
                                                                                                     cGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACATCCAGACCACGGATCATATc                                                <  1:584278/148‑1 (MQ=255)
                                                                                                       cccATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACATCCAGACCACGGATCATATCtt                                              >  2:2913592/1‑148 (MQ=255)
                                                                                                          aTCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACATCCAGACCACGGATCATATCttttt                                           >  2:3004552/1‑148 (MQ=255)
                                                                                                          aTCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACATCCAGACCACGGATCATATCttttt                                           >  2:1358222/1‑148 (MQ=255)
                                                                                                           tCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACATCCAGACCACGGATCATATCTTTTTc                                          >  2:3276008/1‑148 (MQ=255)
                                                                                                           tCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGAGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGACCAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACATCCAGACCACGGATCATATCTTTTTc                                          >  2:666858/1‑148 (MQ=255)
                                                                                                                 catcaAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACATCCAGACCACGGATCATATCTTTTTCACGAcc                                    >  2:1155101/1‑148 (MQ=255)
                                                                                                                  atcaAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACATCCAGACCACGGATCATATCTTTTTCACGACCg                                   >  2:3151455/1‑148 (MQ=255)
                                                                                                                  atcaAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACATCCAGACCACGGATCATATCTTTTTCACGACCg                                   <  1:865667/148‑1 (MQ=255)
                                                                                                                   tcaAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACATCCAGACCACGGATCATATCTTTTTCACGACCGt                                  >  2:3046473/1‑148 (MQ=255)
                                                                                                                   tcaAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACATCCAGACCACGGATCATATCTTTTTCACGACCGt                                  >  2:3292675/1‑148 (MQ=255)
                                                                                                                     aaGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACATCCAGACCACGGATCATATCTTTTTCACGACCGTCg                                <  1:168165/148‑1 (MQ=255)
                                                                                                                     aaGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACATCCAGACCACGGATCATATCTTTTTCACGACCGTCg                                <  2:1001439/148‑1 (MQ=255)
                                                                                                                     aaGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACATCCAGACCACGGATCATATCTTTTTCACGACCGTCg                                >  1:3193325/1‑148 (MQ=255)
                                                                                                                      aGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACATCCAGACCACGGATCATATCTTTTTCACGACCGTCGa                               <  1:1451217/148‑1 (MQ=255)
                                                                                                                       gAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACATCCAGACCACGGATCATATCTTTTTCACGACCGTCGAt                              <  1:347726/148‑1 (MQ=255)
                                                                                                                        aaCCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACATCCAGACCACGGATCATATCTTTTTCACGACCGTCGATa                             >  1:2409098/1‑148 (MQ=255)
                                                                                                                         aCCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACATCCAGACCACGGATCATATCTTTTTCACGACCGTCGATAc                            >  2:129325/1‑148 (MQ=255)
                                                                                                                         aCCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACATCCAGACCACGGATCATATCTTTTTCACGACCGTCGATAc                            >  2:430806/1‑148 (MQ=255)
                                                                                                                          ccTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACATCCAGACCACGGATCATATCTTTTTCACGACCGTCGATAcg                           >  2:2091694/1‑148 (MQ=255)
                                                                                                                           cTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACATCCAGACCACGGATCATATCTTTTTCACGACCGTCGATAcgc                          >  1:1593987/1‑148 (MQ=255)
                                                                                                                                cgTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACATCCAGACCACGGATCATATCTTTTTCACGACCGTCGATACGCGGTTc                     >  2:359612/1‑148 (MQ=255)
                                                                                                                                   cACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACATCCAGACCACGGATCATATCTTTTTCACGACCGTCGATACGCGGTTCGcc                  >  1:1228326/1‑148 (MQ=255)
                                                                                                                                   cACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACATCCAGACCACGGATCATATCTTTTTCACGACCGTCGATACGCGGTTCGcc                  >  1:1772879/1‑148 (MQ=255)
                                                                                                                                    aCGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACATCCAGACCACGGATCATATCTTTTTCACGACCGTCGATACGCGGTTCGCCt                 >  2:146226/1‑148 (MQ=255)
                                                                                                                                       aGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACATCCAGACCACGGATCATATCTTTTTCACGACCGTCGATACGCGGTTCGCCTGcc              >  1:1673185/1‑148 (MQ=255)
                                                                                                                                        gCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACATCCAGACCACGGATCATATCTTTTTCACGACCGTCGATACGCGGTTCGCCTGCCa             >  2:452629/1‑148 (MQ=255)
                                                                                                                                            tACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACATCCAGACCACGGATCATATCTTTTTCACGACCGTCGATACGCGGTTCGCCTGCCAGTAc         >  1:951498/1‑148 (MQ=255)
                                                                                                                                            tACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACATCCAGACCACGGATCATATCTTTTTCACGACCGTCGATACGCGGTTCGCCTGCCAGTAc         >  1:2604947/1‑148 (MQ=255)
                                                                                                                                                cAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACATCCAGACCACGGATCATATCTTTTTCACGACCGTCGATACGCGGTTCGCCTGCCAGTACGCgg     >  2:3192026/1‑148 (MQ=255)
                                                                                                                                                cAGAGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACATCCAGACCACGGATCATATCTTTTTCACGACCGTCGATACGCGGTTCGCCTGCCAGTACGCgg     <  1:1823073/148‑1 (MQ=255)
                                                                                                                                                 aGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACATCCAGACCACGGATCATATCTTTTTCACGACCGTCGATACGCGGTTCGCCTGCCAGTACGCGGc    >  1:1236282/1‑148 (MQ=255)
                                                                                                                                                  gCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACATCCAGACCACGGATCATATCTTTTTCACGACCGTCGATACGCGGTTCGCCTGCCAGTACGCGGCt   >  2:2329616/1‑148 (MQ=255)
                                                                                                                                                  gCGTTGCGGTAACAAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACATCCAGACCACGGATCATATCTGTTTCACGACCGTCGATACGCTGTTCGCCTGCCCGTACGCGGCt   >  2:1721458/1‑148 (MQ=255)
                                                                                                                                                   cGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACATCCAGACCACGGATCATATCTTTTTCACGACCGTCGATACGCGGTTCGCCTGCCAGTACGCGGCTa  >  1:596879/1‑148 (MQ=255)
                                                                                                                                                   cGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACATCCAGACCACGGATCATATCTTTTTCACGACCGTCGATACGCGGTTCGCCTGCCAGTACGCGGCTa  >  2:1939436/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ACCGTGACCAATTTCACGACGCTTCGGAGAACCGACCATGCCGGTTTCGCCTACGGAGTACGGAGGGAAGTTGTAGTGGAACAGGAAGGTATCGGTACGTTCGCCCATCAGTTCATCAAGAACCTGCGCGTCACGAGCAGTACCCAGCGTTGCGGTAACCAGCGCCTGCGTTTCACCACGGGTGAACAGCGCAGAACCGTGAGTACGCGGCAGCACGCCAGTACGCACATCCAGACCACGGATCATATCTTTTTCACGACCGTCGATACGCGGTTCGCCTGCCAGTACGCGGCTA  >  CP009273/3303316‑3303610

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