Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,093,612 C→A 93.5% E380* (GAA→TAA)  gnd ← 6‑phosphogluconate dehydrogenase, decarboxylating

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,093,6120CA93.5% 718.1 / 35.2 247E380* (GAA→TAA) gnd6‑phosphogluconate dehydrogenase, decarboxylating
Reads supporting (aligned to +/- strand):  ref base C (6/10);  new base A (128/103);  total (134/113)
Fisher's exact test for biased strand distribution p-value = 1.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 247 total aligned reads displayed.

CGTAATAGGCAACCGCTGCGGAGAAGGTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGAA  >  CP009273/2093465‑2093757
                                                                                                                                                   |                                                                                                                                                 
cGTAATAGGCAACCGCTGCGGAGAAGGTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAgg                                                                                                                                                 <  2:2027875/150‑1 (MQ=255)
cGTAATAGGCAACCGCTGCGGAGAAGGTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAgg                                                                                                                                                 <  2:1471023/150‑1 (MQ=255)
 gTAATAGGCAACCGCTGCGGAGAAGGTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGc                                                                                                                                                >  2:1406267/1‑150 (MQ=255)
  tAATAGGCAACCGCTGCGGAGAAGGTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCa                                                                                                                                               <  2:3222311/150‑1 (MQ=255)
  tAATAGGCAACCGCTGCGGAGAAGGTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCa                                                                                                                                               <  2:5144469/150‑1 (MQ=255)
     tAGGCAACCGCTGCGGAGAAGGTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATaa                                                                                                                                            >  2:4348816/1‑150 (MQ=255)
     tAGGCAACCGCTGCGGAGAAGGTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGGGGATTTTAGGCATaa                                                                                                                                            >  2:1579238/1‑150 (MQ=255)
      aGGCAACCGCTGCGGAGAAGGTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAg                                                                                                                                           <  1:4308421/150‑1 (MQ=255)
        gCAACCGCTGCGGAGAAGGTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCa                                                                                                                                         >  2:4845667/1‑150 (MQ=255)
        gCAACCGCTGCGGAGAAGGTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCa                                                                                                                                         <  1:2037112/150‑1 (MQ=255)
          aaCCGCTGCGGAGAAGGTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATc                                                                                                                                       >  2:1264945/1‑150 (MQ=255)
             cGCTGCGGAGAAGGTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGt                                                                                                                                    >  2:901167/1‑150 (MQ=255)
             cGCTGCGGAGAAGGTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGt                                                                                                                                    <  2:3639087/150‑1 (MQ=255)
                tGCGGAGAAGGTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCACCAGGTTAGCGATCGGTGGATGTTCGGCCTAAGCAACGGAGAt                                                                                                                                 >  2:5078461/1‑150 (MQ=255)
                   ggAGAAGGTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGAtttt                                                                                                                              >  2:5098112/1‑150 (MQ=255)
                         ggTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAg                                                                                                                        >  1:4513087/1‑150 (MQ=255)
                              gAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTATCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACt                                                                                                                   <  1:3975118/150‑1 (MQ=255)
                               aaCCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTg                                                                                                                  <  2:2818848/150‑1 (MQ=255)
                                aCCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGACAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTgc                                                                                                                 <  1:4649235/150‑1 (MQ=255)
                                   ggAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCTGGAACTGCGCa                                                                                                              >  1:103565/1‑150 (MQ=255)
                                   ggAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCa                                                                                                              >  1:2108979/1‑150 (MQ=255)
                                    gAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCAc                                                                                                             >  1:2691043/1‑150 (MQ=255)
                                    gAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCAc                                                                                                             >  2:2477084/1‑150 (MQ=255)
                                     aaTACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACg                                                                                                            >  2:930459/1‑150 (MQ=255)
                                     aaTACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACg                                                                                                            >  1:557621/1‑150 (MQ=255)
                                      aTACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACgg                                                                                                           >  1:3850590/1‑150 (MQ=255)
                                      aTACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACgg                                                                                                           <  2:4045272/150‑1 (MQ=255)
                                        aCCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGgat                                                                                                         <  2:500484/150‑1 (MQ=255)
                                          cGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGgatga                                                                                                       >  2:2779665/1‑150 (MQ=255)
                                           gTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGgatgat                                                                                                      >  1:751874/1‑150 (MQ=255)
                                             tCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGc                                                                                                    >  2:3002757/1‑150 (MQ=255)
                                              cTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGca                                                                                                   <  1:3619124/150‑1 (MQ=255)
                                              cTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGca                                                                                                   <  1:1172226/150‑1 (MQ=255)
                                                gTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGcagc                                                                                                 <  2:2684322/150‑1 (MQ=255)
                                                 tACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGcagcc                                                                                                <  2:3819306/150‑1 (MQ=255)
                                                  aCTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGcagcca                                                                                               <  2:3687919/150‑1 (MQ=255)
                                                    tGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGcagccagc                                                                                             <  2:2882964/150‑1 (MQ=255)
                                                    tGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGcagccagc                                                                                             <  2:1696022/150‑1 (MQ=255)
                                                    tGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGcagccagc                                                                                             >  1:3248771/1‑150 (MQ=255)
                                                     gCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCa                                                                                            <  1:670553/150‑1 (MQ=255)
                                                      cATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCAc                                                                                           >  1:469719/1‑150 (MQ=255)
                                                      cATAAGCAACGACATCACGCAGCGCCTGCGGGTAGTCAGCGGCAATTTGCTTGAAGAACGGAGCCAGCAACAGGTTAGCGATCGGTGGATTTTAGGCATAAGCATCGGTGAATTTCTCCCGGACCTGCGCACGGCTGATGCAGCCAGCAc                                                                                           >  1:1041373/1‑150 (MQ=255)
                                                       aTAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACg                                                                                          >  1:1786859/1‑150 (MQ=255)
                                                         aaGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGa                                                                                        <  2:2570651/150‑1 (MQ=255)
                                                             aaCGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAAt                                                                                    >  2:4692714/1‑150 (MQ=255)
                                                              aCGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATc                                                                                   >  1:3303467/1‑150 (MQ=255)
                                                               cGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCt                                                                                  <  2:446273/150‑1 (MQ=255)
                                                                  cATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTcg                                                                               >  1:447086/1‑150 (MQ=255)
                                                                  cATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTcg                                                                               >  1:1407624/1‑150 (MQ=255)
                                                                     cACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGa                                                                            >  1:912031/1‑150 (MQ=255)
                                                                      aCGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGAt                                                                           >  1:3896799/1‑150 (MQ=255)
                                                                       cGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGAtt                                                                          <  1:1812403/150‑1 (MQ=255)
                                                                        gCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGAttt                                                                         <  2:4700692/150‑1 (MQ=255)
                                                                        gCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGAttt                                                                         >  1:2946880/1‑150 (MQ=255)
                                                                         cAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTc                                                                        <  1:580743/150‑1 (MQ=255)
                                                                         cAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTc                                                                        >  2:2724788/1‑150 (MQ=255)
                                                                           gcgcCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGc                                                                      >  1:2447581/1‑150 (MQ=255)
                                                                            cgcCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGcc                                                                     <  1:2237980/150‑1 (MQ=255)
                                                                                   tgGTAGTCATCGGCAATGTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTc                                                              <  2:3077302/150‑1 (MQ=255)
                                                                                    ggTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCa                                                             >  2:2154977/1‑150 (MQ=255)
                                                                                    ggTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCa                                                             <  1:1292907/150‑1 (MQ=255)
                                                                                     gTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAg                                                            <  2:3979866/150‑1 (MQ=255)
                                                                                          cATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATccc                                                       >  1:1660876/1‑150 (MQ=255)
                                                                                           aTCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCa                                                      <  2:1010940/150‑1 (MQ=255)
                                                                                               gCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTg                                                  <  2:619560/150‑1 (MQ=255)
                                                                                                      gCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCtt                                           <  1:3824712/150‑1 (MQ=255)
                                                                                                       cTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTc                                          <  2:3409956/150‑1 (MQ=255)
                                                                                                         tGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGACAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAg                                        >  2:397530/1‑150 (MQ=255)
                                                                                                             gTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGacgc                                    <  1:1374446/150‑1 (MQ=255)
                                                                                                              tACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGacgca                                   <  2:1316570/150‑1 (MQ=255)
                                                                                                              tACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTCGTGCAGATCCCAGTTGTACTCTTCAGacgca                                   >  2:3953403/1‑150 (MQ=255)
                                                                                                               aCGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACAGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCCTAGTTCAGATCCCAGTTGTACTCTTCAGacgcag                                  >  1:4153416/1‑150 (MQ=255)
                                                                                                                 ggAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGacgcagca                                >  2:2968631/1‑150 (MQ=255)
                                                                                                                     ccAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGacgcagcacgca                            >  1:1993246/1‑150 (MQ=255)
                                                                                                                     ccAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGacgcagcacgca                            <  1:1271345/150‑1 (MQ=255)
                                                                                                                     ccAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGACCCCAGTTGTACTCTTCAGacgcagcacgca                            <  1:3000504/150‑1 (MQ=255)
                                                                                                                      cAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGacgcagcacgcag                           <  2:2358208/150‑1 (MQ=255)
                                                                                                                      cAGCAACAGGTTAGCGATCTATGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGacgcagcacgcag                           >  1:1960423/1‑150 (MQ=255)
                                                                                                                        gCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCt                         <  1:175503/150‑1 (MQ=255)
                                                                                                                          aaCAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTga                       >  1:727508/1‑150 (MQ=255)
                                                                                                                             aGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCATGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATATTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTgagag                    <  2:2306758/150‑1 (MQ=255)
                                                                                                                             aGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTgagag                    <  1:2470125/150‑1 (MQ=255)
                                                                                                                               gTTAGCGATCTGTGGATTTGAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGaa                  <  2:1756501/150‑1 (MQ=255)
                                                                                                                               gTTAGCGATCTGGGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGaa                  <  1:4928256/150‑1 (MQ=255)
                                                                                                                                ttAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAg                 <  2:3875655/150‑1 (MQ=255)
                                                                                                                                  aGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGcc               <  1:2358829/150‑1 (MQ=255)
                                                                                                                                  aGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGcc               >  2:4195851/1‑150 (MQ=255)
                                                                                                                                    cGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCt             >  1:3049179/1‑150 (MQ=255)
                                                                                                                                      aTCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTgg           <  1:2135007/150‑1 (MQ=255)
                                                                                                                                        cTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGc         >  1:3463749/1‑150 (MQ=255)
                                                                                                                                        cTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGc         >  1:2163082/1‑150 (MQ=255)
                                                                                                                                        cTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGc         <  2:2908745/150‑1 (MQ=255)
                                                                                                                                         tgtgGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCg        >  1:3552968/1‑150 (MQ=255)
                                                                                                                                         tgtgGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCg        <  1:3702973/150‑1 (MQ=255)
                                                                                                                                         tgtgGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCg        <  1:1579305/150‑1 (MQ=255)
                                                                                                                                          gtgGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGt       <  2:1658862/150‑1 (MQ=255)
                                                                                                                                          gtgGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGt       <  1:3505662/150‑1 (MQ=255)
                                                                                                                                          gtgGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGt       >  1:4289605/1‑150 (MQ=255)
                                                                                                                                           tgGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTa      <  1:4845667/150‑1 (MQ=255)
                                                                                                                                             gATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAg    >  1:4223801/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CGTAATAGGCAACCGCTGCGGAGAAGGTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGAA  >  CP009273/2093465‑2093757

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: