Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,388,644 T→A 7.3% D325V (GAC→GTC)  yhdP ← DUF3971‑AsmA2 domains protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,388,6440TA 7.3% 1058.9 / 49.8 338D325V (GAC→GTC) yhdPDUF3971‑AsmA2 domains protein
Reads supporting (aligned to +/- strand):  ref base T (126/187);  new base A (13/12);  total (139/199)
Fisher's exact test for biased strand distribution p-value = 2.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 339 total aligned reads displayed.

GCCGCTCGGTTGTGTGGAGCGCCAAACATCACCCAGTGCAGGTGAAAGTTTCGCGGCCAGCGGGCGTATGCCCTCCAGGCCTGCCAGCTCCAGATTACTGGCGCGAATCCGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGTCCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCGGCCAGGGTTTGCCGTCCATCGTGATCCGTGTATCGGGAATAGAGAACTGCCAGCCCGGATTTTCACGCGTAATATGCGCGGTCAGATTATCCACCGACAGCGTATGCGTT  >  CP009273/3388496‑3388793
                                                                                                                                                    |                                                                                                                                                     
  cGCTCGGTTGTGTGGAGCGCCAAACATCACCCAGTGCAGGTGAAAGTTTCGCGGCCAGCGGGCGTATGCCCTCCAGGCCTGCCAGCTCCAGATTACTGGCGCGAATCCGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGTCCt                                                                                                                                                    >  2:3236012/1‑150 (MQ=255)
   gCTCGGTTGTGTGGAGCGCCAAACATCACCCAGTGCAGGTGAAAGTTTCGCGGCCAGCGGGCGTATGCCCTCCAGGCCTGCCAGCTCCAGATTACTGGCGCGAATCCGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGTCCTg                                                                                                                                                   >  2:367980/1‑150 (MQ=255)
   gCTCGGTTGTGTGGAGCGCCAAACATCACCCAGTGCAGGTGAAAGTTTCGCGGCCAGCGGGCGTATGCCCTCCAGGCCTGCCAGCTCCAGATTACTGGCGCGAATCCGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGTCCTg                                                                                                                                                   <  2:1511120/150‑1 (MQ=255)
   gCTCGGTTGTGTGGAGCGCCAAACATCACCCAGTGCAGGTGAAAGTTTCGCGGCCAGCGGGCGTATGCCCTCCAGGCCTGCCAGCTCCAGATTACTGGCGCGAATCCGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGTCCTg                                                                                                                                                   <  1:832134/150‑1 (MQ=255)
     tCGGTTGTGTGGAGCGCCAAACATCACCCAGTGCAGGTGAAAGTTTCGCGGCCAGCGGGCGTATGCCCTCCAGGCCTGCCAGCTCCAGGTTACTGGCGCGAATCCGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGTCCTGtt                                                                                                                                                 >  1:634247/1‑150 (MQ=255)
      cGGTTGTGTGGAGCGCCAAACATCACCCAGTGCAGGTGAAAGTTTCGCGGCCAGCGGGCGTATGCCCTCCAGGCCTGCCAGCTCCAGATTACTGGCGCGAATCCGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGTCCTGTTc                                                                                                                                                <  2:1251866/150‑1 (MQ=255)
        gTTGTGTGGAGCGCCAAACATCACCCAGTGCAGGTGAAAGTTTCGCGGCCAGCGGGCGTATGCCCTCCAGGCCTGCCAGCTCCAGATTACTGGCGCGAATCCGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGTCCTGTTCCg                                                                                                                                              <  2:1815892/150‑1 (MQ=255)
           gtgtgGAGCGCCAAACATCACCCAGTGCAGGTGAAAGTTTCGCGGCCAGCGGGCGTATGCCCTCCAGGCCTGCCAGCTCCAGATTACTGGCGCGAATCCGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGTCCTGTTCCGGTa                                                                                                                                           >  1:4288294/1‑150 (MQ=255)
              tgGAGCGCCAAACATCACCCAGTGCAGGTGAAAGTTTCGCGGCCAGCGGGCGTATGCCCTCCAGGCCTGCCAGCTCCAGATTACTGGCGCGAATCCGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGTCCTGTTCCGGTATcc                                                                                                                                        <  1:1280905/150‑1 (MQ=255)
               ggAGCGCCAAACATCACCCAGTGCAGGTGAAAGTTTCGCGGCCAGCGGGCGTATGCCCTCCAGGCCTGCCAGCTCCAGATTACTGGCGCGAATCCGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGTCCTGTTCCGGTATCCa                                                                                                                                       >  1:1689716/1‑150 (MQ=255)
               ggAGCGCCAAACATCACCCAGTGCAGGTGAAAGTTTCGCGGCCAGCGGGCGTATGCCCTCCAGGCCTGCCAGCTCCAGATTACTGGCGCGAATCCGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGTCCTGTTCCGGTATCCa                                                                                                                                       >  1:2628509/1‑150 (MQ=255)
                gAGCGCCAAACATCACCCAGTGCAGGTGAAAGTTTCGCGGCCAGCGGGCGTATGCCCTCCAGGCCTGCCAGCTCCAGATTACTGGCGCGAATCCGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGTCCTGTTCCGGTATCCAg                                                                                                                                      >  1:3252876/1‑150 (MQ=255)
                gAGCGCCAAACATCACCCAGTGCAGGTGAAAGTTTCGCGGCCAGCGGGCGTATGCCCTCCAGGCCTGCCAGCTCCAGATTACTGGCGCGAATCCGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGTCCTGTTCCGGTATCCAg                                                                                                                                      <  1:2027930/150‑1 (MQ=255)
                    gcCAAACATCACCCAGTGCAGGTGAAAGTTTCGCGGCCAGCGGGCGTATGCCCTCCAGGCCTGCCAGCTCCAGATTACTGGCGCGAATCCGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGTCCTGTTCCGGTATCCAGGCCa                                                                                                                                  >  1:2137514/1‑150 (MQ=255)
                        aaCATCACCCAGTGCAGGTGAAAGTTTCGCGGCCAGCGGGCGTATGCCCTCCAGGCCTGCCAGCTCCAGATTACTGGCGCGAATCCGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGTCCTGTTCCGGTATCCAGGCCAGCGt                                                                                                                              <  1:2284570/150‑1 (MQ=255)
                        aaCATCACCCAGTGCAGGTGAAAGTTTCGCGGCCAGCGGGCGTATGCCCTCCAGGCCTGCCAGCTCCAGATTACTGGCGCGAATCCGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGTCCTGTTCCGGTATCCAGGCCAGCGt                                                                                                                              >  1:2494770/1‑150 (MQ=255)
                          cATCACCCAGTGCAGGTGAAAGTTTCGCGGCCAGCGGGCGTATGCCCTCCAGGCCTGCCAGCTCCAGATTACTGGCGCGAATCCGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGTCCTGTTCCGGTATCCAGGCCAGCGTCa                                                                                                                            <  1:2584563/150‑1 (MQ=255)
                              aCCCAGTGCAGGTGAAAGTTTCGCGGCCAGCGGGCGTATGCCCTCCAGGCCTGCCAGCTCCAGATTACTGGCGCGAATCCGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGTCCTGTTCCGGTATCCAGGCCAGCGTCAATGc                                                                                                                        <  1:4487711/150‑1 (MQ=255)
                                 cAGTGCAGGTGAAAGTTTCGCGGCCAGCGGGCGTATGCCCTCCAGGCCTGCCAGCTCCAGATTACTGGCGCGAATCCGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGTCCTGTTCCGGTATCCAGGCCAGCGTCAATGCTcc                                                                                                                     <  1:3819200/150‑1 (MQ=255)
                                   gTGCAGGTGAAAGTTTCGCGGCCAGCGGGCGTATGCCCTCCAGGCCTGCCAGCTCCAGATTACTGGCGCGAATCCGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGTCCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGc                                                                                                                   >  2:1045933/1‑150 (MQ=255)
                                   gTGCAGGTGAAAGTTTCGCGGCCAGCGGGCGTATGCCCTCCAGGCCTGCCAGCTCCAGATTACTGGCGCGAATCCGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGACCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGc                                                                                                                   >  2:2131967/1‑150 (MQ=255)
                                    tGCAGGTGAAAGTTTCGCGGCCAGCGGGCGTATGCCCTCCAGGCCTGCCAGCTCCAGATTACTGGCGCGAATCCGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGTCCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCt                                                                                                                  <  2:2952041/150‑1 (MQ=255)
                                    tGCAGGTGAAAGTTTCGCGGCCAGCGGGCGTATGCCCTCCAGGCCTGCCAGCTCCAGATTACTGGCGCGAATCCGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGTCCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCt                                                                                                                  >  1:4813347/1‑150 (MQ=255)
                                    tGCAGGTGAAAGTTTCGCGGCCAGCGGGCGTATGCCCTCCAGGCCTGCCAGCTCCAGATTACTGGCGCGAATCCGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGTCCTGTTCCGGTATCCAGGCAAGCGTCAATGCTCCGCt                                                                                                                  <  1:3721857/150‑1 (MQ=255)
                                      cAGGTGAAAGTTTCGCGGCCAGCGGGCGTATGCCCTCCAGGCCTGCCAGCTCCAGATTACTGGCGCGAATCCGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGTCCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCg                                                                                                                <  1:4145885/150‑1 (MQ=255)
                                       aGGTGAAAGTTTCGCGGCCAGCGGGCGTATGCCCTCCAGGCCTGCCAGCTCCAGATTACTGGCGCGAATCCGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGTCCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCgg                                                                                                               <  2:268746/150‑1 (MQ=255)
                                          tGAAAGTTTCGCGGCCAGCGGGCGTATGCCCTCCAGGCCTGCCAGCTCCAGATTACTGGCGCGAATCCGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGTCCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCGGCCa                                                                                                            >  2:4629373/1‑150 (MQ=255)
                                           gAAAGTTTCGCGGCCAGCGGGCGTATGCCCTCCAGGCCTGCCAGCTCCAGATTACTGGCGCGAATCCGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGTCCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCGGCCAg                                                                                                           <  2:2010133/150‑1 (MQ=255)
                                             aaGTTTCGCGGCCAGCGGGCGTATGCCCTCCAGGCCTGCCAGCTCCAGATTACTGGCGCGAATCCGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGTCCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCGGCCAggg                                                                                                         <  2:474780/150‑1 (MQ=255)
                                              aGTTTCGCGGCCAGCGGGCGTATGCCCTCCAGGCCTGCCAGCTCCAGATTACTGGCGCGAATCCGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGTCCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCGGCCAGGGt                                                                                                        >  2:4441403/1‑150 (MQ=255)
                                                tttCGCGGCCAGCGGGCGTATGCCCTCCAGGCCTGCCAGCTCCAGATTACTGGCGCGAATCCGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGTCCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCGGCCAGGGttt                                                                                                      <  1:2464572/150‑1 (MQ=255)
                                                tttCGCGGCCAGCGGGCGTATGCCCTCCAGGCCTGCCAGCTCCAGATTACTGGCGCGAATCCGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGACCTGCTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCGGCCAGGGttt                                                                                                      <  1:3219905/150‑1 (MQ=255)
                                                 ttCGCGGCCAGCGGGCGTATGCCCTCCAGGCCTGCCAGCTCCAGATTACTGGCGCGAATCCGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGTCCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCGGCCAGGGTTTg                                                                                                     <  1:2933635/150‑1 (MQ=255)
                                                 ttCGCGGCCAGCGGGCGTATGCCCTCCAGGCCTGCCAGCTCCAGATTACTGGCGCGAATCCGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGTCCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCGGCCAGGGTTTg                                                                                                     >  2:583183/1‑150 (MQ=255)
                                                  tCGCGGCCAGCGGGCGTATGCCCTCCAGGCCTGCCAGCTCCAGATTACTGGCGCGAATCCGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGTCCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCGGCCAGGGTTTGc                                                                                                    <  1:472630/150‑1 (MQ=255)
                                                   cgcgGCCAGCGGGCGTATGCCCTCCAGGCCTGCCAGCTCGTGATTACTGGCGCGAATCCGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGTCCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCGGCCAGGGTTTGcc                                                                                                   <  1:4682137/150‑1 (MQ=255)
                                                    gcgGCCAGCGGGCGTATGCCCTCCAGGCCTGCCAGCTCCAGATTACTGGCGCGAATCCGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGTCCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCGGCCAGGGTTTGCCg                                                                                                  <  1:3549764/150‑1 (MQ=255)
                                                    gcgGCCAGCGGGCGTATGCCCTCCAGGCCTGCCAGCTCCAGATTACTGGCGCGAATCCGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGACCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCGGCCAGGGTTTGCCg                                                                                                  >  2:4978676/1‑150 (MQ=255)
                                                    gcgGCCAGCGGGCGTATGCCCTCCAGGCCTGCCAGCTCCAGATTACTGGCGCGAATCCGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGACCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCGGCCAGGGTTTGCCg                                                                                                  >  2:1379260/1‑150 (MQ=255)
                                                     cgGCCAGCGGGCGTATGCCCTCCAGGCCTGCCAGCTCCAGATTACTGGCGCGAATCCGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGTCCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCGGCCAGGGTTTGCCGt                                                                                                 >  2:2007202/1‑150 (MQ=255)
                                                     cgGCCAGCGGGCGTATGCCCTCCAGGCCTGCCAGCTCCAGATTACTGGCGCGAATCCGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGTCCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCGGCCAGGGTTTGCCGt                                                                                                 >  2:1032616/1‑150 (MQ=255)
                                                     cgGCCAGCGGGCGTATGCCCTCCAGGCCTGCCAGCTCCAGATTACTGGCGCGAATCCGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGTCCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCGGCCAGGGTTTGCCGt                                                                                                 >  2:715133/1‑150 (MQ=255)
                                                         cAGCGGGCGTATGCCCTCCAGGCCTGCCAGCTCCAGATTACTGGCGCGAATCCGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGTCCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCGGCCAGGGTTTGCCGTCCAt                                                                                             <  1:4006727/150‑1 (MQ=255)
                                                         cAGCGGGCGTATGCCCTCCAGGCCTGCCAGCTCCAGATTACTGGCGCGAATCCGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGTCCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCGGCCAGGGTTTGCCGTCCAt                                                                                             <  1:2118157/150‑1 (MQ=255)
                                                          aGCGGGCGTATGCCCTCCAGGCCTGCCAGCTCCAGATTACTGGCGCGAATCCGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGACCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCGGCCAGGGTTTGCCGTCCATc                                                                                            >  1:4105613/1‑150 (MQ=255)
                                                           gCGGGCGTATGCCCTCCAGGCCTGCCAGCTCCAGATTACTGGCGCGAATCCGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGTCCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCGGCCAGGGTTTGCCGTCCATCg                                                                                           <  1:55799/150‑1 (MQ=255)
                                                            cGGGCGTATGCCCTCCAGGCCTGCCAGCTCCAGATTACTGGCGCGAATCCGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGTCCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCGGCCAGGGTTTGCCGTCCATCGt                                                                                          <  1:4964596/150‑1 (MQ=255)
                                                             gggCGTATGCCCTCCAGGCCTGCCAGCTCCAGATTACTGGCGCGAATCCGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGTCCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCGGCCAGGGTTTGCCGTCCATCGTg                                                                                         >  2:1860759/1‑150 (MQ=255)
                                                             gggCGTATGCCCTCCAGGCCTGCCAGCTCCAGATTACTGGCGCGAATCCGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGTCCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCGGCCAGGGTTTGCCGTCCATCGTg                                                                                         >  1:3351998/1‑150 (MQ=255)
                                                             gggCGTATGCCCTCCAGGCCTGCCAGCTCCAGATTACTGGCGCGAATCCGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGTCCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCGGCCAGGGTTTGCCGTCCATCGTg                                                                                         >  1:1662121/1‑150 (MQ=255)
                                                               gCGTATGCCCTCCAGGCCTGCCAGCTCCAGATTACTGGCGCGAATCCGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGTCCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCGGCCAGGGTTTGCCGTCCATCGTGAt                                                                                       <  1:4380796/150‑1 (MQ=255)
                                                                  tATGCCCTCCAGGCCTGCCAGCTCCAGATTACTGGCGCGAATCCGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGTCCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCGGCCAGGGTTTGCCGTCCATCGTGATCCg                                                                                    <  1:1865221/150‑1 (MQ=255)
                                                                   aTGCCCTCCAGGCCTGCCAGCTCCAGATTACTGGCGCGAATCCGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGTCCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCGGCCAGGGTTTGCCGTCCATCGTGATCCgt                                                                                   <  2:4597127/150‑1 (MQ=255)
                                                                   aTGCCCTCCAGGCCTGCCAGCTCCAGATTACTGGCGCGAATCCGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGTCCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCGGCCAGGGTTTGCCGTCCATCGTGATCCgt                                                                                   <  1:1050945/150‑1 (MQ=255)
                                                                    tGCCCTCCAGGCCTGCCAGCTCCAGATTACTGGCGCGAATCCGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGTCCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCGGCCAGGGTTTGCCGTCCATCGTGATCCgtg                                                                                  >  2:955728/1‑150 (MQ=255)
                                                                        cTCCAGGCCTGCCAGCTCCAGATTACTGGCGCGAATCCGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGTCCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCGGCCAGGGTTTGCCGTCCATCGTGATCCGTGTATc                                                                              >  1:501543/1‑150 (MQ=255)
                                                                           cAGGCCTGCCAGCTCCAGATTACTGGCGCGAATCCGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGTCCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCGGCCAGGGTTTGCCGTCCATCGTGATCCGTGTATCggg                                                                           <  2:339428/150‑1 (MQ=255)
                                                                            aGGGCTGCCAGCTCCAGATTACTGGAGCTAATCCGGAGTTCGTCACTGCGTATATTGTATTTGCCGAAAAAGTCCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCTGCCAGGGTGTTCCGTCCATCGTGATCAGTGTATCgggc                                                                          >  1:1748645/1‑149 (MQ=255)
                                                                            aGGCCTGCCAGCTCCAGATTACTGGCGCGAATCCGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGTCCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCGGCCAGGGTTTGCCGTCCATCGTGATCCGTGTATCGGGa                                                                          <  1:955728/150‑1 (MQ=255)
                                                                            aGGCCTGCCAGCTCCAGATTACTGGCGCGAATCCGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGACCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCGGCCAGGGTTTGCCGTCCATCGTGATCCGTGTATCGGGa                                                                          <  1:283537/150‑1 (MQ=255)
                                                                              gCCTGCCAGCTCCAGATTACTGGCGCGAATCCGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGTCCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCGGCCAGGGTTTGCCGTCCATCGTGATCCGTGTATCGGGAAt                                                                        <  1:2160781/150‑1 (MQ=255)
                                                                                cTGCCAGCTCCAGATTACTGGCGCGAATCCGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGTCCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCGGCCAGGGTTTGCCGTCCATCGTGATCCGTGTATCGGGAATag                                                                      >  2:22110/1‑150 (MQ=255)
                                                                                 tGCCAGCTCCAGATTACTGGCGCGAATCCGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGTCCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCGGCCAGGGTTTGCCGTCCATCGTGATCCGTGTATCGGGAATaga                                                                     <  2:331356/150‑1 (MQ=255)
                                                                                   ccAGCTCCAGATTACTGGCGCGAATCCGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGTCCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCGGCCAGGGTTTGCCGTCCATCGTGATCCGTGTATCGGGAATAgaga                                                                   <  2:2325654/150‑1 (MQ=255)
                                                                                     aGCTCCAGGTTACTGGCGCGAATCCGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGTCCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCGGCCAGGGTTTGCCGTCCATCGTGATCCGTGTATCGGGAATAGAGAAc                                                                 <  2:634247/150‑1 (MQ=255)
                                                                                         ccAGATTACTGGCGCGAATCCGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGTCCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCGGCCAGGGTTTGCCGTCCATCGTGATCCGTGTATCGGGAATAGAGAACTGcc                                                             <  2:376235/150‑1 (MQ=255)
                                                                                           aGATTACTGGCGCGAATCCGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGTCCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCGGCCAGGGTTTGCCGTCCATCGTGATCCGTGTATCGGGAATAGAGAACTGCCAg                                                           >  1:4072803/1‑150 (MQ=255)
                                                                                              ttACTGGCGCGAATCCGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGTCCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCGGCCAGGGTTTGCCGTCCATCGTGATCCGTGTATCGGGAATAGAGAACTGCCAGccc                                                        <  2:2599820/150‑1 (MQ=255)
                                                                                              ttACTGGCGCGAATCCGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGTCCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCGGCCAGGGTTTGCCGTCCATCGTGATCCGTGTATCGGGAATAGAGAACTGCCAGccc                                                        <  2:3951361/150‑1 (MQ=255)
                                                                                               tACTGGCGCGAATCCGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGTCCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCGGCCAGGGTTTGCCGTCCATCGTGATCCGTGTATCGGGAATAGAGAACTGCCAGCCCg                                                       <  1:4092514/150‑1 (MQ=255)
                                                                                                aCTGGCGCGAATCCGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGTCCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCGGCCAGGGTTTGCCGTCCATCGTGATCCGTGTATCGGGAATAGAGAACTGCCAGCCCgg                                                      <  2:1662236/150‑1 (MQ=255)
                                                                                                aCTGGCGCGAATCCGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGTCCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCGGCCAGGGTTTGCCGTCCATCGTGATCCGTGTATCGGGAATAGAGAACTGCCAGCCCgg                                                      >  1:3058493/1‑150 (MQ=255)
                                                                                                aCTGGCGCGAATCCGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGTCCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCGGCCAGGGTTTGCCGTCCATCGTGATCCGTGTATCGGGAATAGAGAACTGCCAGCCCgg                                                      >  1:1478162/1‑150 (MQ=255)
                                                                                                    gcgcGAATCCGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGTCCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCGGCCAGGGTTTGCCGTCCATCGTGATCCGTGTATCGGGAATAGAGAACTGCCAGCCCGGAttt                                                  <  1:162139/150‑1 (MQ=255)
                                                                                                    gcgcGAATCCGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGTCCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCGGCCAGGGTTTGCCGTCCATCGTGATCCGTGTATCGGGAATAGAGAACTGCCAGCCCGGAttt                                                  <  2:4774228/150‑1 (MQ=255)
                                                                                                    gcgcGAATCCGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGTCCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCGGCCAGGGTTTGCCGTCCATCGTGATCCGTGTATCGGGAATAGAGAACTGCCAGCCCGGAttt                                                  <  2:107982/150‑1 (MQ=255)
                                                                                                     cgcgAATCCGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGTCCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCGGCCAGGGTTTGCCGTCCATCGTGATCCGTGTATCGGGAATAGAGAACTGCCAGCCCGGAtttt                                                 <  2:4160910/150‑1 (MQ=255)
                                                                                                     cgcgAATCCGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGTCCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCGGCCAGGGTTTGCCGTCCATCGTGATCCGTGTATCGGGAATAGAGAACTGCCAGCCCGGAtttt                                                 >  2:5064112/1‑150 (MQ=255)
                                                                                                       cgAATCCGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGTCCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCGGCCAGGGTTTGCCGTCCATCGTGATCCGTGTATCGGGAATAGAGAACTGCCAGCCCGGATTTTCa                                               >  1:4772933/1‑150 (MQ=255)
                                                                                                         aaTCCGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGTCCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCGGCCAGGGTTTGCCGTCCATCGTGATCCGTGTATCGGGAATAGAGAACTGCCAGCCCGGATTTTCAcg                                             <  2:4552323/150‑1 (MQ=255)
                                                                                                          aTCCGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGTCCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCGGCCAGGGTTTGCCGTCCATCGTGATCCGTGTATCGGGAATAGAGAACTGCCAGCCCGGATTTTCAcgc                                            >  2:1139912/1‑150 (MQ=255)
                                                                                                            ccGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGTCCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCGGCCAGGGTTTGCCGTCCATCGTGATCCGTGTATCGGGAATAGAGAACTGCCAGCCCGGATTTTCACGCGt                                          >  2:888343/1‑150 (MQ=255)
                                                                                                             cGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGTCCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCGGCCAGGGTTTGCCGTCCATCGTGATCCGTGTATCGGGAATAGAGAACTGCCAGCCCGGATTTTCACGCGTa                                         <  1:3630738/150‑1 (MQ=255)
                                                                                                             cGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGACCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCGGCCAGGGTTTGCCGTCCATCGTGATCCGTGTATCGGGAATAGAGAACTGCCAGCCCGGATTTTCACGCGTa                                         >  1:4944520/1‑150 (MQ=255)
                                                                                                              ggAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGTCCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCGGCCAGGGTTTGCCGTCCATCGTGATCCGTGTATCGGGAATAGAGAACTGCCAGCCCGGATTTTCACGCGTaa                                        <  2:2780383/150‑1 (MQ=255)
                                                                                                                    cGTCACTGCGTTTATTGTCTTTGCCGCCAACGTCCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCGGCCAGGGTTTGCCGTCCATCGTGATCCGTGTATCGGGAATAGAGAACTGCCAGCCCGGATTTTCACGCGTAATATgcg                                  <  2:289304/150‑1 (MQ=255)
                                                                                                                      tCACTGCGTTTATTGTCTTTGCCGCCAACGTCCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCGGCCAGGGTTTGCCGTCCATCGTGATCCGTGTATCGGGAATAGAGAACTGCCAGCCCGGATTTTCACGCGTAATATGcgcg                                <  2:4691354/150‑1 (MQ=255)
                                                                                                                        aCTGCGTTTATTGTCTTTGCCGCCAACGTCCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCGGCCAGGGTTTGCCGTCCATCGTGATCCGTGTATCGGGAATAGAGAACTGCCAGCCCGGATTTTCACGCGTAATATGCGCGGt                              <  2:3077380/150‑1 (MQ=255)
                                                                                                                          tGCGTTTATTGTCTTTGCCGCCAACGTCCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCGGCCAGGGTTTGCCGTCCATCGTGATCCGTGTATCGGGAATAGAGAACTGCCAGCCCGGATTTTCACGCGTAATATGCGCGGTCa                            <  1:3278584/150‑1 (MQ=255)
                                                                                                                            cGTTTATTGTCTTTGCCGCCAACGTCCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCGGCCAGGGTTTGCCGTCCATCGTGATCCGTGTATCGGGAATAGAGAACTGCCAGCCCGGATTTTCACGCGTAATATGCGCGGTCAGa                          <  2:1846528/150‑1 (MQ=255)
                                                                                                                            cGTTTATTGTCTTTGCCGCCAACGACCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCGGCCAGGGTTTGCCGTCCATCGTGATCCGTGTATCGGGAATAGAGAACTGCCAGCCCGGATTTTCACGCGTAATATGCGCGGTCAGa                          <  2:4104519/150‑1 (MQ=255)
                                                                                                                             gTTTATTGTCTTTGCCGCCAACGTCCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCGGCCAGGGTTTGCCGTCCATCGTGATCCGTGTATCGGGAATAGAGAACTGCCAGCCCGGATTTTCACGCGTAATATGCGCGGTCAGAt                         <  1:300395/150‑1 (MQ=255)
                                                                                                                                   tGTCTTTGCCGCCAACGTCCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCGGCCAGGGTTTGCCGTCCATCGTGATCCGTGTATCGGGAATAGAGAACTGCCAGCCCGGATTTTCACGCGTAATATGCGCGGTCAGATTATCCa                   >  1:1465761/1‑150 (MQ=255)
                                                                                                                                    gTCTTTGCCGCCAACGTCCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCGGCCAGGGTTTGCCGTCCATCGTGATCCGTGTATCGGGAATAGAGAACTGCCAGCCCGGATTTTCACGCGTAATATGCGCGGTCAGATTATCCAc                  >  2:1636310/1‑150 (MQ=255)
                                                                                                                                     tCTTTGCCGCCAACGTCCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCGGCCAGGGTTTGCCGTCCATCGTGATCCGTGTATCGGGAATAGAGAACTGCCCGCCCGGATTTTCACGCGTAATATGCGCGGTCAGATTATCCAcc                 >  1:1178176/1‑150 (MQ=255)
                                                                                                                                     tCTTTGCCGCCAACGTCCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCGGCCAGGGTTTGCCGTCCATCGTGATCCGTGTATCGGGAATAGAGAACTGCCAGCCCGGATTTTCACGCGTAATATGCGCGGTCAGATTATCCAcc                 <  1:590844/150‑1 (MQ=255)
                                                                                                                                          gccgccAACGTCCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCGGCCAGGGTTTGCCGTCCATCGTGATCCGTGTATCGGGAATAGAGAACTGCCAGCCCGGATTTTCACGCGTAATATGCGCGGTCAGATTATCCACCGACAg            <  1:4847312/150‑1 (MQ=255)
                                                                                                                                                aaCGTCCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCGGCCAGGGTTTGCCGTCCATCGTGATCCGTGTATCGGGAATAGAGAACTGCCAGCCCGGATTTTCACGCGTAATATGCGCGGTCAGATTATCCACCGACAGCGTATg      >  1:5106450/1‑150 (MQ=255)
                                                                                                                                                   gTCCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCGGCCAGGGTTTGCCGTCCATCGTGATCCGTGTATCGGGAATAGAGAACTGCCAGCCCGGATTTTCACGCGTAATATGCGCGGTCAGATTATCCACCGACAGCGTATGCGt   >  2:1166910/1‑150 (MQ=255)
                                                                                                                                                   gTCCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCGGCCAGGGTTTGCCGTCCATCGTGATCCGTGTATCGGGAATAGAGAACTGCCAGCCCGGATTTTCACGCGTAATATGCGCGGTCAGATTATCCACCGACAGCGTATGCGt   >  1:1851233/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCCGCTCGGTTGTGTGGAGCGCCAAACATCACCCAGTGCAGGTGAAAGTTTCGCGGCCAGCGGGCGTATGCCCTCCAGGCCTGCCAGCTCCAGATTACTGGCGCGAATCCGGAGTTCGTCACTGCGTTTATTGTCTTTGCCGCCAACGTCCTGTTCCGGTATCCAGGCCAGCGTCAATGCTCCGCTCGGCCAGGGTTTGCCGTCCATCGTGATCCGTGTATCGGGAATAGAGAACTGCCAGCCCGGATTTTCACGCGTAATATGCGCGGTCAGATTATCCACCGACAGCGTATGCGTT  >  CP009273/3388496‑3388793

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: ‑