Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 784,094 G→T 5.2% A67E (GCG→GAG)  galM ← aldose 1‑epimerase; type‑1 mutarotase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273784,0940GT 5.2% 1161.5 / 24.4 300A67E (GCG→GAG) galMaldose 1‑epimerase; type‑1 mutarotase
Reads supporting (aligned to +/- strand):  ref base G (150/134);  new base T (11/4);  total (161/139)
Fisher's exact test for biased strand distribution p-value = 1.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

Only 100 of 304 total aligned reads displayed.

TGCCAGCGACGTTTGTCGAACCCTTCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCGGAAAGTAAAGTCGCACCCCAGTCCATCAGCGTGACTACCATCCCTGCGTTGTTACGCAAAGTTAAC  >  CP009273/783947‑784243
                                                                                                                                                   |                                                                                                                                                     
 gCCAGCGACGTTTGTCGAACCCTTCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCgg                                                                                                                                                    <  2:205635/150‑1 (MQ=255)
   cAGCGACGTTTGTCGAACCCTTCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGcc                                                                                                                                                  <  1:3803449/150‑1 (MQ=255)
   cAGCGACGTTTGTCGAACCCTTCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGcc                                                                                                                                                  <  1:529134/150‑1 (MQ=255)
   cAGCGACGTTTGTCGAACCCTTCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGcc                                                                                                                                                  >  1:843469/1‑150 (MQ=255)
     gCGACGTTTGTCGAACCCTTCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGcctg                                                                                                                                                <  1:1448214/150‑1 (MQ=255)
      cGACGTTTGTCGAACCCTTCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGcctga                                                                                                                                               >  1:1467607/1‑150 (MQ=255)
       gACGTTTGTCGAACCCTTCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGcctgat                                                                                                                                              <  1:1640188/150‑1 (MQ=255)
          gTTTGTCGAACCCTTCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGcctgatcct                                                                                                                                           >  2:4350467/1‑150 (MQ=255)
          gTTTGTCGAACCCTTCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAAGAGAGGCCCCCAGAAACTCGGcctgatcct                                                                                                                                           <  1:1711000/150‑1 (MQ=255)
            ttGTCGAACCCTTCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGcctgatcctga                                                                                                                                         >  1:661407/1‑150 (MQ=255)
            ttGTCGAACCCTTCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGcctgatcctga                                                                                                                                         >  2:306508/1‑150 (MQ=255)
            ttGTCGAACCCTTCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGcctgatcctga                                                                                                                                         >  2:3047488/1‑150 (MQ=255)
             tGTCGAACCCTTCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGcctgatcctgat                                                                                                                                        <  1:1607261/150‑1 (MQ=255)
             tGTCGAACCCTTCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGcctgatcctgat                                                                                                                                        >  1:3565664/1‑150 (MQ=255)
             tGTCGAACCCTTCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGcctgatcctgat                                                                                                                                        >  1:3000148/1‑150 (MQ=255)
                 gAACCCTTCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCa                                                                                                                                    >  2:4571507/1‑150 (MQ=255)
                     ccTTCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTcc                                                                                                                                <  1:1728882/150‑1 (MQ=255)
                     ccTTCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTcc                                                                                                                                >  2:1754185/1‑150 (MQ=255)
                         ccGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGc                                                                                                                            >  2:2325602/1‑150 (MQ=255)
                         ccGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGc                                                                                                                            <  1:1288687/150‑1 (MQ=255)
                           ggCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTg                                                                                                                          <  1:1616236/150‑1 (MQ=255)
                           ggCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTg                                                                                                                          >  1:3742373/1‑150 (MQ=255)
                             cccGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTGGc                                                                                                                        >  2:3117527/1‑150 (MQ=255)
                               cgccgTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTGGcac                                                                                                                      <  2:4241500/150‑1 (MQ=255)
                                gccgTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTGGcaca                                                                                                                     <  2:1381413/150‑1 (MQ=255)
                                gccgTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTGGcaca                                                                                                                     >  1:1509057/1‑150 (MQ=255)
                                       aGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGagc                                                                                                              <  2:703225/150‑1 (MQ=255)
                                        ggtgtttAACGCCATTACTTGGCGAAAGCGTCTCGGTTTCACCGTAAAAGTTATAACGTCTATTGGCTATAGGGTTGGAATAACTAAAAATAGAGGCCCCCATAAACGCGGCATGATCCTGATAGCATTCCGGGCTGGAACAGCCGagca                                                                                                             <  2:1677966/145‑1 (MQ=255)
                                          tGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGagcagc                                                                                                           >  2:4126268/1‑150 (MQ=255)
                                          tGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGagcagc                                                                                                           >  1:2308875/1‑150 (MQ=255)
                                           ggTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAgcg                                                                                                          <  1:1657966/150‑1 (MQ=255)
                                             ttAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGcc                                                                                                        >  2:1906293/1‑150 (MQ=255)
                                              tAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCt                                                                                                       >  1:4349046/1‑150 (MQ=255)
                                               aaCGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCATAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTc                                                                                                      <  1:4217807/150‑1 (MQ=255)
                                               aaCGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTc                                                                                                      <  2:3242995/150‑1 (MQ=255)
                                                aCGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTcg                                                                                                     <  1:2670418/150‑1 (MQ=255)
                                                  gCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTcgcg                                                                                                   <  1:415675/150‑1 (MQ=255)
                                                   cccTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCgg                                                                                                  >  1:851421/1‑150 (MQ=255)
                                                     cTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGAc                                                                                                <  1:232754/150‑1 (MQ=255)
                                                        aCTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCt                                                                                             <  2:3572654/150‑1 (MQ=255)
                                                         cTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTg                                                                                            <  1:2952578/150‑1 (MQ=255)
                                                             gCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCAt                                                                                        <  1:3113679/150‑1 (MQ=255)
                                                              cGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATc                                                                                       <  1:3528339/150‑1 (MQ=255)
                                                               gAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCg                                                                                      >  2:3721299/1‑150 (MQ=255)
                                                               gAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCg                                                                                      >  2:4072953/1‑150 (MQ=255)
                                                               gAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCg                                                                                      >  1:1805254/1‑150 (MQ=255)
                                                               gAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCg                                                                                      >  2:4534836/1‑150 (MQ=255)
                                                                 aaGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGa                                                                                    >  2:1400933/1‑150 (MQ=255)
                                                                  aGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGaa                                                                                   >  1:2602342/1‑150 (MQ=255)
                                                                     gTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGc                                                                                <  1:68877/150‑1 (MQ=255)
                                                                       cACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAATCGCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCgg                                                                              <  2:3947615/150‑1 (MQ=255)
                                                                       cACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCgg                                                                              >  2:3452925/1‑150 (MQ=255)
                                                                        aCGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGa                                                                             >  1:1339250/1‑150 (MQ=255)
                                                                          ggTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAAt                                                                           <  1:1427615/150‑1 (MQ=255)
                                                                           gTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATa                                                                          >  1:995767/1‑150 (MQ=255)
                                                                            tttCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATAc                                                                         >  1:190186/1‑150 (MQ=255)
                                                                                 ccGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCg                                                                    <  1:4219982/150‑1 (MQ=255)
                                                                                  cGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCgg                                                                   <  1:3049081/150‑1 (MQ=255)
                                                                                     cAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCGGaaa                                                                <  2:2159791/150‑1 (MQ=255)
                                                                                     cAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCGGaaa                                                                <  2:2531527/150‑1 (MQ=255)
                                                                                      aaaGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCGGaaag                                                               >  2:3639543/1‑150 (MQ=255)
                                                                                      aaaGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCGGaaag                                                               >  2:3521471/1‑150 (MQ=255)
                                                                                      aaaGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCGGaaag                                                               <  1:2116033/150‑1 (MQ=255)
                                                                                       aaGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCGGaaagt                                                              <  2:4039635/150‑1 (MQ=255)
                                                                                         ggTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACTCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCGGaaagtaa                                                            >  1:2437890/1‑150 (MQ=255)
                                                                                         ggTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCGGaaagtaa                                                            <  1:4571510/150‑1 (MQ=255)
                                                                                          gTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCGGaaagtaaa                                                           >  1:3963923/1‑150 (MQ=255)
                                                                                            ataACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCGGaaagtaaagt                                                         <  2:515189/150‑1 (MQ=255)
                                                                                            ataACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCGGaaagtaaagt                                                         <  2:3677895/150‑1 (MQ=255)
                                                                                                cGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCGGAAAGTAAAGTCGCa                                                     <  2:1056543/150‑1 (MQ=255)
                                                                                                   cTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCGGAAAGTAAAGTCGCAccc                                                  >  1:726004/1‑150 (MQ=255)
                                                                                                   cTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCGGAAAGTAAAGTCGCAccc                                                  >  2:3300317/1‑150 (MQ=255)
                                                                                                     aTTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCGGAAAGTAAAGTCGCACCCCa                                                <  2:1851069/150‑1 (MQ=255)
                                                                                                       tGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCGGAAAGTAAAGTCGCACCCCAGt                                              <  1:4033638/150‑1 (MQ=255)
                                                                                                           gATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCGGAAAGTAAAGTCGCACCCCAGTCCAt                                          <  2:655506/150‑1 (MQ=255)
                                                                                                            aTACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCGGAAAGTAAAGTCGCACCCCAGTCCATc                                         <  1:4265213/150‑1 (MQ=255)
                                                                                                             tACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCGGAAAGTAAAGTCGCACCCCAGTCCATCa                                        <  1:4319066/150‑1 (MQ=255)
                                                                                                               cGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACTCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCGGAAAGTAAAGTCGCACCCCAGTCCATCAGc                                      >  2:2388075/1‑150 (MQ=255)
                                                                                                                ggTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCGGAAAGTAAAGTCGCACCCCAGTCCATCAGCg                                     >  2:1718636/1‑150 (MQ=255)
                                                                                                                  ttGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCGGAAAGTAAAGTCGCACCCCAGTCCATCAGCGTg                                   <  2:5855/150‑1 (MQ=255)
                                                                                                                      cATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCGGAAAGTAAAGTCGCACCCCAGTCCATCAGCGTGACTa                               <  2:2774004/150‑1 (MQ=255)
                                                                                                                       aTAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCGGAAAGTAAAGTCGCACCCCAGTCCATCAGCGTGACTAc                              >  2:1412219/1‑150 (MQ=255)
                                                                                                                          aCGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCGGAAAGTAAAGTCGCACCCCAGTCCATCAGCGTGACTACCAt                           <  1:3903041/150‑1 (MQ=255)
                                                                                                                            gACCAATAGAGGCCCCCAGAAACTCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCGGAAAGTAAAGTCGCACCCCAGTCCATCAGCGTGACTACCATcc                         >  2:350836/1‑150 (MQ=255)
                                                                                                                              ccAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCGGAAAGTAAAGTCGCACCCCAGTCCATCAGCGTGACTACCATCCCt                       <  1:4126268/150‑1 (MQ=255)
                                                                                                                              ccAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCGGAAAGTAAAGTCGCACCCCAGTCCATCAGCGTGACTACCATCCCt                       <  2:2968000/150‑1 (MQ=255)
                                                                                                                              ccAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCGGAAAGTAAAGTCGCACCCCAGTCCATCAGCGTGACTACCATCCCt                       <  1:2143946/150‑1 (MQ=255)
                                                                                                                               cAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCGGAAAGTAAAGTCGCACCCCAGTCCATCAGCGTGACTACCATCCCTg                      >  2:1418149/1‑150 (MQ=255)
                                                                                                                                aaTAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCGGAAAGTAAAGTCGCACCCCAGTCCATCAGCGTGACTACCATCCCTGc                     >  1:3822099/1‑150 (MQ=255)
                                                                                                                                aaTAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCGGAAAGTAAAGTCGCACCCCAGTCCATCAGCGTGACTACCATCCCTGc                     >  1:2509081/1‑150 (MQ=255)
                                                                                                                                aaTAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCGGAAAGTAAAGTCGCACCCCAGTCCATCAGCGTGACTACCATCCCTGc                     >  2:1649874/1‑150 (MQ=255)
                                                                                                                                   agagGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGGATACGGGCGGAAAGTAAAGTCGCACCCCAGTCCATCAGCGTGACTACCATCCCTGCgtt                  >  1:1105712/1‑150 (MQ=255)
                                                                                                                                    gagGCCCCCAGAAACTCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCGGAAAGTAAAGTCGCACCCCAGTCCATCAGCGTGACTACCATCCCTGCgttg                 <  1:3353590/150‑1 (MQ=255)
                                                                                                                                    gagGCCCCCAGAAACTCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCGGAAAGTAAAGTCGCACCCCAGTCCATCAGCGTGACTACCATCCCTGCgttg                 >  2:1816300/1‑150 (MQ=255)
                                                                                                                                     agGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCGGAAAGTAAAGTCGCACCCCAGTCCATCAGCGTGACTACCATCCCTGCgttgt                <  1:166570/150‑1 (MQ=255)
                                                                                                                                         ccacAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGACGAGCAGCGCCTAGCGTACGCAGCCATCGGAAAGCGGAATACGGGCGGAAAGTAAAGTCTCACCCCAGTCCATAAGAGTGAATACCATCCCTGCGTTGTTACg            >  1:1975554/4‑150 (MQ=255)
                                                                                                                                              gAAACGCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCGGAAAGTAAAGTCGCACCCCAGTCCATCAGCGTGACTACCATCCCTGCGTTGTTACGCAAAg       <  1:1050595/150‑1 (MQ=255)
                                                                                                                                               aaaCGCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCGGAAAGTAAAGTCGCACCCCAGTCCATCAGCGTGACTACCATCCCTGCGTTGTTACGCAAAGt      >  1:4491897/1‑150 (MQ=255)
                                                                                                                                                aaCGCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCGGAAAGTAAAGTCGCACCCCAGTCCATCAGCGTGACTACCATCCCTGCGTTGTTACGCAAAGtt     <  2:387044/150‑1 (MQ=255)
                                                                                                                                                   gcgGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCGGAAAGTAAAGTCGCACCCCAGTCCATCAGCGTGACTACCATCCCTGCGTTGTTACGCAAAGTTAAc  <  2:1242224/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TGCCAGCGACGTTTGTCGAACCCTTCCGGCCCGCCGTGCAGCTGGTTAACGCCCTGACTTGGCGAAAGCGTCACGGTTTCACCGTCAAAGGTATAACGGCTATTGGCGATACGGTTGGCATAACGACCAATAGAGGCCCCCAGAAACGCGGCCTGATCCTGATAGCATTCCGGGCTGGCACAGCCGAGCAGCGCCTCGCGGACGCTGCCATCGGAAAGCGGAATACGGGCGGAAAGTAAAGTCGCACCCCAGTCCATCAGCGTGACTACCATCCCTGCGTTGTTACGCAAAGTTAAC  >  CP009273/783947‑784243

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: