Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,913,019 C→A 100% 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,0190CA96.2% 547.0 / 12.0 183A363S (GCT→TCT) gudXglucarate dehydratase‑related protein, substrate unknown
Reads supporting (aligned to +/- strand):  ref base C (4/3);  new base A (84/92);  total (88/95)
Fisher's exact test for biased strand distribution p-value = 7.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 183 total aligned reads displayed.

TCATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGCGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAGTCGTCGCAAAGCTGCGCCACACGGACTGCACCGGAAAGCGTCCAGAA  >  CP009273/2912874‑2913163
                                                                                                                                                 |                                                                                                                                                
tCATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAgg                                                                                                                                                <  1:76549/148‑1 (MQ=255)
 cATGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGt                                                                                                                                               >  2:1103017/1‑148 (MQ=255)
  aTGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTa                                                                                                                                              <  1:272334/148‑1 (MQ=255)
  aTGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTa                                                                                                                                              <  1:303027/148‑1 (MQ=255)
  aTGTGCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTa                                                                                                                                              <  2:176293/148‑1 (MQ=255)
    gtgCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAgg                                                                                                                                            <  2:1987601/148‑1 (MQ=255)
    gtgCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAgg                                                                                                                                            <  2:1527039/148‑1 (MQ=255)
    gtgCCTTTTGTACCTGTTCCCAGTCCAGGTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAgg                                                                                                                                            <  1:699429/148‑1 (MQ=255)
     tgCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGa                                                                                                                                           >  1:2059449/1‑148 (MQ=255)
      gCCTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGAt                                                                                                                                          <  2:1567318/148‑1 (MQ=255)
        cTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTa                                                                                                                                        >  2:550994/1‑148 (MQ=255)
        cTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTa                                                                                                                                        <  1:2858947/148‑1 (MQ=255)
        cTTTTGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTa                                                                                                                                        <  1:1078052/148‑1 (MQ=255)
         ttttGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTAc                                                                                                                                       <  1:2505538/148‑1 (MQ=255)
         ttttGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTAc                                                                                                                                       >  1:2257255/1‑148 (MQ=255)
         ttttGTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTAc                                                                                                                                       <  1:2505564/148‑1 (MQ=255)
             gTACCTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCgg                                                                                                                                   <  1:1185360/148‑1 (MQ=255)
                ccTGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGCGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGCGGTAGGATTACCCGGTgc                                                                                                                                >  1:279508/1‑148 (MQ=255)
                  tGTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATAGATAGAGGTAGGATTACCCGGTgccg                                                                                                                              >  2:2262768/1‑148 (MQ=255)
                   gTTCCCAGTCCAGTTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTgccgc                                                                                                                             <  1:1338915/148‑1 (MQ=255)
                        cagtccagtTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTCTCGATAGAGGTAGGATTACCCGGTGCCGCCgcgc                                                                                                                        >  2:2587407/1‑148 (MQ=255)
                        cagtccagtTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCgcgc                                                                                                                        >  2:2587357/1‑148 (MQ=255)
                             cagtTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCcaca                                                                                                                   <  1:2323789/148‑1 (MQ=255)
                              agtTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACAt                                                                                                                  >  2:2999383/1‑148 (MQ=255)
                              agtTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACAt                                                                                                                  >  1:3353760/1‑148 (MQ=255)
                              agtTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACAt                                                                                                                  >  2:3001145/1‑148 (MQ=255)
                               gtTCCACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATg                                                                                                                 >  1:2589638/1‑148 (MQ=255)
                                   cACGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTa                                                                                                             >  1:1599706/1‑148 (MQ=255)
                                    aCGCCCAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTaa                                                                                                            <  2:3116937/148‑1 (MQ=255)
                                       cccAGACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACa                                                                                                         <  1:2570521/148‑1 (MQ=255)
                                           gACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTCTCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGc                                                                                                     >  2:1966788/1‑148 (MQ=255)
                                           gACCGGGCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGc                                                                                                     >  2:3432894/1‑148 (MQ=255)
                                                ggCGCATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCagag                                                                                                >  2:580575/1‑148 (MQ=255)
                                                  cgcATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCagagag                                                                                              <  1:2661189/148‑1 (MQ=255)
                                                  cgcATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCagagag                                                                                              >  1:2976485/1‑148 (MQ=255)
                                                  cgcATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCagagag                                                                                              <  2:229654/148‑1 (MQ=255)
                                                  cgcATCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCagagag                                                                                              <  1:3414335/148‑1 (MQ=255)
                                                     aTCAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGata                                                                                           <  2:3507233/148‑1 (MQ=255)
                                                       cAGGAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATc                                                                                         >  2:718935/1‑148 (MQ=255)
                                                         ggAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGa                                                                                       <  1:3039185/148‑1 (MQ=255)
                                                         ggAACGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGa                                                                                       >  2:3164133/1‑148 (MQ=255)
                                                           aaCGGCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGaaa                                                                                     <  1:1495309/148‑1 (MQ=255)
                                                              ggCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATgg                                                                                  <  1:1747753/148‑1 (MQ=255)
                                                               gCAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGt                                                                                 <  2:1040341/148‑1 (MQ=255)
                                                                cAATTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGtt                                                                                >  1:2928029/1‑148 (MQ=255)
                                                                  aTTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGttat                                                                              >  2:3270708/1‑148 (MQ=255)
                                                                  aTTTTTCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGttat                                                                              >  2:3270575/1‑148 (MQ=255)
                                                                   tttttCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGttatt                                                                             <  2:779510/148‑1 (MQ=255)
                                                                      ttCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGa                                                                          <  2:3188956/148‑1 (MQ=255)
                                                                       tCCGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGaa                                                                         <  1:1879701/148‑1 (MQ=255)
                                                                        ccGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAAt                                                                        >  1:2056885/1‑148 (MQ=255)
                                                                        ccGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATAGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAAt                                                                        <  1:2072973/148‑1 (MQ=255)
                                                                        ccGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATAGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAAt                                                                        <  1:1938217/148‑1 (MQ=255)
                                                                         cGTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATg                                                                       <  2:2095878/148‑1 (MQ=255)
                                                                          gTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATgg                                                                      >  1:706055/1‑148 (MQ=255)
                                                                          gTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATgg                                                                      <  2:2176831/148‑1 (MQ=255)
                                                                          gTTTTTAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATgg                                                                      >  1:706051/1‑148 (MQ=255)
                                                                           tttttAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGc                                                                     <  2:459776/148‑1 (MQ=255)
                                                                            ttttAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCa                                                                    <  2:659398/148‑1 (MQ=255)
                                                                            ttttAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCa                                                                    <  2:659418/148‑1 (MQ=255)
                                                                             tttAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAg                                                                   >  2:454999/1‑148 (MQ=255)
                                                                               tAATCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGcc                                                                 <  2:3475506/148‑1 (MQ=255)
                                                                                aaTCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGccc                                                                <  2:1349691/148‑1 (MQ=255)
                                                                                 aTCTCCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGcccc                                                               >  1:1547509/1‑148 (MQ=255)
                                                                                  tctcCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCa                                                              >  1:1865957/1‑148 (MQ=255)
                                                                                  tctcCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCa                                                              <  1:100835/148‑1 (MQ=255)
                                                                                  tctcCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCAGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCa                                                              >  2:286751/1‑148 (MQ=255)
                                                                                   ctcCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAg                                                             >  2:1937421/1‑148 (MQ=255)
                                                                                    tcCAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAgg                                                            >  2:2776734/1‑148 (MQ=255)
                                                                                     ccAGCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGt                                                           <  2:3418523/148‑1 (MQ=255)
                                                                                        gCGGATTTTGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAg                                                        <  1:3141475/148‑1 (MQ=255)
                                                                                              tttGGGTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCa                                                  >  1:1967027/1‑148 (MQ=255)
                                                                                                 gggTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAgtc                                               <  1:1633860/148‑1 (MQ=255)
                                                                                                 gggTCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAgtc                                               <  1:257350/148‑1 (MQ=255)
                                                                                                    tCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGCGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAgtcgtc                                            >  2:1013915/1‑148 (MQ=255)
                                                                                                    tCAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGCGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAgtcgtc                                            >  2:2794798/1‑148 (MQ=255)
                                                                                                     cAGGCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAGtcgtcg                                           <  1:812524/148‑1 (MQ=255)
                                                                                                       ggCGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAGTCGTCGCa                                         >  1:2821905/1‑148 (MQ=255)
                                                                                                         cGACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAGTCGTCGCaaa                                       >  1:3080346/1‑148 (MQ=255)
                                                                                                          gACAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAGTCGTCGCAAAg                                      <  2:568297/148‑1 (MQ=255)
                                                                                                            cAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAGTCGTCGCAAAGCt                                    <  2:1850964/148‑1 (MQ=255)
                                                                                                            cAATCGCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAGTCGTCGCAAAGCt                                    <  2:632911/148‑1 (MQ=255)
                                                                                                                 gCCCTCCTGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAGTCGTCGCAAAGCTGCGcc                               <  1:867998/148‑1 (MQ=255)
                                                                                                                     tcctGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAGTCGTCGCAAAGCTGCGCCacac                           >  1:592876/1‑148 (MQ=255)
                                                                                                                       ctGCCAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCTAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAGTCGTCGCAAAGCTGCGCCACACgg                         <  2:1599545/148‑1 (MQ=255)
                                                                                                                        tGCCAAATCCAGTGGGTATAGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAGTCGTCGCAAAGCTGCGCCACACGGa                        >  2:945140/1‑148 (MQ=255)
                                                                                                                          ccAAATCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAGTCGTCGCAAAGCTGCGCCACACGGACt                      >  1:1865330/1‑148 (MQ=255)
                                                                                                                              aTCCAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAGTCGTCGCAAAGCTGCGCCACACGGACTGCAc                  >  1:1480031/1‑148 (MQ=255)
                                                                                                                                 cAGTGGGTATCGATAGCGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAGTCGTCGCAAAGCTGCGCCACACGGACTGCACCgg               <  1:3363624/148‑1 (MQ=255)
                                                                                                                                 cAGTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAGTCGTCGCAAAGCTGCGCCACACGGACTGCACCgg               <  2:2589638/148‑1 (MQ=255)
                                                                                                                                   gTGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAGTCGTCGCAAAGCTGCGCCACACGGACTGCACCGGaa             <  1:2454177/148‑1 (MQ=255)
                                                                                                                                    tGGGTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAGTCGTCGCAAAGCTGCGCCACACGGACTGCACCGGaaa            <  1:455591/148‑1 (MQ=255)
                                                                                                                                     gggTATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAGTCGTCGCAAAGCTGCGCCACACGGACTGCACCGGAAAg           <  2:1600461/148‑1 (MQ=255)
                                                                                                                                        tATCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAGTCGTCGCAAAGCTGCGCCACACGGACTGCACCGGAAAGCGt        >  2:3119971/1‑148 (MQ=255)
                                                                                                                                         atagATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCATGTCAGCCCCCAGTCGTCGCAAAGCTGCGCCACACGGACTGCACCGGAAAGCGTc       >  1:529593/4‑148 (MQ=255)
                                                                                                                                         aTCGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAGTCGTCGCAAAGCTGCGCCACACGGACTGCACCGGAAAGCGTc       <  1:2015459/148‑1 (MQ=255)
                                                                                                                                           cGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAGTCGTCGCAAAGCTGCGCCACACGGACTGCACCGGAAAGCGTCCa     >  1:367282/1‑148 (MQ=255)
                                                                                                                                           cGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAGTCGTCGCAAAGCTGCGCCACACGGACTGCACCGGAAAGCGTCCa     >  2:1224947/1‑148 (MQ=255)
                                                                                                                                           cGATAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAGTCGTCGCAAAGCTGCGCCACACGGACTGCACCGGAAAGCGTCCa     >  1:1348427/1‑148 (MQ=255)
                                                                                                                                             aTAGAGGTAGGATTACCCGGTGCCGCCGCGCCCACATGGGTAAACATCGCCAGAGAGATATCGAAATGGTTATTAGAATGGCAGCCCCAGGTCAGCCCCCAGTCGTCGCAAAGCTGCGCCACACGGACTGCACCGGAAAGCGTCCAGa   >  2:874362/1‑148 (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: