Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,960,251 C→A 7.7% E515* (GAA→TAA)  ptsP ← fused PTS enzyme: PEP‑protein phosphotransferase (enzyme I)/GAF domain containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,960,2510CA 7.7% 1335.4 / 55.4 350E515* (GAA→TAA) ptsPfused PTS enzyme: PEP‑protein phosphotransferase (enzyme I)/GAF domain containing protein
Reads supporting (aligned to +/- strand):  ref base C (157/166);  new base A (15/12);  total (172/178)
Fisher's exact test for biased strand distribution p-value = 5.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

Only 100 of 358 total aligned reads displayed.

TTCATCGAGGCTTGTGACCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGA  >  CP009273/2960102‑2960400
                                                                                                                                                     |                                                                                                                                                     
ttCATCGAGGCTTGTGACCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTc                                                                                                                                                       >  1:1447237/1‑150 (MQ=255)
  cATCGAGGCTTGTGACCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGc                                                                                                                                                     >  2:3418807/1‑150 (MQ=255)
   aTCGAGGCTTGTGACCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCt                                                                                                                                                    <  1:1800879/150‑1 (MQ=255)
       aGGCTTGTGACCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATc                                                                                                                                                <  2:3037675/150‑1 (MQ=255)
        ggCTTGTGACCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTAGCTGATCg                                                                                                                                               <  1:852480/150‑1 (MQ=255)
           ttGTGACCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCa                                                                                                                                            >  1:1004723/1‑150 (MQ=255)
            tgtgACCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCAt                                                                                                                                           <  1:3750922/150‑1 (MQ=255)
            tgtgACCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCAt                                                                                                                                           <  2:4019423/150‑1 (MQ=255)
               gACCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTa                                                                                                                                        >  2:1339462/1‑150 (MQ=255)
                 ccATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAg                                                                                                                                      >  1:1105154/1‑150 (MQ=255)
                  cATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTAGCTGATCGGCATGTAAgg                                                                                                                                     <  2:2817979/150‑1 (MQ=255)
                   aTCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGAGCGGCATGTAAGGc                                                                                                                                    >  2:3873507/1‑150 (MQ=255)
                      ggCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAgc                                                                                                                                 >  1:1602009/1‑150 (MQ=255)
                      ggCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAgc                                                                                                                                 <  1:4699188/150‑1 (MQ=255)
                      ggCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTAGCTGATCGGCATGTAAGGCAgc                                                                                                                                 >  2:3557928/1‑150 (MQ=255)
                         aaCAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAgctgc                                                                                                                              <  1:2344863/150‑1 (MQ=255)
                          aCAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAgctgct                                                                                                                             <  2:1150990/150‑1 (MQ=255)
                          aCAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAgctgct                                                                                                                             <  2:2198926/150‑1 (MQ=255)
                            aGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTa                                                                                                                           >  2:4635653/1‑150 (MQ=255)
                             gAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTAt                                                                                                                          >  1:3224268/1‑150 (MQ=255)
                             gAAGATTCAGGTTGCCCGTAGCGGCATTAGCACTCAGCATCGCGCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTATTGCGTATCCCACGCCAACCCATGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTAt                                                                                                                          <  1:4773672/150‑1 (MQ=255)
                              aaTATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATc                                                                                                                         <  1:2158043/150‑1 (MQ=255)
                              aaTATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATc                                                                                                                         <  2:4808863/150‑1 (MQ=255)
                                  ttCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCt                                                                                                                     <  2:499277/150‑1 (MQ=255)
                                     aGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCg                                                                                                                  <  2:1895541/150‑1 (MQ=255)
                                        gtGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGAGTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACa                                                                                                               >  1:3587777/2‑150 (MQ=255)
                                            ccGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCa                                                                                                           <  1:395003/150‑1 (MQ=255)
                                             cGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAg                                                                                                          >  1:3766689/1‑150 (MQ=255)
                                             cGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAg                                                                                                          >  2:3489025/1‑150 (MQ=255)
                                               tAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCg                                                                                                        >  1:4435294/1‑150 (MQ=255)
                                                  cGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTAc                                                                                                     <  2:2663473/150‑1 (MQ=255)
                                                   ggCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACg                                                                                                    <  2:832115/150‑1 (MQ=255)
                                                   ggCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACg                                                                                                    >  1:1145325/1‑150 (MQ=255)
                                                   ggCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACg                                                                                                    >  1:1992904/1‑150 (MQ=255)
                                                        tAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAgg                                                                                               >  1:937447/1‑150 (MQ=255)
                                                         aGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGt                                                                                              >  1:631098/1‑150 (MQ=255)
                                                           cACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGa                                                                                            >  1:286503/1‑150 (MQ=255)
                                                            aCGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGAc                                                                                           >  1:4113634/1‑150 (MQ=255)
                                                              gcagcaTCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACgg                                                                                         <  1:4533532/150‑1 (MQ=255)
                                                              gcagcaTCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACgg                                                                                         >  2:2014131/1‑150 (MQ=255)
                                                              gcagcaTCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGAGCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACgg                                                                                         >  1:4626631/1‑150 (MQ=255)
                                                               cagcaTCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACggg                                                                                        <  1:158412/150‑1 (MQ=255)
                                                                agcaTCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGt                                                                                       >  1:1224517/1‑150 (MQ=255)
                                                                 gcaTCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGtt                                                                                      <  1:781602/150‑1 (MQ=255)
                                                                 gcaTCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGtt                                                                                      <  1:304725/150‑1 (MQ=255)
                                                                  caTCGCCCGCACCTGGATCAAGACGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGttt                                                                                     <  2:2519641/150‑1 (MQ=255)
                                                                      gCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCa                                                                                 <  1:333121/150‑1 (MQ=255)
                                                                       cccGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCAt                                                                                <  2:4097788/150‑1 (MQ=255)
                                                                         cGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTg                                                                              <  1:3108814/150‑1 (MQ=255)
                                                                                ggATCTAGATGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCTACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACAttt                                                                       <  2:376125/150‑1 (MQ=255)
                                                                                ggATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACAttt                                                                       <  1:588259/150‑1 (MQ=255)
                                                                                 gATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTg                                                                      >  2:930479/1‑150 (MQ=255)
                                                                                  aTCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTgc                                                                     <  1:1749750/150‑1 (MQ=255)
                                                                                   tCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTgca                                                                    <  1:1349496/150‑1 (MQ=255)
                                                                                     aagaagATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTgcagc                                                                  <  1:1173922/150‑1 (MQ=255)
                                                                                     aagaagATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTgcagc                                                                  <  2:4384662/150‑1 (MQ=255)
                                                                                     aagaagATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTgcagc                                                                  <  1:2486942/150‑1 (MQ=255)
                                                                                      agaagaTCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTgcagca                                                                 >  2:3405679/1‑150 (MQ=255)
                                                                                        aagaTCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATc                                                               <  1:686849/150‑1 (MQ=255)
                                                                                          gaTCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATccc                                                             >  1:3896502/1‑150 (MQ=255)
                                                                                           aTCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATcccc                                                            >  1:4855505/1‑150 (MQ=255)
                                                                                           aTCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTAGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATcccc                                                            >  2:2590375/1‑150 (MQ=255)
                                                                                            tctcCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCCGCGTACGCACGGTGACGGGGTTATCATTGACCATTTGCAGCATCCCCt                                                           >  2:3150527/1‑150 (MQ=255)
                                                                                            tctcCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTAGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCt                                                           <  2:726573/150‑1 (MQ=255)
                                                                                                cGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATa                                                       >  2:3861225/1‑150 (MQ=255)
                                                                                                cGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATa                                                       >  2:786217/1‑150 (MQ=255)
                                                                                                 ggCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATAc                                                      <  1:2952495/150‑1 (MQ=255)
                                                                                                   cTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTg                                                    >  2:2901229/1‑150 (MQ=255)
                                                                                                   cTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGAATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTg                                                    <  1:3962239/150‑1 (MQ=255)
                                                                                                   cTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTAGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTg                                                    >  1:4265602/1‑150 (MQ=255)
                                                                                                    tGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTgc                                                   <  2:2278212/150‑1 (MQ=255)
                                                                                                    tGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTgc                                                   <  2:156048/150‑1 (MQ=255)
                                                                                                      aTCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTgcgc                                                 >  1:2707511/1‑150 (MQ=255)
                                                                                                                 tGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGttctt                                      <  1:1020268/150‑1 (MQ=255)
                                                                                                                  gCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGttcttc                                     >  2:1259205/1‑150 (MQ=255)
                                                                                                                    gAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGttcttctt                                   <  2:2682216/150‑1 (MQ=255)
                                                                                                                    gAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGttcttctt                                   <  1:2746346/150‑1 (MQ=255)
                                                                                                                         ccACGCCAACCCAGGCATGGATTCTCTTAGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGAc                              <  1:4193126/150‑1 (MQ=255)
                                                                                                                          cACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACg                             <  2:644617/150‑1 (MQ=255)
                                                                                                                          cACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACg                             >  2:4940083/1‑150 (MQ=255)
                                                                                                                           aCGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACgg                            >  1:3479773/1‑150 (MQ=255)
                                                                                                                            cGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACggg                           >  2:1710863/1‑150 (MQ=255)
                                                                                                                             gCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGa                          <  2:1531744/150‑1 (MQ=255)
                                                                                                                               cAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGaaa                        >  2:1991493/1‑150 (MQ=255)
                                                                                                                                aaCCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAAc                       >  2:3846917/1‑150 (MQ=255)
                                                                                                                                aaCCCAGGCATGGATTCTCTTAGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAAc                       >  2:2859506/1‑150 (MQ=255)
                                                                                                                                 aCCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAAcc                      <  2:1447237/150‑1 (MQ=255)
                                                                                                                                 aCCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAAcc                      <  1:54936/150‑1 (MQ=255)
                                                                                                                                 aCCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAAcc                      >  2:1813046/1‑150 (MQ=255)
                                                                                                                                 aCCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCCGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAAcc                      <  2:4702203/150‑1 (MQ=255)
                                                                                                                                       gCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTg                >  1:5024871/1‑150 (MQ=255)
                                                                                                                                        cATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTgc               >  1:1511096/1‑150 (MQ=255)
                                                                                                                                         aTGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTgca              >  2:3014209/1‑150 (MQ=255)
                                                                                                                                           ggATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTgcagc            >  2:4257070/1‑150 (MQ=255)
                                                                                                                                            gATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTgcagca           <  2:2145645/150‑1 (MQ=255)
                                                                                                                                             aTTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCAGCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCAt          <  2:2446120/150‑1 (MQ=255)
                                                                                                                                                 tctTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAAt      <  2:3373501/150‑1 (MQ=255)
                                                                                                                                                 tctTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAAt      <  2:3601902/150‑1 (MQ=255)
                                                                                                                                                 tctTAGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAAt      >  2:1323465/1‑150 (MQ=255)
                                                                                                                                                     cGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGa  <  1:912643/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTCATCGAGGCTTGTGACCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGA  >  CP009273/2960102‑2960400

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: