Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,977,650 A→T 48.7% intergenic (+97/+62) aslB → / ← aslA putative AslA‑specific sulfatase‑maturating enzyme/putative Ser‑type periplasmic non‑aryl sulfatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,977,6500AT48.7% 28.6 / inf 320intergenic (+97/+62)aslB/aslAputative AslA‑specific sulfatase‑maturating enzyme/putative Ser‑type periplasmic non‑aryl sulfatase
Reads supporting (aligned to +/- strand):  ref base A (86/78);  new base T (64/92);  total (150/170)
Fisher's exact test for biased strand distribution p-value = 4.43e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 320 total aligned reads displayed.

CGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATATTTTTTCAGGATCTCCATATACGCGTGCATTTCGGTCTGTAGCGGTACACCC  >  CP009273/3977503‑3977796
                                                                                                                                                   |                                                                                                                                                  
cGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTTgg                                                                                                                                                  <  1:4330230/150‑1 (MQ=255)
cGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTTgg                                                                                                                                                  <  1:3582235/150‑1 (MQ=255)
  cccGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTTGGCt                                                                                                                                                <  1:2664611/150‑1 (MQ=255)
  cccGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCt                                                                                                                                                <  1:2550627/150‑1 (MQ=255)
     ggCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGa                                                                                                                                             <  2:607823/150‑1 (MQ=255)
             aCATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGc                                                                                                                                     >  2:4235683/1‑150 (MQ=255)
              cATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCg                                                                                                                                    <  1:1002334/150‑1 (MQ=255)
               aTTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTTGGCTGGATAAGATGCGt                                                                                                                                   >  1:4630114/1‑150 (MQ=255)
               aTTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTTGGCTGGATAAGATGCGt                                                                                                                                   <  1:3479837/150‑1 (MQ=255)
                ttATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTTGGCTGGATAAGATGCGTc                                                                                                                                  <  2:2760883/150‑1 (MQ=255)
                ttATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTc                                                                                                                                  <  2:1917477/150‑1 (MQ=255)
                 tATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCa                                                                                                                                 <  1:1196260/150‑1 (MQ=255)
                  atgcatgcGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATCAGATGCGTCAg                                                                                                                                >  2:4343732/1‑150 (MQ=255)
                  atgcatgcGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAg                                                                                                                                <  1:2839819/150‑1 (MQ=255)
                     catgcGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTTGGCTGGATAAGATGCGTCAgcat                                                                                                                             <  1:2340091/150‑1 (MQ=255)
                     catgcGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTTGGCTGGATAAGATGCGTCAgcat                                                                                                                             <  1:3053670/150‑1 (MQ=255)
                      atgcGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAgcatc                                                                                                                            <  1:3227025/150‑1 (MQ=255)
                      atgcGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAgcatc                                                                                                                            <  1:3226666/150‑1 (MQ=255)
                       tgcGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTTGGCTGGATAAGATGCGTCAgcatcg                                                                                                                           <  1:3261271/150‑1 (MQ=255)
                       tgcGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAgcatcg                                                                                                                           <  2:353364/150‑1 (MQ=255)
                       tgcGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAgcatcg                                                                                                                           >  1:1336766/1‑150 (MQ=255)
                         cgcATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAgcatcgca                                                                                                                         >  1:3842409/1‑150 (MQ=255)
                               tGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGc                                                                                                                   >  2:3471284/1‑150 (MQ=255)
                                 cTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCaa                                                                                                                 >  1:106858/1‑150 (MQ=255)
                                  tggtggTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCaaa                                                                                                                >  2:1138048/1‑150 (MQ=255)
                                             aagtaGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCa                                                                                                     <  1:1272792/150‑1 (MQ=255)
                                              agtaGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTTGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAg                                                                                                    >  2:2111459/1‑150 (MQ=255)
                                              agtaGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAg                                                                                                    <  1:1486723/150‑1 (MQ=255)
                                              agtaGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAAGTGCAGCGTTCTGTTGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAg                                                                                                    <  2:3963315/150‑1 (MQ=255)
                                                taGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTTGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCg                                                                                                  >  2:1488198/1‑150 (MQ=255)
                                                 aGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGa                                                                                                 >  2:3009655/1‑150 (MQ=255)
                                                    aaTCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTTGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAg                                                                                              <  1:3385227/150‑1 (MQ=255)
                                                    aaTCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAg                                                                                              <  1:405248/150‑1 (MQ=255)
                                                    aaTCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAg                                                                                              <  2:2771943/150‑1 (MQ=255)
                                                    aaTCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAg                                                                                              <  2:1985043/150‑1 (MQ=255)
                                                     aTCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTTGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAgc                                                                                             >  2:3185816/1‑150 (MQ=255)
                                                        ggcggcCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTTGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGcc                                                                                          <  2:303371/150‑1 (MQ=255)
                                                        ggcggcCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTTGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGcc                                                                                          >  2:2614743/1‑150 (MQ=255)
                                                        ggcggcCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGcc                                                                                          <  1:1720119/150‑1 (MQ=255)
                                                            gccgccTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTTGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCt                                                                                      <  1:2210756/150‑1 (MQ=255)
                                                              cgccTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTa                                                                                    <  2:432941/150‑1 (MQ=255)
                                                               gccTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTTGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAg                                                                                   >  2:3815930/1‑150 (MQ=255)
                                                               gccTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAg                                                                                   >  2:2945491/1‑150 (MQ=255)
                                                                ccTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGt                                                                                  <  1:1285701/150‑1 (MQ=255)
                                                                  tGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCa                                                                                >  2:2315649/1‑150 (MQ=255)
                                                                    cGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGa                                                                              >  1:717340/1‑150 (MQ=255)
                                                                     ggTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTTGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGAt                                                                             >  2:1847884/1‑150 (MQ=255)
                                                                      gTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTTGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGAtt                                                                            >  1:3159278/1‑150 (MQ=255)
                                                                       ttGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTTGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGAttt                                                                           <  2:2849292/150‑1 (MQ=255)
                                                                       ttGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTTGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGAttt                                                                           <  2:1446623/150‑1 (MQ=255)
                                                                       ttGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGAttt                                                                           <  2:4567130/150‑1 (MQ=255)
                                                                          aTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTTGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAAt                                                                        >  1:4488393/1‑150 (MQ=255)
                                                                           ttGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTTGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATc                                                                       <  2:450273/150‑1 (MQ=255)
                                                                            tGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTTGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCt                                                                      <  1:2287468/150‑1 (MQ=255)
                                                                               cGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTTGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTgcg                                                                   <  1:3941215/150‑1 (MQ=255)
                                                                               cGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTTGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTgcg                                                                   <  1:2961249/150‑1 (MQ=255)
                                                                                ggATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTgcgc                                                                  >  2:1694379/1‑150 (MQ=255)
                                                                                ggATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTgcgc                                                                  >  2:52500/1‑150 (MQ=255)
                                                                                 gATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTgcgcg                                                                 <  2:2914009/150‑1 (MQ=255)
                                                                                      ggcgTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTTGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCgtgg                                                            >  2:3383011/1‑150 (MQ=255)
                                                                                       gcgTAAACGCCTTCTCCGGCCTACATGATCGTGCAAATTCACTAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCgtggt                                                           <  1:2468431/150‑1 (MQ=255)
                                                                                        cgTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTTGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCgtggtg                                                          <  1:1422536/150‑1 (MQ=255)
                                                                                        cgTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCgtggtg                                                          >  1:1849825/1‑150 (MQ=255)
                                                                                          tAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGa                                                        <  1:3774857/150‑1 (MQ=255)
                                                                                             aCGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGatat                                                     <  2:2194784/150‑1 (MQ=255)
                                                                                               gCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTTGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATAttt                                                   <  2:2472801/150‑1 (MQ=255)
                                                                                               gCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATAttt                                                   >  1:1517763/1‑150 (MQ=255)
                                                                                                      ccGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATATTTTTTCAgg                                            >  2:901777/1‑150 (MQ=255)
                                                                                                      ccGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATATTTTTTCAgg                                            >  2:974462/1‑150 (MQ=255)
                                                                                                        ggCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTTGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATATTTTTTCAGGAt                                          >  1:1636001/1‑150 (MQ=255)
                                                                                                        ggCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTTGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATATTTTTTCAGGAt                                          >  1:4221019/1‑150 (MQ=255)
                                                                                                        ggCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTTGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATATTTTTTCAGGAt                                          <  2:3625574/150‑1 (MQ=255)
                                                                                                         gCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTTGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATATTTTTTCAGGAtc                                         >  2:3822291/1‑150 (MQ=255)
                                                                                                          ccTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATATTTTTTCAGGAtct                                        >  2:1368686/1‑150 (MQ=255)
                                                                                                             aCATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTTGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATATTTTTTCAGGATCTCCa                                     <  2:3706396/150‑1 (MQ=255)
                                                                                                             aCATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATATTTTTTCAGGATCTCCa                                     <  1:4483965/150‑1 (MQ=255)
                                                                                                               aTGATCGTGCAAATTCACTAAATTGCAGCGTTCTGTTGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATATTTTTTCAGGATCTCCata                                   <  1:851779/150‑1 (MQ=255)
                                                                                                               aTGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTTGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATATTTTTTCAGGATCTCCata                                   >  1:1887020/1‑150 (MQ=255)
                                                                                                                  aTCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATATTTTTTCAGGATCTCCATATAc                                <  1:1501144/150‑1 (MQ=255)
                                                                                                                     gTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATATTTTTTCAGGATCTCCATATAcgcg                             >  2:2638388/1‑150 (MQ=255)
                                                                                                                          aaTTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATATTTTTTCAGGATCTCCATATACGCGTGCAt                        >  2:164620/1‑150 (MQ=255)
                                                                                                                          aaTTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATATTTTTTCAGGATCTCCATATACGCGTGCAt                        >  1:3580632/1‑150 (MQ=255)
                                                                                                                            ttCAATAAATTGCAGCGTTCTGTTGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATATTTTTTCAGGATCTCCATATACGCGTGCAttt                      <  2:1166094/150‑1 (MQ=255)
                                                                                                                             tCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATATTTTTTCAGGATCTCCATATACGCGTGCATTTc                     <  1:1633066/150‑1 (MQ=255)
                                                                                                                             tCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATATTTTTTCAGGATCTCCATATACGCGTGCATTTc                     <  1:2767516/150‑1 (MQ=255)
                                                                                                                             tCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATATTTTTTCAGGATCTCCATATACGCGTGCATTTc                     >  2:1006056/1‑150 (MQ=255)
                                                                                                                               aaTAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATATTTTTTCAGGATCTCCATATACGCGTGCATTTCgg                   >  2:24057/1‑150 (MQ=255)
                                                                                                                                aTAAATTGCAGCGTTCTGTTGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATATTTTTTCAGGATCTCCATATACGCGTGCATTTCGGt                  <  1:612345/150‑1 (MQ=255)
                                                                                                                                aTAAATTGCAGCGTTCTGTTGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATATTTTTTCAGGATCTCCATATACGCGTGCATTTCGGt                  <  1:3521832/150‑1 (MQ=255)
                                                                                                                                aTAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATATTTTTTCAGGATCTCCATATACGCGTGCATTTCGGt                  <  2:4460339/150‑1 (MQ=255)
                                                                                                                                  aaaTTGCAGCGTTCTGTTGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATATTTTTTCAGGATCTCCATATACGCGTGCATTTCGGTCt                >  1:537980/1‑150 (MQ=255)
                                                                                                                                    aTTGCAGCGTTCTGTTGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATATTTTTTCAGGATCTCCATATACGCGTGCATTTCGGTCTGt              <  1:4293211/150‑1 (MQ=255)
                                                                                                                                      tGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATATTTTTTCAGGATCTCCATATACGCGTGCATTTCGGTCTGTAg            >  1:2421340/1‑150 (MQ=255)
                                                                                                                                          gCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATATTTTTTCAGGATCTCCATATACGCGTGCATTTCGGTCTGTAGCGGt        <  1:563373/150‑1 (MQ=255)
                                                                                                                                           cGTTCTGTTGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATATTTTTTCAGGATCTCCATATACGCGTGCATTTCGGTCTGTAGCGGTa       >  2:3259994/1‑150 (MQ=255)
                                                                                                                                            gTTCTGTTGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATATTTTTTCAGGATCTCCATATACGCGTGCATTTCGGTCTGTAGCGGTac      >  2:3236695/1‑150 (MQ=255)
                                                                                                                                               cTGTTGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATATTTTTTCAGGATCTCCATATACGCGTGCATTTCGGTCTGTAGCGGTACAcc   >  2:1511565/1‑150 (MQ=255)
                                                                                                                                               cTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATATTTTTTCAGGATCTCCATATACGCGTGCATTTCGGTCTGTAGCGGTACAcc   >  2:128605/1‑150 (MQ=255)
                                                                                                                                                tGTTGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATATTTTTTCAGGATCTCCATATACGCGTGCATTTCGGTCTGTAGCGGTACAccc  <  1:4553008/150‑1 (MQ=255)
                                                                                                                                                tGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATATTTTTTCAGGATCTCCATATACGCGTGCATTTCGGTCTGTAGCGGTACAccc  >  2:2157804/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATATTTTTTCAGGATCTCCATATACGCGTGCATTTCGGTCTGTAGCGGTACACCC  >  CP009273/3977503‑3977796

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: