Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,224,419 C→A 7.4% R725S (CGT→AGT)  ygjK → alpha‑glucosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,224,4190CA 7.4% 965.6 / 33.8 261R725S (CGT→AGT) ygjKalpha‑glucosidase
Reads supporting (aligned to +/- strand):  ref base C (118/123);  new base A (10/9);  total (129/132)
Fisher's exact test for biased strand distribution p-value = 8.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.88e-01

Only 100 of 261 total aligned reads displayed.

GCTCGATCCTAAAGAGTTCAACACCTTTGTCCCGCTGGGAACGGCGGCGTTAACCAATCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCATTTGTATA  >  CP009273/3224271‑3224566
                                                                                                                                                    |                                                                                                                                                   
 cTCGATCCTAAAGAGTTCAACACCTTTGTCCCGCTGGGAACGGCGGCGTTAACCAATCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGt                                                                                                                                                   >  1:589614/1‑150 (MQ=255)
 cTCGATCCTAAAGAGTTCAACACCTTTGTCCCGCTGGGAACGGCGGCGTTAACCAATCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGt                                                                                                                                                   >  2:3055429/1‑150 (MQ=255)
    gATCCTAAAGAGTTCAACACCTTTGTCCCGCTGGGAACGGCGGCGTTAACCAATCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGtgat                                                                                                                                                >  1:2924672/1‑150 (MQ=255)
    gATCCTAAAGAGTTCAACACCTTTGTCCCGCTGGGAACGGCGGCGTTAACCAATCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGtgat                                                                                                                                                >  2:2346213/1‑150 (MQ=255)
    gATCCTAAAGAGTTCAACACCTTTGTCCCGCTGGGAACGGCGGCGTTAACCAATCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGtgat                                                                                                                                                >  2:3684381/1‑150 (MQ=255)
      tCCTAAAGAGTTCAACACCTTTGTCCCGCTGGGAACGGCGGCGTTAACCAATCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGtgatga                                                                                                                                              <  2:841504/150‑1 (MQ=255)
      tCCTAAAGAGTTCAACACCTTTGTCCCGCTGGGAACGGCGGCGTTAACCAATCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGtgatga                                                                                                                                              >  2:2935039/1‑150 (MQ=255)
       ccTAAAGAGTTCAACACCTTTGTCCCGCTGGGAACGGCGGCGTTAACCAATCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTgatgat                                                                                                                                             >  2:3572595/1‑150 (MQ=255)
         tAAAGAGTTCAACACCTTTGTCCCGCTGGGAACGGCGGCGTTAACCAATCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGc                                                                                                                                           >  1:3394875/1‑150 (MQ=255)
          aaaGAGTTCAACACCTTTGTCCCGCTGGGAACGGCGGCGTTAACCAATCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGcc                                                                                                                                          <  1:2167097/150‑1 (MQ=255)
          aaaGAGTTCAACACCTTTGTCCCGCTGGGAACGGCGGCGTTAACCAATCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGcc                                                                                                                                          <  1:1103872/150‑1 (MQ=255)
           aaGAGTTCAACACCTTTGTCCCGCTGGGAACGGCGGCGTTAACCAATCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATGCAGATGATGccc                                                                                                                                         >  2:3677440/1‑150 (MQ=255)
            agagTTCAACACCTTTGTCCCGCTGGGAACGGCGGCGTTAACCAATCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCt                                                                                                                                        <  1:187430/150‑1 (MQ=255)
              agTTCAACACCTTTGTCCCGCTGGGAACGGCGGCGTTAACCAATCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGa                                                                                                                                      >  2:732836/1‑150 (MQ=255)
                ttCAACACCTTTGTCCCGCTGGGAACGGCGGCGTTAACCAATCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAg                                                                                                                                    >  2:4283973/1‑150 (MQ=255)
                   aaCACCTTTGTCCCGCTGGGAACGGCGGCGTTAACCAATCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTg                                                                                                                                 >  2:43801/1‑150 (MQ=255)
                   aaCACCTTTGTCCCGCTGGGAACGGCGGCGTTAACCAATCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTg                                                                                                                                 >  1:322817/1‑150 (MQ=255)
                        cTTTGTCCCGCTGGGAACGGCGGCGTTAACCAATCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGa                                                                                                                            <  1:3596512/150‑1 (MQ=255)
                          ttGTCCCGCTGGGAACGGCGGCGTTAACCAATCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATa                                                                                                                          >  1:3131037/1‑150 (MQ=255)
                           tGTCCCGCTGGGAACGGCGGCGTTAACCAATCCGGCTTTTTGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGTAGCTTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATAc                                                                                                                         >  2:3263460/1‑150 (MQ=255)
                              cccGCTGGGAACGGCGGCGTTAACCAATCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGtt                                                                                                                      <  2:4318844/150‑1 (MQ=255)
                              cccGCTGGGAACGGCGGCGTTAACCAATCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGtt                                                                                                                      <  2:2073046/150‑1 (MQ=255)
                                cGCTGGGAACGGCGGCGTTAACCAATCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGttct                                                                                                                    >  2:1068993/1‑150 (MQ=255)
                                cGCTGGGAACGGCGGCGTTAACCAATCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGttct                                                                                                                    <  2:1268197/150‑1 (MQ=255)
                                 gCTGGGAACGGCGGCGTTAACCAATCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGttctt                                                                                                                   >  1:2907399/1‑150 (MQ=255)
                                  cTGGGAACGGCGGCGTTAACCAATCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGttcttc                                                                                                                  >  1:4067442/1‑150 (MQ=255)
                                   tGGGAACGGCGGCGTTAACCAATCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATAGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTcc                                                                                                                 <  1:2987234/150‑1 (MQ=255)
                                   tGGGAACGGCGGCGTTAACCAATCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATAGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTcc                                                                                                                 <  1:708671/150‑1 (MQ=255)
                                   tGGGAACGGCGGCGTTAACCAATCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATAGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTcc                                                                                                                 <  1:3436749/150‑1 (MQ=255)
                                       aaCGGCGGCGTTAACCAATCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCa                                                                                                             >  2:2554118/1‑150 (MQ=255)
                                               cgTTAACCAATCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAg                                                                                                     >  1:1592808/1‑150 (MQ=255)
                                                gTTAACCAATCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAgg                                                                                                    >  2:79280/1‑150 (MQ=255)
                                                   aaCCAATCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGtt                                                                                                 <  1:4421144/150‑1 (MQ=255)
                                                       aaTCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAAcc                                                                                             >  1:1930049/1‑150 (MQ=255)
                                                       aaTCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAAcc                                                                                             >  1:3642975/1‑150 (MQ=255)
                                                           cGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAAccgccg                                                                                         >  1:453266/1‑150 (MQ=255)
                                                              cTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATg                                                                                      >  2:668241/1‑150 (MQ=255)
                                                                   ggCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCa                                                                                 <  1:520811/150‑1 (MQ=255)
                                                                      gcTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAAtt                                                                              >  2:1258882/1‑150 (MQ=255)
                                                                      gcTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAAtt                                                                              <  2:133045/150‑1 (MQ=255)
                                                                       cTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTc                                                                             <  1:4165799/150‑1 (MQ=255)
                                                                        tGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCa                                                                            >  2:3717457/1‑150 (MQ=255)
                                                                         gATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAg                                                                           <  2:895544/150‑1 (MQ=255)
                                                                          atatCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAgg                                                                          <  1:1646375/150‑1 (MQ=255)
                                                                          atatCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAgg                                                                          <  2:493270/150‑1 (MQ=255)
                                                                           tatCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGa                                                                         >  1:1913735/1‑150 (MQ=255)
                                                                               tACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAAt                                                                     <  1:186508/150‑1 (MQ=255)
                                                                                aCTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAAtt                                                                    <  1:3677980/150‑1 (MQ=255)
                                                                                  tGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTAc                                                                  >  1:2052490/1‑150 (MQ=255)
                                                                                   ggCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACa                                                                 >  2:4426964/1‑150 (MQ=255)
                                                                                   ggCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACa                                                                 <  2:4058232/150‑1 (MQ=255)
                                                                                   ggCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACa                                                                 <  2:1886713/150‑1 (MQ=255)
                                                                                         ggCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGc                                                           >  2:1861291/1‑150 (MQ=255)
                                                                                          gcgcGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCt                                                          <  2:1206155/150‑1 (MQ=255)
                                                                                           cgcgTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTg                                                         <  2:4395434/150‑1 (MQ=255)
                                                                                           cgcgTATGGGTGGATCAGTTCTGGTTTGGTCGGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTg                                                         <  2:2152537/150‑1 (MQ=255)
                                                                                              gTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACa                                                      >  2:688939/1‑150 (MQ=255)
                                                                                               tATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATAGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAg                                                     >  2:1534365/1‑150 (MQ=255)
                                                                                                aTGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAgg                                                    >  1:3946262/1‑150 (MQ=255)
                                                                                                aTGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAgg                                                    >  2:2139844/1‑150 (MQ=255)
                                                                                                 tGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGgc                                                   >  1:1594502/1‑150 (MQ=255)
                                                                                                 tGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGgc                                                   >  1:1955155/1‑150 (MQ=255)
                                                                                                   ggTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGgcgc                                                 <  1:1011903/150‑1 (MQ=255)
                                                                                                    gTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGca                                                >  1:772633/1‑150 (MQ=255)
                                                                                                       gATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAg                                             <  1:1069240/150‑1 (MQ=255)
                                                                                                        aTCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGc                                            <  1:1718285/150‑1 (MQ=255)
                                                                                                        aTCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGc                                            <  1:1262549/150‑1 (MQ=255)
                                                                                                           aGTACTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAg                                         >  1:1673704/1‑150 (MQ=255)
                                                                                                            gTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAgg                                        >  2:2812324/1‑150 (MQ=255)
                                                                                                              tCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGgcg                                      <  2:1537522/150‑1 (MQ=255)
                                                                                                              tCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGgcg                                      <  2:3798697/150‑1 (MQ=255)
                                                                                                                tGGTTTGGTCTGATAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCa                                    <  2:2653406/150‑1 (MQ=255)
                                                                                                                tGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCa                                    <  2:2653654/150‑1 (MQ=255)
                                                                                                                tGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCa                                    <  2:3206839/150‑1 (MQ=255)
                                                                                                                 ggTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCAc                                   >  1:2783611/1‑150 (MQ=255)
                                                                                                                    ttGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCaa                                >  1:3081537/1‑150 (MQ=255)
                                                                                                                         cTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTtc                           <  1:1187218/150‑1 (MQ=255)
                                                                                                                          tGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTtct                          >  1:2090465/1‑150 (MQ=255)
                                                                                                                          tGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTtct                          <  1:419414/150‑1 (MQ=255)
                                                                                                                          tGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTtct                          <  1:371328/150‑1 (MQ=255)
                                                                                                                          tGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTtct                          >  1:1567969/1‑150 (MQ=255)
                                                                                                                          tGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTtct                          >  1:1569062/1‑150 (MQ=255)
                                                                                                                          tGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTAAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAAtttat                          >  1:3323222/1‑148 (MQ=255)
                                                                                                                          tGAAAGGGATGAAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTtct                          <  2:3601656/150‑1 (MQ=255)
                                                                                                                            aaaGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTcc                        <  2:3532379/150‑1 (MQ=255)
                                                                                                                             aaGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCt                       <  1:3572595/150‑1 (MQ=255)
                                                                                                                             aaGGGATGGAGCGTTACGGTTATAGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCt                       >  2:1648820/1‑150 (MQ=255)
                                                                                                                             aaGGGATGGAGCGTTACGGTTATAGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCt                       >  2:526776/1‑150 (MQ=255)
                                                                                                                              aGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTg                      <  1:2430774/150‑1 (MQ=255)
                                                                                                                                  aTGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGt                  <  1:433128/150‑1 (MQ=255)
                                                                                                                                  aTGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGt                  <  1:427472/150‑1 (MQ=255)
                                                                                                                                     gAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGcc               <  2:2239965/150‑1 (MQ=255)
                                                                                                                                       gCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGGTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATGTCTCCTGGAGTGCcgc             >  2:3788272/1‑150 (MQ=255)
                                                                                                                                             cGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCAttt       <  1:688939/150‑1 (MQ=255)
                                                                                                                                              ggTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCATTTg      <  1:1353203/150‑1 (MQ=255)
                                                                                                                                                ttATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCATTTGta    >  1:3527953/1‑150 (MQ=255)
                                                                                                                                                ttATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCATTTGta    >  1:1255936/1‑150 (MQ=255)
                                                                                                                                                 tATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCATTTGtat   >  1:3904610/1‑150 (MQ=255)
                                                                                                                                                 tATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCATTTGtat   <  1:2285076/150‑1 (MQ=255)
                                                                                                                                                  aTCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCATTTGtata  >  1:792704/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GCTCGATCCTAAAGAGTTCAACACCTTTGTCCCGCTGGGAACGGCGGCGTTAACCAATCCGGCTTTTGGCGCTGATATCTACTGGCGTGGGCGCGTATGGGTGGATCAGTTCTGGTTTGGTCTGAAAGGGATGGAGCGTTACGGTTATCGTGATGATGCCCTGAAGCTGGCGGATACGTTCTTCCGGCACGCCAAAGGGTTAACCGCCGATGGACCAATTCAGGAAAATTACAACCCGCTGACAGGCGCACAGCAAGGCGCACCAAATTTCTCCTGGAGTGCCGCGCATTTGTATA  >  CP009273/3224271‑3224566

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: