Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,268,967 A→T 10.6% S187S (TCA→TCT yjcD → inner membrane putative guanine permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,268,9670AT10.6% 1227.3 / 108.1 378S187S (TCA→TCTyjcDinner membrane putative guanine permease
Reads supporting (aligned to +/- strand):  ref base A (163/175);  new base T (16/24);  total (179/199)
Fisher's exact test for biased strand distribution p-value = 4.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.42e-01

Only 100 of 382 total aligned reads displayed.

CTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTTCTCCGGCGTTTTCGCCATGCCTTCATTGAGCGA  >  CP009273/4268819‑4269116
                                                                                                                                                    |                                                                                                                                                     
 tCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACt                                                                                                                                                     <  1:4303613/150‑1 (MQ=255)
   aCGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTgg                                                                                                                                                   >  2:2117989/1‑150 (MQ=255)
     ggTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTa                                                                                                                                                 <  2:1304270/150‑1 (MQ=255)
        gtgGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGt                                                                                                                                              <  2:2771809/150‑1 (MQ=255)
           gcgcACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTg                                                                                                                                           <  2:3392410/150‑1 (MQ=255)
            cgcACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTgg                                                                                                                                          <  2:785997/150‑1 (MQ=255)
              cacggcacggGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCTCTGGTAGGTCTGGCg                                                                                                                                        >  2:1767322/1‑150 (MQ=255)
                cggcacggGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGt                                                                                                                                      <  1:38671/150‑1 (MQ=255)
                 ggcacggGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTg                                                                                                                                     <  2:2565285/150‑1 (MQ=255)
                    acggGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGatc                                                                                                                                  <  2:4288068/150‑1 (MQ=255)
                       gggATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGatcatc                                                                                                                               <  2:1038387/150‑1 (MQ=255)
                       gggATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGatcatc                                                                                                                               <  2:566341/150‑1 (MQ=255)
                       gggATTGGTATCGGTCTGTTCCTGCTGCTCAGTGCCGCTAACGGTGTCGGTCTGGTGATTAACAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGatcatc                                                                                                                               <  2:3477503/150‑1 (MQ=255)
                        ggATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCg                                                                                                                              <  1:1884013/150‑1 (MQ=255)
                        ggATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCg                                                                                                                              >  1:4548537/1‑150 (MQ=255)
                            tGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCt                                                                                                                          >  2:2884925/1‑150 (MQ=255)
                               tATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGa                                                                                                                       >  1:2503959/1‑150 (MQ=255)
                                  cGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGaaaa                                                                                                                    >  1:2850779/1‑150 (MQ=255)
                                    gTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGCGATCATCGGCCTGGAAAAAc                                                                                                                  <  1:2634894/150‑1 (MQ=255)
                                     tCTGTTACTGCTGCTCATTGCCGCTACCGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCGTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACt                                                                                                                 <  1:3321196/150‑1 (MQ=255)
                                      cTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTg                                                                                                                <  2:4017091/150‑1 (MQ=255)
                                          tCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAg                                                                                                            <  2:4131242/150‑1 (MQ=255)
                                            ctgctgCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTc                                                                                                          >  2:3496011/1‑150 (MQ=255)
                                               ctgctCAGTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCGTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCt                                                                                                       <  1:155437/150‑1 (MQ=255)
                                                tgctCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCtg                                                                                                      <  1:2351018/150‑1 (MQ=255)
                                                   tCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCtggtg                                                                                                   >  1:2412909/1‑150 (MQ=255)
                                                     aTTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGc                                                                                                 <  1:3149456/150‑1 (MQ=255)
                                                       tGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCAt                                                                                               >  1:1316582/1‑150 (MQ=255)
                                                          cGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTct                                                                                            <  2:2789536/150‑1 (MQ=255)
                                                              aaCGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCTCTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTctgctg                                                                                        <  1:780576/150‑1 (MQ=255)
                                                              aaCGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTctgctg                                                                                        <  2:321208/150‑1 (MQ=255)
                                                               aCGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGa                                                                                       <  2:1923735/150‑1 (MQ=255)
                                                                   tgtCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCAt                                                                                   >  1:179210/1‑150 (MQ=255)
                                                                      cGGTCTGTTGATTAAAAACCCGATTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCTCTGGTAGGTCTGGCGGTGATCATCGGCATGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTAt                                                                                <  1:1781054/150‑1 (MQ=255)
                                                                      cGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTAt                                                                                >  1:3554276/1‑150 (MQ=255)
                                                                      cGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTAt                                                                                <  2:1948915/150‑1 (MQ=255)
                                                                       ggTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATc                                                                               <  1:2236892/150‑1 (MQ=255)
                                                                       ggTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATc                                                                               <  2:713537/150‑1 (MQ=255)
                                                                       ggTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATc                                                                               >  1:1199702/1‑150 (MQ=255)
                                                                       ggTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATc                                                                               <  1:2055887/150‑1 (MQ=255)
                                                                       ggTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATc                                                                               >  1:4387584/1‑150 (MQ=255)
                                                                       ggTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATc                                                                               <  2:2324130/150‑1 (MQ=255)
                                                                       ggTCTGGTGATTAAAAACCCGCTAGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGAACATTaac                                                                               <  1:2864977/150‑3 (MQ=255)
                                                                        gTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCg                                                                              <  2:91816/150‑1 (MQ=255)
                                                                          cTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGt                                                                            <  1:482906/150‑1 (MQ=255)
                                                                           tGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTCTGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTa                                                                           <  1:594805/150‑1 (MQ=255)
                                                                            ggTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTAt                                                                          <  1:4666459/150‑1 (MQ=255)
                                                                            ggTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTAt                                                                          <  2:4771534/150‑1 (MQ=255)
                                                                            ggTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTAt                                                                          <  2:1574126/150‑1 (MQ=255)
                                                                               gATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTAtctc                                                                       <  2:4176263/150‑1 (MQ=255)
                                                                                   aaaaaCCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCTCTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAAtt                                                                   >  2:3139984/1‑150 (MQ=255)
                                                                                   aaaaaCCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAAtt                                                                   <  1:1737528/150‑1 (MQ=255)
                                                                                    aaaaCCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCTCTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTg                                                                  <  2:86106/150‑1 (MQ=255)
                                                                                    aaaaCCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTg                                                                  <  2:1852988/150‑1 (MQ=255)
                                                                                      aaCCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTc                                                                >  1:1079044/1‑150 (MQ=255)
                                                                                      aaCCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTc                                                                >  1:1356720/1‑150 (MQ=255)
                                                                                       aCCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCg                                                               <  2:4159111/150‑1 (MQ=255)
                                                                                          cGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGtt                                                            >  1:3767708/1‑150 (MQ=255)
                                                                                            cTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTg                                                          >  1:2705831/1‑150 (MQ=255)
                                                                                            cTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTg                                                          <  1:973709/150‑1 (MQ=255)
                                                                                                aTGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCt                                                      >  1:347882/1‑150 (MQ=255)
                                                                                                aTGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCt                                                      <  1:88413/150‑1 (MQ=255)
                                                                                                  ggTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTc                                                    >  2:2526969/1‑150 (MQ=255)
                                                                                                  ggTCTGCCCGTTGCGCTGGGAGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGACTGGCGGTGATCATCTGCCTGGAAAAACTGAAAGTCCCTGGTGGCATACTGCTGACCATTATCGGTATCTCAATTGTCGGTGTGATCTTc                                                    >  2:170152/1‑150 (MQ=255)
                                                                                                     cTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGAt                                                 >  2:2470111/1‑150 (MQ=255)
                                                                                                      tGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATc                                                <  2:64811/150‑1 (MQ=255)
                                                                                                      tGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATc                                                >  2:3775845/1‑150 (MQ=255)
                                                                                                         ccGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTa                                             <  2:2379879/150‑1 (MQ=255)
                                                                                                          cGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTaa                                            <  2:1261477/150‑1 (MQ=255)
                                                                                                              gcgcTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTc                                        <  2:1532816/150‑1 (MQ=255)
                                                                                                               cgcTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTcc                                       >  1:3542222/1‑150 (MQ=255)
                                                                                                                  tGGGTGATTTCGCGACCTTCCCGGTGATTATGTCTCTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCAtt                                    <  2:74887/150‑1 (MQ=255)
                                                                                                                     gTGATTTCGCGACCTTCCCGGTGATTATGTCTCTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTtct                                 >  2:1992794/1‑150 (MQ=255)
                                                                                                                     gTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTtct                                 <  2:2269133/150‑1 (MQ=255)
                                                                                                                          ttCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTTCTCCGGc                            <  2:1653974/150‑1 (MQ=255)
                                                                                                                           tCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTTCTCCGGCg                           >  2:1061862/1‑150 (MQ=255)
                                                                                                                           tCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTTCTCCGGCg                           >  1:2215735/1‑150 (MQ=255)
                                                                                                                           tCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCAACGGCCTGGAAAAACTGAAAGTCCCTGGTGGCAATCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCAAACGTCCATTTCTCCGGCg                           >  2:15694/1‑150 (MQ=255)
                                                                                                                             gcgACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTTCTCCGGCGtt                         <  1:1414823/150‑1 (MQ=255)
                                                                                                                              cgACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTTCTCCGGCGttt                        <  2:3602894/150‑1 (MQ=255)
                                                                                                                              cgACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTTCTCCGGCGttt                        >  2:44304/1‑150 (MQ=255)
                                                                                                                               gACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTTCTCCGGCGtttt                       >  2:1454657/1‑150 (MQ=255)
                                                                                                                                  cTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTTCTCCGGCGTTTTCGc                    >  2:4658317/1‑150 (MQ=255)
                                                                                                                                  cTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTTCTCCGGCGTTTTCGc                    >  2:4543024/1‑150 (MQ=255)
                                                                                                                                   ttCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTTCTCCGGCGTTTTCGcc                   >  2:3044724/1‑150 (MQ=255)
                                                                                                                                   ttCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTTCTCCGGCGTTTTCGcc                   >  1:3362679/1‑150 (MQ=255)
                                                                                                                                     cccGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTTCTCCGGCGTTTTCGCCAt                 <  1:687040/150‑1 (MQ=255)
                                                                                                                                       cGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTTCTCCGGCGTTTTCGCCATGc               >  2:2637310/1‑150 (MQ=255)
                                                                                                                                        ggTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTTCTCCGGCGTTTTCGCCATGcc              <  1:3091980/150‑1 (MQ=255)
                                                                                                                                        ggTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTTCTCCGGCGTTTTCGCCATGcc              <  1:1047364/150‑1 (MQ=255)
                                                                                                                                          tGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTTCTCCGGCGTTTTCGCCATGCCtt            <  1:2462548/150‑1 (MQ=255)
                                                                                                                                           gATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTTCTCCGGCGTTTTCGCCATGCCTTc           >  2:2972512/1‑150 (MQ=255)
                                                                                                                                            aTTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTTCTCCGGCGTTTTCGCCATGCCTTCa          <  1:2687429/150‑1 (MQ=255)
                                                                                                                                             ttATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTTCTCCGGCGTTTTCGCCATGCCTTCAt         >  2:4409746/1‑150 (MQ=255)
                                                                                                                                              tATGTCACTGGTAGGTCTGGCGGTTATCATCGGCCTTGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTTCTCCGGCGTTTTCGCCATGCCTTCAtt        <  1:4141143/150‑1 (MQ=255)
                                                                                                                                              tATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTTCTCCGGCGTTTTCGCAATGCCTTCAtt        <  1:2265886/150‑1 (MQ=255)
                                                                                                                                               aTGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTTCTCCGGCGTTTTCGCCATGCCTTCATTg       >  2:676507/1‑150 (MQ=255)
                                                                                                                                               aTGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTTCTCCGGCGTTTTCGCCATGCCTTCATTg       >  1:3902190/1‑150 (MQ=255)
                                                                                                                                                   cACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTTCTCCGGCGTTTTCGCCATGCCTTCATTGAGCg   >  2:46351/1‑150 (MQ=255)
                                                                                                                                                    aCTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTTCTCCGGCGTTTTCGCCATGCCTTCATTGAGCGa  <  1:2117989/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTTCTCCGGCGTTTTCGCCATGCCTTCATTGAGCGA  >  CP009273/4268819‑4269116

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: