Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,889,399 A→T 21.8% D113V (GAT→GTT)  pabB → aminodeoxychorismate synthase, subunit I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,889,3990AT21.8% 756.1 / 213.6 316D113V (GAT→GTT) pabBaminodeoxychorismate synthase, subunit I
Reads supporting (aligned to +/- strand):  ref base A (120/127);  new base T (35/34);  total (155/161)
Fisher's exact test for biased strand distribution p-value = 7.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 319 total aligned reads displayed.

ACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCGGATATGGCAGTGGGTATCTACGATTGGGCGCTCATTGTCGACCACCAGCGTCATACAGTTTCTTTGCTGAGTCATAATGATGTCAATGCCC  >  CP009273/1889251‑1889548
                                                                                                                                                    |                                                                                                                                                     
 ccGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATc                                                                                                                                                     >  2:3845475/1‑150 (MQ=255)
  cGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGTTCt                                                                                                                                                    <  2:3373092/150‑1 (MQ=255)
   gttgttAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTg                                                                                                                                                   <  1:2235555/150‑1 (MQ=255)
    ttgttAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGTTCTgg                                                                                                                                                  <  2:1114628/150‑1 (MQ=255)
    ttgttAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTgg                                                                                                                                                  <  1:1055070/150‑1 (MQ=255)
          gTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGTTCTGGgccgcc                                                                                                                                            <  2:377269/150‑1 (MQ=255)
          gTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGgccgcc                                                                                                                                            >  1:1831248/1‑150 (MQ=255)
           tGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGccgccg                                                                                                                                           <  1:414822/150‑1 (MQ=255)
             aaaGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGtt                                                                                                                                         >  1:269197/1‑150 (MQ=255)
               aGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGtttt                                                                                                                                       <  2:2678153/150‑1 (MQ=255)
               aGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGtttt                                                                                                                                       <  2:2025990/150‑1 (MQ=255)
                 cGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGTTCTGGGCCGCCGTTTTGa                                                                                                                                     <  1:3617621/150‑1 (MQ=255)
                 cGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGTTCTGGGCCGCCGTTTTGa                                                                                                                                     <  1:2651410/150‑1 (MQ=255)
                  gAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAg                                                                                                                                    <  2:461915/150‑1 (MQ=255)
                  gAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAg                                                                                                                                    <  1:1323768/150‑1 (MQ=255)
                  gAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAg                                                                                                                                    <  2:3037723/150‑1 (MQ=255)
                  gAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAg                                                                                                                                    >  1:2610083/1‑150 (MQ=255)
                   aaaaaCGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGt                                                                                                                                   >  2:4099287/1‑150 (MQ=255)
                   aaaaaCGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGt                                                                                                                                   >  1:4484856/1‑150 (MQ=255)
                   aaaaaCGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGt                                                                                                                                   >  2:857331/1‑150 (MQ=255)
                    aaaaCGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTc                                                                                                                                  >  2:4179854/1‑150 (MQ=255)
                     aaaCGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCa                                                                                                                                 <  2:1621732/150‑1 (MQ=255)
                      aaCGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCAc                                                                                                                                <  2:2818021/150‑1 (MQ=255)
                       aCGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACt                                                                                                                               >  2:2221032/1‑150 (MQ=255)
                        cGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCTCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGTTCTGGGCCGCCGTTTTGAGTCACTg                                                                                                                              >  2:2496579/1‑150 (MQ=255)
                        cGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGTTCTGGGCCGCCGTTTTGAGTCACTg                                                                                                                              >  2:1409255/1‑150 (MQ=255)
                         gCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGc                                                                                                                             >  1:1168535/1‑150 (MQ=255)
                          cacaACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTGTTGAGTCACTGcc                                                                                                                            >  1:4672710/1‑150 (MQ=255)
                           acaACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCa                                                                                                                           <  1:332599/150‑1 (MQ=255)
                            caACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGTTCTGGGCCGCCGTTTTGAGTCACTGCCAg                                                                                                                          <  1:1232798/150‑1 (MQ=255)
                            caACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAg                                                                                                                          >  2:3325609/1‑150 (MQ=255)
                             aaCGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGa                                                                                                                         >  1:754868/1‑150 (MQ=255)
                             aaCGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGa                                                                                                                         >  2:3853222/1‑150 (MQ=255)
                               cGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGaaa                                                                                                                       >  1:4739004/1‑150 (MQ=255)
                                gACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAAt                                                                                                                      <  1:4579359/150‑1 (MQ=255)
                                    aCTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCg                                                                                                                  >  2:440434/1‑150 (MQ=255)
                                        atgaCCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAAc                                                                                                              <  1:2916309/150‑1 (MQ=255)
                                           aCCCGCTACAGGTGCCCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGAGCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAg                                                                                                           >  2:1764663/1‑150 (MQ=255)
                                            cccGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGa                                                                                                          >  2:4517089/1‑150 (MQ=255)
                                              cGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGata                                                                                                        <  1:1199949/150‑1 (MQ=255)
                                              cGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGata                                                                                                        >  1:4527213/1‑150 (MQ=255)
                                              cGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGata                                                                                                        <  1:283197/150‑1 (MQ=255)
                                              cGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGata                                                                                                        >  1:3272356/1‑150 (MQ=255)
                                                cTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTCCGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATc                                                                                                      >  2:2435978/1‑150 (MQ=255)
                                                cTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATc                                                                                                      <  1:2427469/150‑1 (MQ=255)
                                                          tCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCgg                                                                                            <  1:1460205/150‑1 (MQ=255)
                                                          tCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCgg                                                                                            <  2:3083039/150‑1 (MQ=255)
                                                               cagGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCGGATATg                                                                                       <  1:3974368/150‑1 (MQ=255)
                                                               cagGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCGGATATg                                                                                       >  1:1622780/1‑150 (MQ=255)
                                                                agGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCGGATATgg                                                                                      <  2:34105/150‑1 (MQ=255)
                                                                 ggTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCGGATATGGc                                                                                     <  2:1524576/150‑1 (MQ=255)
                                                                  gTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGTTCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCGGATATGGCa                                                                                    <  2:654552/150‑1 (MQ=255)
                                                                    gCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCGGATATGGCAGt                                                                                  >  1:2367010/1‑150 (MQ=255)
                                                                      tGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGTTCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCGGATATGGCAGTgg                                                                                >  1:2979430/1‑150 (MQ=255)
                                                                       ggATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCGGATATGGCAGTggg                                                                               >  2:4045069/1‑150 (MQ=255)
                                                                        gATCGCGCAGACATTCGCCCAACGCATAACGAAGATGTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCGGATATGGCAGTGGGt                                                                              <  1:1904599/150‑1 (MQ=255)
                                                                            gcgcAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCGGATATGGCAGTGGGTATCt                                                                          <  1:2254103/150‑1 (MQ=255)
                                                                              gcAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGTTCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCGGATATGGCAGTGGGTATCTAc                                                                        <  2:473744/150‑1 (MQ=255)
                                                                               cAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGTTCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCGGATATGGCAGTGGGTATCTACg                                                                       <  1:418646/150‑1 (MQ=255)
                                                                                    aTTAGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGTTCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGATCTGCCGGATATGGCAGTGGGTCTCTCAGATTgg                                                                  >  2:4239633/1‑150 (MQ=255)
                                                                                            aaCGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCGGATATGGCAGTGGGTATCTACGATTGGGCGCTCAt                                                          <  1:103386/150‑1 (MQ=255)
                                                                                               gCATAACGAAGATTTGCCATTTCAGGGTGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCGGATATGGCAGTGGGTATCTCCGATTGGGCGCTCATTGt                                                       >  2:3267022/1‑150 (MQ=255)
                                                                                                cATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGTTCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCGGATATGGCAGTGGGTATCTACGATTGGGCGCTCATTGTc                                                      <  1:1411619/150‑1 (MQ=255)
                                                                                                cATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCGGATATGGCCGTGGGTATCTACGATTGGGCGCTCATTGTc                                                      >  1:960136/1‑150 (MQ=255)
                                                                                                cATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCGGATATGGCAGTGGGTATCTACGATTGGGCGCTCATTGTc                                                      >  1:3159185/1‑150 (MQ=255)
                                                                                                    aCGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCGGATATGGCAGTGGGTATCTACGATTGGGCGCTCATTGTCGacc                                                  >  1:1393421/1‑150 (MQ=255)
                                                                                                    aCGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCGGATATGGCAGTGGGTATCTACGATTGGGCGCTCATTGTCGacc                                                  >  1:2789886/1‑150 (MQ=255)
                                                                                                     cGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCGGATATGGCAGTGGGTATCTACGATTGGGCGCTCATTTTCGacca                                                 >  1:483573/1‑150 (MQ=255)
                                                                                                      gAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGTTCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCGGATATGGCAGTGGGTATCTACGATTGGGCGCTCATTGTCGaccac                                                >  2:376109/1‑150 (MQ=255)
                                                                                                        aGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCGGATATGGCAGTGGGTATCTACGATTGGGCGCTCATTGTCGAccacca                                              <  1:1921602/150‑1 (MQ=255)
                                                                                                          aTTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCGGATATGGCAGTGGGTATCTACGATTGGGCGCTCATTGTCGACCACCAGc                                            >  2:2373424/1‑150 (MQ=255)
                                                                                                          aTTTGCCAGTTCAGGGCGGCGCACTGGGGTGGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCGGATATGGCAGGGGGTATCTACGATTGGGCGCTCATTGTCGACCACCAGc                                            >  2:18457/1‑150 (MQ=255)
                                                                                                            ttGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGTTCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCGGATATGGCAGTGGGTATCTACGATTGGGCGCTCATTGTCGACCACCAGCGt                                          >  2:4531566/1‑150 (MQ=255)
                                                                                                             tGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCGGATATGGCAGTGGGTATCTACGATTGGGCGCTCATTGTCGACCACCAGCGTc                                         >  1:3885928/1‑150 (MQ=255)
                                                                                                               ccATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCGGATATGGCAGTGGGTATCTACGATTGGGCGCTCATTGTCGACCACCAGCGTCAt                                       >  1:2095960/1‑150 (MQ=255)
                                                                                                                cATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCGGATATGGCAGTGGGTATCTACGATTGGGCGCTCATTGTCGACCACCAGCGTCATa                                      <  1:440434/150‑1 (MQ=255)
                                                                                                                cATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCGGATATGGCAGTGGGTATCTACGATTGGGCGCTCATTGTCGACCACCAGCGTCATa                                      >  1:2011912/1‑150 (MQ=255)
                                                                                                                  tttCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTTAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCGGATATGGCAGTGGGTATCTACGATTGGGCGCTCATTGTCGACCACCAGCGTCATACa                                    >  1:1237672/1‑150 (MQ=255)
                                                                                                                   ttCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCGGATATGGCAGTGGGTATCTACGATTGGGCGCTCATTGTCGACCACCAGCGTCATACAg                                   >  1:4872560/1‑150 (MQ=255)
                                                                                                                      aGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCGGATATGGCAGTGGGTATCTACGATTGGGCGCTCATTGTCGACCACCAGCGTCATACAGttt                                <  2:3723035/150‑1 (MQ=255)
                                                                                                                        ggcggcGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCGGATATGGCAGTGGGTATCTACGATTGGGCGCTCATTGTCGACCACCAGCGTCATACAGTTTCt                              <  2:912778/150‑1 (MQ=255)
                                                                                                                         gcggcgCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCGGATATGGCAGTGGGTATCTACGATTGGGCGCTCATTGTCGACCACCAGCGTCATACAGTTTCtt                             <  1:3583228/150‑1 (MQ=255)
                                                                                                                          cggcgCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCGGATATGGCAGTGGGTATCTACGATTGGGCGCTCATTGTCGACCACCAGCGTCATACAGTTTCttt                            <  1:4259593/150‑1 (MQ=255)
                                                                                                                          cggcgCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCGGATATGGCAGTGGGTATCTACGATTGGGCGCTCATTGTCGACCACCAGCGTCATACAGTTTCttt                            <  1:1471437/150‑1 (MQ=255)
                                                                                                                          cggcgCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCGGATATGGCAGTGGGTATCTACGATTGGGCGCTCATTGTCGACCACCAGCGTCATACAGTTTCttt                            <  2:3617914/150‑1 (MQ=255)
                                                                                                                           ggcgCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCGGATATGGCAGTGGGTATCTACGATTGGGCGCTCATTGTCGACCACCAGCGTCATACAGTTTCTTTg                           <  2:1963808/150‑1 (MQ=255)
                                                                                                                           ggcgCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCGGATATGGCAGTGGGTATCTACGATTGGGCGCTCATTGTCGACCACCAGCGTCATACAGTTTCTTTg                           >  1:3492531/1‑150 (MQ=255)
                                                                                                                                aCTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCGGATATGGCAGTGGGTATCTACGATTGGGCGCTCATTGTCGACCACCAGCGTCATACAGTTTCTTTGCTGAg                      >  2:3468160/1‑150 (MQ=255)
                                                                                                                                   ggggTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCATAAATTGCTGAACAAGATATCGTTCTGCCGGATATGGCAGTGGGTATCTACGATTGGGCGCTCATTGTCGACCACCAGCGTCATACAGTTTCTTTGCTGAGTCa                   <  1:3530437/150‑1 (MQ=255)
                                                                                                                                      gttgttTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCGGATATGGCAGTGGGTATCTACGATTGGGCGCTCATTGTCGACCACCAGCGTCATACAGTTTCTTTGCTGAGTCATaa                >  2:4347090/1‑150 (MQ=255)
                                                                                                                                      gttgttTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCGGATATGGCAGTGGGTATCTACGATTGGGCGCTCATTGTCGACCACCAGCGTCATACAGTTTCTTTGCTGAGTCATaa                >  2:2874151/1‑150 (MQ=255)
                                                                                                                                        tgttTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCGGATATGGCAGTGGGTATCTACGATTGGGCGCTCATTGTCGACCACCAGCGTCATACAGTTTCTTTGCTGAGTCATAatg              <  2:3595865/150‑1 (MQ=255)
                                                                                                                                         gttTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCGGATATGGCAGTGGGTATCTACGATTGGGCGCTCATTGTCGACCACCAGCGTCATACAGTTTCTTTGCTGAGTCATAatga             >  1:3718564/1‑150 (MQ=255)
                                                                                                                                          tttGGCTACGTTCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCGGATATGGCAGTGGGTATCTACGATTGGGCGCTCATTGTCGACCACCAGCGTCATACAGTTTCTTTGCTGAGTCATAatgat            <  1:3862543/150‑1 (MQ=255)
                                                                                                                                          tttGGCTACGTTCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCGGATATGGCAGTGGGTATCTACGATTGGGCGCTCATTGTCGACCACCAGCGTCATACAGTTTCTTTGCTGAGTCATAatgat            <  2:3431696/150‑1 (MQ=255)
                                                                                                                                            tGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCGGATATGGCAGTGGGTATCTACGATTGGGCGCTCATTGTCGACCACCAGCGTCATACAGTTTCTTTGCTGAGTCATAATGATGt          >  1:1861462/1‑150 (MQ=255)
                                                                                                                                              gCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCGGATATGGCAGTGGGTATCTACGATTGGGCGCTCATTGTCGACCACCAGCGTCATACAGTTTCTTTGCTGAGTCATAATGATGTCa        <  2:2610083/150‑1 (MQ=255)
                                                                                                                                                tACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCGGATATGGCAGTGGGTATCTACGATTGGGCGCTCATTGTCGACCACCAGCGTCATACAGTTTCTTTGCTGAGTCATAATGATGTCAAt      <  1:4126853/150‑1 (MQ=255)
                                                                                                                                                  cGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCGGATATGGCAGTGGGTATCTACGATTGGGCGCTCATTGTCGACCACCAGCGTCATACAGTTTCTTTGCTGAGTCATAATGATGTCAATGc    <  2:487744/150‑1 (MQ=255)
                                                                                                                                                    aTCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCGGATATGGCAGTGGGTATCTACGATTGGGCGCTCATTGTCGACCACCAGCGTCATACAGTTTCTTTGCTGAGTCATAATGATGTCAATGccc  <  1:663825/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
ACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCGGATATGGCAGTGGGTATCTACGATTGGGCGCTCATTGTCGACCACCAGCGTCATACAGTTTCTTTGCTGAGTCATAATGATGTCAATGCCC  >  CP009273/1889251‑1889548

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: