Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,913,019 C→A 88.4% A363S (GCT→TCT)  gudX ← glucarate dehydratase‑related protein, substrate unknown

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,913,0190CA88.4% 419.1 / 55.1 173A363S (GCT→TCT) gudXglucarate dehydratase‑related protein, substrate unknown
Reads supporting (aligned to +/- strand):  ref base C (8/12);  new base A (82/71);  total (90/83)
Fisher's exact test for biased strand distribution p-value = 3.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.08e-01

Only 100 of 173 total aligned reads displayed.

TCATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGCGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAGTCGTCGCAAAGCTGCGCCACACGGACTGCACCGGAAAGCGTCCAGAA  >  CP009273/2912874‑2913163
                                                                                                                                                 |                                                                                                                                                
 cATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGt                                                                                                                                               <  1:2773637/148‑1 (MQ=255)
 cATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATAGATAGAGGt                                                                                                                                               >  1:2424990/1‑148 (MQ=255)
     tgCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGa                                                                                                                                           <  1:1122005/148‑1 (MQ=255)
      gCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGAt                                                                                                                                          <  1:255192/148‑1 (MQ=255)
       ccTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATAGATAGAGGTAGGAtt                                                                                                                                         >  1:568995/1‑148 (MQ=255)
          tttGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGCGGTAGGATTAcc                                                                                                                                      >  2:105373/1‑148 (MQ=255)
          tttGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTGCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATAGATAGAGGTAGGATTAcc                                                                                                                                      <  2:1140422/148‑1 (MQ=255)
              tACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATAGATAGAGGTAGGATTACCCGGt                                                                                                                                  >  1:2918357/1‑148 (MQ=255)
               aCCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTg                                                                                                                                 >  2:2584984/1‑148 (MQ=255)
               aCCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTg                                                                                                                                 >  2:1304096/1‑148 (MQ=255)
                ccTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTgc                                                                                                                                >  2:467189/1‑148 (MQ=255)
                ccTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTgc                                                                                                                                >  2:3009825/1‑148 (MQ=255)
                ccTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTgc                                                                                                                                >  2:469387/1‑148 (MQ=255)
                   gTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGGATAGAGAGAGGTAGGATTACCCGGTgccgc                                                                                                                             >  2:1880405/1‑148 (MQ=255)
                     tCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGccgccg                                                                                                                           >  2:2886980/1‑148 (MQ=255)
                       ccagtccagtTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGTCGACAATCGCCCTCCTGCCAAATCAAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCcgcg                                                                                                                         >  2:948857/1‑148 (MQ=255)
                         agtccagtTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGcc                                                                                                                       >  1:293238/1‑148 (MQ=255)
                            ccagtTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCcac                                                                                                                    >  2:2396555/1‑148 (MQ=255)
                            ccagtTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCcac                                                                                                                    <  1:2521897/148‑1 (MQ=255)
                             cagtTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCcaca                                                                                                                   <  1:965621/148‑1 (MQ=255)
                                ttCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATgg                                                                                                                >  1:862021/1‑148 (MQ=255)
                                ttCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATgg                                                                                                                >  1:711583/1‑148 (MQ=255)
                                ttCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATgg                                                                                                                >  1:3091071/1‑148 (MQ=255)
                                ttCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATAGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATgg                                                                                                                >  1:2939126/1‑148 (MQ=255)
                                    aCGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTaa                                                                                                            <  1:2234337/148‑1 (MQ=255)
                                    aCGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATAGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTaa                                                                                                            <  1:535245/148‑1 (MQ=255)
                                     cGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTaaa                                                                                                           >  1:737538/1‑148 (MQ=255)
                                     cGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATAGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTaaa                                                                                                           <  2:21489/148‑1 (MQ=255)
                                      gCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAAc                                                                                                          <  2:1635127/148‑1 (MQ=255)
                                        ccAGACCGTGCGCATCAGGATAGGCAATTTTTGCGTTTTTAATCTCCAGCTGATTTTGGGTCAGGCGACAATCGCCCTCGTGCCAAATCGAGTGGGTATCGATAGAGGTAGGATTACCCGTTGCCGCCGCGCCCACTTGGGTATActt                                                                                                        <  2:3096560/148‑3 (MQ=255)
                                        ccAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGCGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACAt                                                                                                        <  1:1976470/148‑1 (MQ=255)
                                        ccAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACAt                                                                                                        <  1:2271211/148‑1 (MQ=255)
                                          aGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATAGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCg                                                                                                      <  1:2217465/148‑1 (MQ=255)
                                            aCCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGcc                                                                                                    >  1:2199921/1‑148 (MQ=255)
                                                ggCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCagag                                                                                                <  1:1695027/148‑1 (MQ=255)
                                                ggCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCagag                                                                                                <  1:1695061/148‑1 (MQ=255)
                                                ggCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCagag                                                                                                >  2:2677617/1‑148 (MQ=255)
                                                    cATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGat                                                                                            <  2:2767461/148‑1 (MQ=255)
                                                    cATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGat                                                                                            <  1:2000611/148‑1 (MQ=255)
                                                     aTCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGata                                                                                           >  1:1805240/1‑148 (MQ=255)
                                                         ggAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGa                                                                                       >  2:1186387/1‑148 (MQ=255)
                                                         ggAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATAGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGa                                                                                       >  2:641929/1‑148 (MQ=255)
                                                          gAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGaa                                                                                      >  2:305014/1‑148 (MQ=255)
                                                          gAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGaa                                                                                      >  1:272473/1‑148 (MQ=255)
                                                             cgtcAATTTGTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGTCGACAATCGCGCTCCTGCAAAATCCTGTTGGTATCGATAGCGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATg                                                                                   <  1:2591042/145‑1 (MQ=38)
                                                                cAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGCGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGtt                                                                                <  2:545973/148‑1 (MQ=255)
                                                                cAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGtt                                                                                <  1:428415/148‑1 (MQ=255)
                                                                cAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGtt                                                                                <  1:178430/148‑1 (MQ=255)
                                                                 aaTTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGtta                                                                               <  1:445028/148‑1 (MQ=255)
                                                                 aaTTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGtta                                                                               <  1:2002082/148‑1 (MQ=255)
                                                                 aaTTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATAGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGtta                                                                               <  2:3261479/148‑1 (MQ=255)
                                                                  aTTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATAGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGttat                                                                              <  2:2606250/148‑1 (MQ=255)
                                                                   tttttCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGttatt                                                                             <  2:2097108/148‑1 (MQ=255)
                                                                   tttttCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGttatt                                                                             <  2:144349/148‑1 (MQ=255)
                                                                    ttttCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGCGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGttatta                                                                            >  1:2995295/1‑148 (MQ=255)
                                                                     tttCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGCGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAg                                                                           <  1:1231976/148‑1 (MQ=255)
                                                                      ttCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGa                                                                          <  1:2098076/148‑1 (MQ=255)
                                                                       tCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATAGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGaa                                                                         <  2:1690186/148‑1 (MQ=255)
                                                                        ccGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAAt                                                                        >  1:2713749/1‑148 (MQ=255)
                                                                          gTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATgg                                                                      <  1:2440768/148‑1 (MQ=255)
                                                                          gTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATgg                                                                      <  2:2623584/148‑1 (MQ=255)
                                                                          gTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATgg                                                                      >  2:2186832/1‑148 (MQ=255)
                                                                           tttttAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGc                                                                     <  1:511805/148‑1 (MQ=255)
                                                                                 aTCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGcccc                                                               >  2:1840162/1‑148 (MQ=255)
                                                                                 aTCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGcccc                                                               >  2:1924711/1‑148 (MQ=255)
                                                                                   ctcCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAg                                                             >  1:1736705/1‑148 (MQ=255)
                                                                                   ctcCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAg                                                             >  1:1738248/1‑148 (MQ=255)
                                                                                      cAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTc                                                          >  2:923167/1‑148 (MQ=255)
                                                                                        gCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAg                                                        >  1:2001956/1‑148 (MQ=255)
                                                                                          ggATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGcc                                                      <  1:1764694/148‑1 (MQ=255)
                                                                                            aTTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGcccc                                                    <  2:293238/148‑1 (MQ=255)
                                                                                            aTTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATAGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGcccc                                                    <  2:2016845/148‑1 (MQ=255)
                                                                                                tGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGCGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAgt                                                <  1:307245/148‑1 (MQ=255)
                                                                                                  ggTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGCGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAgtcg                                              <  2:2016340/148‑1 (MQ=255)
                                                                                                    tCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAgtcgtc                                            >  2:2908269/1‑148 (MQ=255)
                                                                                                           aCAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAGTCGTCGCAAAGc                                     <  1:929504/148‑1 (MQ=255)
                                                                                                           aCAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAGTCGTCGCAAAGc                                     <  1:2584984/148‑1 (MQ=255)
                                                                                                            cAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAGTCGTCGCAAAGCt                                    >  2:773946/1‑148 (MQ=255)
                                                                                                                cGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAGTCGTCGCAAAGCTgcgc                                <  2:2726192/148‑1 (MQ=255)
                                                                                                                  cccTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAGTCGTCGCAAAGCTGCGCca                              >  1:1502624/1‑148 (MQ=255)
                                                                                                                   cctcctGCCAAATCCAGTGGGTATCGATAGCGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAGTCGTCGCAAAGCTGCGCcac                             >  1:1226229/1‑148 (MQ=255)
                                                                                                                        tGCCAAATCCAGTGGGTATCGATAGCGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAGTCGTCGCAAAGCTGCGCCACACGGa                        >  1:1928712/1‑148 (MQ=255)
                                                                                                                        tGCCAAATCCAGTGGGTATCGATAGCGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAGTCGTCGCAAAGCTGCGCCACACGGa                        >  1:258478/1‑148 (MQ=255)
                                                                                                                        tGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAGTCGTCGCAAAGCTGCGCCACACGGa                        <  1:3023252/148‑1 (MQ=255)
                                                                                                                        tGCCAAATCCAGTGGGTATAGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAGTCGTCGCAAAGCTGCGCCACACGGa                        <  2:2424990/148‑1 (MQ=255)
                                                                                                                         gCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAGTCGTCGCAAAGCTGCGCCACACGGAc                       <  1:1604734/148‑1 (MQ=255)
                                                                                                                            aaaTCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAGTCGTCGCAAAGCTGCGCCACACGGACTGc                    >  2:1738312/1‑148 (MQ=255)
                                                                                                                                   gTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAGTCGTCGCAAAGCTGCGCCACACGGACTGCACCGGaa             <  2:1683902/148‑1 (MQ=255)
                                                                                                                                      ggTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAGTCGTCGCAAAGCTGCGCCACACGGACTGCACCGGAAAGc          >  2:2423682/1‑148 (MQ=255)
                                                                                                                                       gTATCGATAGCGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAGTCGTCGCAAAGCTGCGCCACACGGACTGCACCGGAAAGCg         <  1:105373/148‑1 (MQ=255)
                                                                                                                                         aTCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAGTCGTCGCAAAGCTGCGCCACACGGACTGCACCGGACAGCGTc       >  2:2746842/1‑148 (MQ=255)
                                                                                                                                         aTCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAGTCGTCGCAAAGCTGCGCCACACGGACTGCACCGGAAAGCGTc       >  2:983896/1‑148 (MQ=255)
                                                                                                                                           cGATAGCGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAGTCGTCGCAAAGCTGCGCCACACGGACTGCACCGGAAAGCGTCCa     >  2:2416927/1‑148 (MQ=255)
                                                                                                                                           cGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAGTCGTCGCAAAGCTGCGCCACACGGACTGCACCGGAAAGCGTCCa     >  2:336522/1‑148 (MQ=255)
                                                                                                                                           agATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAGTCGTCGCAAAGCTGCGCCACACGGACTGCACCGGAAAGCGTCCa     >  2:2219972/2‑148 (MQ=255)
                                                                                                                                            gATAGCGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAGTCGTCGCAAAGCTGCGCCACACGGACTGCACCGGAAAGCGTCCAg    <  2:1559795/148‑1 (MQ=255)
                                                                                                                                            gATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAGTCGTCGCAAAGCTGCGCCACACGGACTGCACCGGAAAGCGTCCAg    >  2:3095379/1‑148 (MQ=255)
                                                                                                                                            gATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAGTCGTCGCAAAGCTGCGCCACACGGACTGCACCGGAAAGCGTCCAg    >  1:1940518/1‑148 (MQ=255)
                                                                                                                                             aTAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAGTCGTCGCAAAGCTGCGCCACACGGACTGCACCGGAAAGCGTCCAGa   >  1:105954/1‑148 (MQ=255)
                                                                                                                                              tAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAGTCGTCGCAAAGCTGCGCCACACGGACTGCACCGGAAAGCGTCCAGaa  <  1:2643043/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
TCATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGCGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAGTCGTCGCAAAGCTGCGCCACACGGACTGCACCGGAAAGCGTCCAGAA  >  CP009273/2912874‑2913163

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: