Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,143,375 A→G 8.7% intergenic (‑349/+249) ppc ← / ← argE phosphoenolpyruvate carboxylase/acetylornithine deacetylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,143,3750AG 8.7% 678.8 / 27.6 253intergenic (‑349/+249)ppc/argEphosphoenolpyruvate carboxylase/acetylornithine deacetylase
Reads supporting (aligned to +/- strand):  ref base A (70/153);  new base G (19/2);  total (93/160)
Fisher's exact test for biased strand distribution p-value = 1.28e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 255 total aligned reads displayed.

CGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAG  >  CP009273/4143227‑4143524
                                                                                                                                                    |                                                                                                                                                     
cGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAg                                                                                                                                                      <  1:3878363/150‑1 (MQ=255)
 gTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGAGCCGGATGCGGCGCGAGCGCCTTATCCGGCCAACAGt                                                                                                                                                     >  1:3814342/1‑150 (MQ=255)
  tAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGtt                                                                                                                                                    >  2:3292784/1‑150 (MQ=255)
      ttATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGc                                                                                                                                                <  1:3303463/150‑1 (MQ=255)
        aTGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGcct                                                                                                                                              >  1:3122599/1‑148 (MQ=255)
        aTGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCAt                                                                                                                                              <  1:4304515/150‑1 (MQ=255)
         tGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAGCAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTGGGGcggc                                                                                                                                             >  2:1270684/1‑147 (MQ=255)
         tGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATc                                                                                                                                             <  1:2101041/150‑1 (MQ=255)
         tGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATc                                                                                                                                             >  2:2375668/1‑150 (MQ=255)
          ggAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCg                                                                                                                                            <  1:358577/150‑1 (MQ=255)
            aaaTGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACAGTCTTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTCGGCCGGCCACCCGTAGGGCCTCgat                                                                                                                                          >  1:2282957/1‑148 (MQ=18)
              aTGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGtgtg                                                                                                                                        >  2:1383133/1‑150 (MQ=255)
              aTGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTg                                                                                                                                        <  1:2662502/150‑1 (MQ=255)
               tGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGa                                                                                                                                       <  1:2615486/150‑1 (MQ=255)
                gTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAg                                                                                                                                      >  1:905084/1‑150 (MQ=255)
                  tAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTc                                                                                                                                    <  1:997611/150‑1 (MQ=255)
                     aaaaTCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACt                                                                                                                                 <  1:4804753/150‑1 (MQ=255)
                     aaaaTCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACt                                                                                                                                 <  1:2873706/150‑1 (MQ=255)
                      aaaTCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTg                                                                                                                                >  1:2300190/1‑150 (MQ=255)
                        aTCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTgtc                                                                                                                              <  2:213718/150‑1 (MQ=255)
                         tCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCCGCGTTTGAGTCACTgtcg                                                                                                                             >  2:1543853/1‑150 (MQ=255)
                         tCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTgtcg                                                                                                                             <  2:3067581/150‑1 (MQ=255)
                         tCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTgtcg                                                                                                                             >  2:479059/1‑150 (MQ=255)
                          cGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTgtcgg                                                                                                                            <  2:4038650/150‑1 (MQ=255)
                           gCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTgtcggt                                                                                                                           >  2:3111724/1‑150 (MQ=255)
                                aaGTAACACCACAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGAt                                                                                                                      <  2:4269818/150‑1 (MQ=255)
                                aaGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGCGCGCCTTATCCGGCCTACAGGTGGGCATCGTTTGAGTCACTGTCGGTCGGAt                                                                                                                      >  1:577023/1‑150 (MQ=255)
                                aaGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGAt                                                                                                                      <  2:4269023/150‑1 (MQ=255)
                                aaGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGAt                                                                                                                      <  1:1740307/150‑1 (MQ=255)
                                    aaCACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGGTTGAGTCACTGTCGGTCGGATAAGa                                                                                                                  >  1:4445011/1‑150 (MQ=255)
                                     acacCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGAt                                                                                                                 <  2:382823/150‑1 (MQ=255)
                                      cacCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATg                                                                                                                <  1:4949039/150‑1 (MQ=255)
                                        ccAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATgcg                                                                                                              >  1:1808012/1‑150 (MQ=38)
                                          aaaGGTGTAGGTCGGATAAGATGCGCAAGTGTCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCa                                                                                                            >  1:4538755/1‑150 (MQ=255)
                                          aaaGGTGTAGGTCGGATAAGATGCGCAAGTATCGCCTCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCa                                                                                                            <  1:837809/150‑1 (MQ=255)
                                          aaaGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCa                                                                                                            >  1:4538954/1‑150 (MQ=40)
                                           aaGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCaa                                                                                                           >  2:2971642/1‑150 (MQ=40)
                                            aGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAg                                                                                                          >  1:333848/1‑150 (MQ=39)
                                             ggTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACCGTTGGGCAGCGGGTGAGTCACTGTCGGTCGGGGAAGATgcgcgcgt                                                                                                         >  2:4387852/1‑146 (MQ=11)
                                             ggTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGt                                                                                                         <  2:1209024/150‑1 (MQ=39)
                                                 tAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCg                                                                                                     <  1:4125919/150‑1 (MQ=38)
                                                   ggTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCa                                                                                                   <  1:500510/150‑1 (MQ=37)
                                                     tCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATc                                                                                                 <  2:2936344/150‑1 (MQ=37)
                                                      cGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATcc                                                                                                <  2:254309/150‑1 (MQ=37)
                                                       ggATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCg                                                                                               >  2:1145892/1‑150 (MQ=37)
                                                       ggATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCg                                                                                               >  2:4644731/1‑150 (MQ=37)
                                                        gagaaGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGa                                                                                              <  2:4914219/147‑1 (MQ=37)
                                                          tAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACa                                                                                            <  2:268702/150‑1 (MQ=37)
                                                           aaGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACat                                                                                           >  1:3505020/1‑150 (MQ=37)
                                                            aGATGCTCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCCACAGTGCCGTATGCGGCGCGTGCGCCTTATCCTGCCTACAGTTGGGAAGGGTTTGAGTCCCTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACatt                                                                                          <  2:3988452/150‑1 (MQ=11)
                                                             gATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACatta                                                                                         <  2:3928812/150‑1 (MQ=37)
                                                             gATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACatta                                                                                         <  2:3330748/150‑1 (MQ=37)
                                                               ggcgcAAGGATCGCATCCGACATTAGTGCGGCACTGGAGTTTGGCACCAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACattatt                                                                                       <  2:3636485/149‑1 (MQ=255)
                                                                gcgcAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTg                                                                                      <  2:554917/150‑1 (MQ=37)
                                                                gcgcAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTg                                                                                      <  1:5018387/150‑1 (MQ=37)
                                                                    aaGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGAATGGATTGAGTCAGTGTCGGTCGGATGAGTTGCGCAAGTATCGCATCCGACAGTATTGCGGc                                                                                  >  2:1084955/1‑150 (MQ=11)
                                                                     aGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCa                                                                                 >  1:1050658/1‑150 (MQ=37)
                                                                       tATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACt                                                                               >  2:1898945/1‑150 (MQ=37)
                                                                       tATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACt                                                                               <  1:2737833/150‑1 (MQ=37)
                                                                         tCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTgg                                                                             >  1:4762303/1‑150 (MQ=17)
                                                                          cGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGa                                                                            <  2:812161/150‑1 (MQ=37)
                                                                                cGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTgg                                                                      <  1:1964034/150‑1 (MQ=37)
                                                                                   cATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCaa                                                                   <  1:485585/150‑1 (MQ=37)
                                                                                   cATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTAGCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCaa                                                                   <  1:3067702/150‑1 (MQ=37)
                                                                                      tattGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCTAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAg                                                                >  1:3207436/1‑150 (MQ=37)
                                                                                        ttGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAgcg                                                              >  2:248168/1‑148 (MQ=11)
                                                                                           ctgCACTGTATATTGGCTACAGTTCCTGATGCGGCTCGAGCTCCATATCCTTCCTACAGTTGGGCATCTACTGAGTCACTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCg                                                           <  2:1511358/148‑1 (MQ=11)
                                                                                                cTGGAGTTTGGCTACCGTGCCGGTTGCGGCGCGAGCGCCTTATCCGTCATACAGTTGGGCATCGATTGACTAACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATgc                                                      <  2:2491376/150‑1 (MQ=14)
                                                                                                  ggAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATgcgg                                                    <  1:2664515/150‑1 (MQ=39)
                                                                                                     gTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGgcgc                                                 <  1:3842912/150‑1 (MQ=39)
                                                                                                      tttGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGcgcg                                                <  1:517216/150‑1 (MQ=39)
                                                                                                       ttGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGa                                               <  2:3239277/150‑1 (MQ=39)
                                                                                                       ttGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGa                                               >  2:4139804/1‑150 (MQ=39)
                                                                                                         ggCTACTGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGGTGGGCAGCGTGTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAgc                                             <  2:224796/150‑1 (MQ=14)
                                                                                                         ggCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAgc                                             >  1:1204114/1‑150 (MQ=39)
                                                                                                            aaCAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGcc                                          <  1:931638/150‑1 (MQ=39)
                                                                                                                gTGCCGGATGCGGCGCGAGCGCCTTATCAGTCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTTTCGCATCCGACATTATTGAGGCACTTGAGTTTGGCTACAGTGCCGGATGCGGCGCGAGCGCCTTAt                                      >  1:4736526/1‑150 (MQ=14)
                                                                                                                  gCCGTATGCGGCGCGTGCGCCTTATGCGGCCTTCAGTTGGGCATCTTTTGATTCCCTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATcc                                    <  2:503522/150‑1 (MQ=14)
                                                                                                                   ccGGATGCGGCTCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCg                                   <  1:2569143/150‑1 (MQ=255)
                                                                                                                     ggATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGc                                 >  1:4195148/1‑150 (MQ=39)
                                                                                                                      gATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTAt                                      <  2:908772/144‑1 (MQ=39)
                                                                                                                      gATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTAt                                      >  1:908772/1‑144 (MQ=39)
                                                                                                                             cgcgAGCGCCTTATCCGGCCTACCGTTGGGCATCGGTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGGATCGCATCCGACATTAGGGCGGCGCTGGAGGTTGGCAACCGTGCCGGCTGCGGCGCGAGCGCCTTATCCGGCCTACgggt                         >  2:1828283/1‑148 (MQ=11)
                                                                                                                              gcgAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTg                        <  2:3443211/150‑1 (MQ=255)
                                                                                                                               cgAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTgg                       <  1:1871221/150‑1 (MQ=255)
                                                                                                                               cgAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTgg                       <  1:1871241/150‑1 (MQ=255)
                                                                                                                               cgAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGAATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGAGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTgg                       >  2:1195039/1‑150 (MQ=255)
                                                                                                                                   cgcCTTATCCGGCCTACAGTTGGGCATCGATTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCAt                   <  1:1163444/150‑1 (MQ=255)
                                                                                                                                    gcCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCCGTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATc                  <  2:3982897/150‑1 (MQ=255)
                                                                                                                                    gcCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCCTc                  >  2:2063180/1‑150 (MQ=255)
                                                                                                                                     cgctaggcGGCCTACAGGCGGGCCGCGTGTGCGTCGCGGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCACCAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCg                 <  2:3749079/143‑1 (MQ=14)
                                                                                                                                       ttATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGAATCgga               >  2:4357878/1‑148 (MQ=255)
                                                                                                                                            cGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGt          <  1:1101268/150‑1 (MQ=255)
                                                                                                                                              gCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGGGTGAGTCa        >  2:11308/1‑150 (MQ=255)
                                                                                                                                              gCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCt                               <  2:2162570/127‑1 (MQ=25)
                                                                                                                                              gCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGAGCCGGATGCGGCGCGAGCGCCTTATCCGGCCt                               >  1:2162570/1‑127 (MQ=38)
                                                                                                                                               ccTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCAc       <  1:2971642/150‑1 (MQ=255)
                                                                                                                                                   cAGTTGGGCATCTTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTa   <  1:691730/150‑1 (MQ=255)
                                                                                                                                                   cAGTTGGGCATCGTTTGAGTCTCTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTa   <  1:2773428/150‑1 (MQ=255)
                                                                                                                                                    aGTTGGGCATCTTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAg  <  2:1707593/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAG  >  CP009273/4143227‑4143524

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: