Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,093,612 C→A 68.9% 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,6120CA68.9% 181.6 / 202.3 193E380* (GAA→TAA) gnd6‑phosphogluconate dehydrogenase, decarboxylating
Reads supporting (aligned to +/- strand):  ref base C (32/28);  new base A (58/75);  total (90/103)
Fisher's exact test for biased strand distribution p-value = 2.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

Only 100 of 194 total aligned reads displayed.

CGTAATAGGCAACCGCTGCGGAGAAGGTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGA  >  CP009273/2093465‑2093756
                                                                                                                                                   |                                                                                                                                                
cGTAATAGGCAACCGCTGCGGAGAAGGTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTc                                                                                                                                                  <  2:419479/148‑1 (MQ=255)
  tAATAGGCAACCGCTGCGGAGAAGGTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCgg                                                                                                                                                >  2:2676577/1‑148 (MQ=255)
   aaTAGGCAACCGCTGCGGAGAAGGTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGc                                                                                                                                               >  1:1274640/1‑148 (MQ=255)
   aaTAGGCAACCGCTGCGGAGAAGGTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGc                                                                                                                                               <  1:130568/148‑1 (MQ=255)
    aTAGGCAACCGCTGCGGAGAAGGTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCa                                                                                                                                              >  1:1580376/1‑148 (MQ=255)
    aTAGGCAACCGCTGCGGAGAAGGTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCa                                                                                                                                              <  1:2527468/148‑1 (MQ=255)
      aGGCAACCGCTGCGGAGAAGGTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATa                                                                                                                                            >  1:282662/1‑148 (MQ=255)
      aGGCAACCGCTGCGGAGAAGGTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATa                                                                                                                                            >  1:2085409/1‑148 (MQ=255)
        gCAACCGCTGCGGAGAAGGTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAg                                                                                                                                          <  2:2615713/148‑1 (MQ=255)
            ccGCTGCGGAGAAGGTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATc                                                                                                                                      <  2:954432/148‑1 (MQ=255)
            ccGCTGCGGAGAAGGTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATc                                                                                                                                      >  1:674099/1‑148 (MQ=255)
                tGCGGAGAAGGTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTg                                                                                                                                  <  2:1688993/148‑1 (MQ=255)
                 gCGGAGAAGGTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGa                                                                                                                                 <  1:123949/148‑1 (MQ=255)
                      gaAGGTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTc                                                                                                                            >  1:3457172/1‑148 (MQ=255)
                       aaGGTCGGAACCGGAATACCGGTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCt                                                                                                                           >  1:266553/1‑148 (MQ=255)
                         ggTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGc                                                                                                                         >  2:3078822/1‑148 (MQ=255)
                         ggTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGc                                                                                                                         <  2:323037/148‑1 (MQ=255)
                         ggTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGc                                                                                                                         <  2:1650700/148‑1 (MQ=255)
                         ggTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGc                                                                                                                         <  2:1650660/148‑1 (MQ=255)
                           tCGGAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAg                                                                                                                       <  1:3040930/148‑1 (MQ=255)
                             ggAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGa                                                                                                                     >  1:311470/1‑148 (MQ=255)
                             ggAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGa                                                                                                                     >  2:752979/1‑148 (MQ=255)
                             ggAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGa                                                                                                                     <  1:752979/148‑1 (MQ=255)
                              gAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGaa                                                                                                                    >  2:97610/1‑148 (MQ=255)
                               aaCCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAAc                                                                                                                   <  1:880020/148‑1 (MQ=255)
                                 ccGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTg                                                                                                                 <  2:1957326/148‑1 (MQ=255)
                                  cGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTgc                                                                                                                <  1:1140509/148‑1 (MQ=255)
                                  cGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTgc                                                                                                                <  1:2982003/148‑1 (MQ=255)
                                    gAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTgcgc                                                                                                              <  1:2611111/148‑1 (MQ=255)
                                     aaTACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCa                                                                                                             <  2:1274640/148‑1 (MQ=255)
                                     aaTACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCa                                                                                                             >  2:310556/1‑148 (MQ=255)
                                      aTACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCAc                                                                                                            >  2:1518072/1‑148 (MQ=255)
                                        aCCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCACCAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACgg                                                                                                          >  2:2193240/1‑148 (MQ=255)
                                        aCCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACgg                                                                                                          <  1:1957770/148‑1 (MQ=255)
                                           gTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGgatg                                                                                                       <  1:120221/148‑1 (MQ=255)
                                            ttCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGgatga                                                                                                      >  2:1146875/1‑148 (MQ=255)
                                             tCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGgatgat                                                                                                     <  1:97610/148‑1 (MQ=255)
                                             tCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGgatgat                                                                                                     >  1:2784381/1‑148 (MQ=255)
                                             tCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGgatgat                                                                                                     >  1:446851/1‑148 (MQ=255)
                                              cTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGatgatg                                                                                                    <  1:3045240/148‑1 (MQ=255)
                                              cTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGatgatg                                                                                                    <  1:2706746/148‑1 (MQ=255)
                                                gTACTGCATAAGCAACTACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGatgatgta                                                                                                  >  1:2353329/1‑146 (MQ=255)
                                                  aCTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGcagc                                                                                                <  2:248936/148‑1 (MQ=255)
                                                   cTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTCGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGcagcc                                                                                               >  2:1831834/1‑148 (MQ=255)
                                                   cTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGtagcc                                                                                               <  2:1432711/148‑1 (MQ=255)
                                                      cATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGcagccagc                                                                                            <  1:1341433/148‑1 (MQ=255)
                                                        tAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCAc                                                                                          >  2:2306925/1‑148 (MQ=255)
                                                            cAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGaa                                                                                      >  1:3200823/1‑148 (MQ=255)
                                                                 aCATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCCGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCt                                                                                 <  2:2252893/148‑1 (MQ=255)
                                                                  cATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCtt                                                                                >  1:945775/1‑148 (MQ=255)
                                                                      aCGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTcgcg                                                                            <  1:2456349/148‑1 (MQ=255)
                                                                       cGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGa                                                                           <  1:2322849/148‑1 (MQ=255)
                                                                        gCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGAt                                                                          <  2:266553/148‑1 (MQ=255)
                                                                         cAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGAtt                                                                         <  1:2101926/148‑1 (MQ=255)
                                                                              ccTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGcc                                                                    >  1:1786861/1‑148 (MQ=255)
                                                                                tgctgGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCCCGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCgg                                                                  >  2:2045762/1‑147 (MQ=255)
                                                                                   tgGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGTAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGt                                                               <  1:2566947/148‑1 (MQ=255)
                                                                                    ggTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGtt                                                              <  1:1949711/148‑1 (MQ=255)
                                                                                      tAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCa                                                            >  1:1053634/1‑148 (MQ=255)
                                                                                        gTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGa                                                          >  1:1509447/1‑148 (MQ=255)
                                                                                         tCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGAt                                                         >  1:2501199/1‑148 (MQ=255)
                                                                                            tCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATccc                                                      <  2:724410/148‑1 (MQ=255)
                                                                                             cGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCa                                                     >  1:1058819/1‑148 (MQ=255)
                                                                                              ggCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAg                                                    <  2:2085317/148‑1 (MQ=255)
                                                                                                cAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGtt                                                  >  1:970813/1‑148 (MQ=255)
                                                                                                cAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGtt                                                  >  2:3378054/1‑148 (MQ=255)
                                                                                                 aaTTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTg                                                 <  2:2390390/148‑1 (MQ=255)
                                                                                                  aTTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGt                                                >  1:2616855/1‑148 (MQ=255)
                                                                                                   tttGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTa                                               <  1:2011788/148‑1 (MQ=255)
                                                                                                      gCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTActc                                            <  1:1448472/148‑1 (MQ=255)
                                                                                                      gCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTActc                                            <  1:764235/148‑1 (MQ=255)
                                                                                                      gCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTActc                                            <  1:2465021/148‑1 (MQ=255)
                                                                                                       cTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGTAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTActct                                           >  1:1691212/1‑148 (MQ=255)
                                                                                                         tGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGTAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTc                                         >  1:1868235/1‑148 (MQ=255)
                                                                                                         tGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTc                                         >  1:3147564/1‑148 (MQ=255)
                                                                                                            aGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCTCGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGa                                      <  2:457116/148‑1 (MQ=255)
                                                                                                             gTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGac                                     <  1:3248770/148‑1 (MQ=255)
                                                                                                              tACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGTAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGacg                                    >  2:345492/1‑148 (MQ=255)
                                                                                                              tACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGacg                                    >  1:2377794/1‑148 (MQ=255)
                                                                                                                     ccAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGacgcagcacg                             >  1:2198325/1‑148 (MQ=255)
                                                                                                                      cAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGacgcagcacgc                            <  1:1670962/148‑1 (MQ=255)
                                                                                                                        gCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGacgcagcacgcag                          <  1:1820365/148‑1 (MQ=255)
                                                                                                                        gCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGacgcagcacgcag                          <  2:19237/148‑1 (MQ=255)
                                                                                                                         cAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGacgcagcacg                             >  1:636828/1‑144 (MQ=255)
                                                                                                                         cAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGacgcagcacg                             <  2:636828/144‑1 (MQ=255)
                                                                                                                         cAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGacgcagcacgcagc                         <  2:1114263/148‑1 (MQ=255)
                                                                                                                            cAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTga                      >  1:46421/1‑148 (MQ=255)
                                                                                                                             aGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTgag                     <  2:674099/148‑1 (MQ=255)
                                                                                                                             aGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTgag                     <  1:332884/148‑1 (MQ=255)
                                                                                                                               gTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTgagag                   <  2:1729667/148‑1 (MQ=255)
                                                                                                                                 tAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGaa                 <  1:653227/148‑1 (MQ=255)
                                                                                                                                  aGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGTAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAg                <  1:504781/148‑1 (MQ=255)
                                                                                                                                   gCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGTAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGc               >  2:1004470/1‑148 (MQ=255)
                                                                                                                                   gCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGc               >  2:1787842/1‑148 (MQ=255)
                                                                                                                                    cGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGTAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGcc              <  1:3078822/148‑1 (MQ=255)
                                                                                                                                     gATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGTAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGccc             <  1:3385959/148‑1 (MQ=255)
                                                                                                                                            ggATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGt      <  2:824348/148‑1 (MQ=255)
                                                                                                                                             gATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTa     <  2:282662/148‑1 (MQ=255)
                                                                                                                                             gATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTa     <  2:2085409/148‑1 (MQ=255)
                                                                                                                                                tttAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGa  <  1:981710/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
CGTAATAGGCAACCGCTGCGGAGAAGGTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGA  >  CP009273/2093465‑2093756

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