Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 630,742 C→A 41.3% E31* (GAA→TAA)  ybdM ← Spo0J family protein, ParB‑like nuclease domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273630,7420CA41.3% 348.1 / inf 304E31* (GAA→TAA) ybdMSpo0J family protein, ParB‑like nuclease domain
Reads supporting (aligned to +/- strand):  ref base C (85/93);  new base A (52/74);  total (137/167)
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 309 total aligned reads displayed.

AATGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCTTTTATCCGTTCGTTATAACGCTGATAGTTCTTCGATTTTGTGGG  >  CP009273/630593‑630891
                                                                                                                                                     |                                                                                                                                                     
 aTGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCa                                                                                                                                                      >  1:3213326/1‑150 (MQ=255)
  tGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCAt                                                                                                                                                     >  2:2235035/1‑150 (MQ=255)
       ttttGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTAATTAATc                                                                                                                                                <  1:2644924/150‑1 (MQ=255)
                agTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTAATTAATCGCTTTGATg                                                                                                                                       <  2:3613935/150‑1 (MQ=255)
                 gTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCCTAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGAGTTCATTAATCGCGTTGATGc                                                                                                                                      >  2:2776423/1‑150 (MQ=255)
                 gTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTAATTAATCGCTTTGATGc                                                                                                                                      >  2:1960858/1‑150 (MQ=255)
                    ttttttCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGat                                                                                                                                   <  2:637025/150‑1 (MQ=255)
                      ttttCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGatca                                                                                                                                 <  2:749810/150‑1 (MQ=255)
                      ttttCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTAATTAATCGCTTTGATGCGatca                                                                                                                                 <  1:450688/150‑1 (MQ=255)
                         tCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCt                                                                                                                              <  2:291711/150‑1 (MQ=255)
                         tCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTAATTAATCGCTTTGATGCGATCATCt                                                                                                                              <  2:911510/150‑1 (MQ=255)
                           gggTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTc                                                                                                                            <  1:3135230/150‑1 (MQ=255)
                            ggtggtGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTcc                                                                                                                           <  1:2027822/150‑1 (MQ=255)
                            ggtggtGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTcc                                                                                                                           >  1:2481430/1‑150 (MQ=255)
                            ggtggtGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTAATTAATCGCTTTGATGCGATCATCTTcc                                                                                                                           <  1:4013698/150‑1 (MQ=255)
                             gtggtgCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCg                                                                                                                          <  2:4057457/150‑1 (MQ=255)
                              tggtgCCTCGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTAATTAATCGCTTTGATGCGATCATCTTCCgg                                                                                                                         <  1:3550005/150‑1 (MQ=255)
                                 tgCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCaa                                                                                                                      >  1:1499519/1‑150 (MQ=255)
                                 tgCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTAATTAATCGCTTTGATGCGATCATCTTCCGGCaa                                                                                                                      >  1:823124/1‑150 (MQ=255)
                                 tgCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTAATTAATCGCTTTGATGCGATCATCTTCCGGCaa                                                                                                                      <  1:591773/150‑1 (MQ=255)
                                  gCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCacg                                                                                                                     >  2:122453/1‑148 (MQ=255)
                                    cACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCt                                                                                                                   >  2:3658203/1‑150 (MQ=255)
                                     aCGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTAATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTg                                                                                                                  <  1:3772062/150‑1 (MQ=255)
                                     aCGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTAATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTg                                                                                                                  <  1:471713/150‑1 (MQ=255)
                                      cGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTgg                                                                                                                 <  2:578937/150‑1 (MQ=255)
                                         tgttTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCga                                                                                                              <  1:275531/150‑1 (MQ=255)
                                            ttGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTAATTAATCGCTTTGATGTGATCATCTTCCGGCAAGCTGGCGAGaa                                                                                                           <  1:2455737/150‑1 (MQ=255)
                                              gggTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAAt                                                                                                         <  1:1556843/150‑1 (MQ=255)
                                                  tGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCg                                                                                                     >  1:1154225/1‑150 (MQ=255)
                                                   gTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGt                                                                                                    <  2:4306696/150‑1 (MQ=255)
                                                    tAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGtt                                                                                                   <  2:3124593/150‑1 (MQ=255)
                                                    tAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGtt                                                                                                   <  1:1759130/150‑1 (MQ=255)
                                                     aaTCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTa                                                                                                  <  1:3379341/150‑1 (MQ=255)
                                                     aaTCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTa                                                                                                  >  2:3253298/1‑150 (MQ=255)
                                                       tCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTaaa                                                                                                >  1:2671026/1‑150 (MQ=255)
                                                         aTTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATc                                                                                              <  2:1255611/150‑1 (MQ=255)
                                                          ttAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTAATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATcc                                                                                             <  1:4362712/150‑1 (MQ=255)
                                                           tAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCt                                                                                            >  2:118954/1‑150 (MQ=255)
                                                            aGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTg                                                                                           >  1:2190605/1‑150 (MQ=255)
                                                            aGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTAATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTg                                                                                           >  1:90494/1‑150 (MQ=255)
                                                              gCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCg                                                                                         >  1:154116/1‑150 (MQ=255)
                                                                      tGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCgt                                                                                 <  2:2276755/150‑1 (MQ=255)
                                                                        gCTGTTTTTAACCCACAGAACGCAATCCAATGGCTCTTCACGAAAGTGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCgttg                                                                               >  2:91661/1‑150 (MQ=255)
                                                                         cTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCgttgt                                                                              <  1:1668491/150‑1 (MQ=255)
                                                                          tGTTTTTAACGCACAGAACGCAATCCACTGGCTCTTCACGAATGGTGATAACCTGATTGATCGCCATGTGGATTTAATTAATCGCTTTGATGCGATCATCTTCCGTCAAGCTGGCGATAAATTGCGTTAAATCCTGCGTTAATCgttttt                                                                             <  2:3792690/150‑1 (MQ=255)
                                                                          tGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCgttgtt                                                                             <  2:1299448/150‑1 (MQ=255)
                                                                          tGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTAATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCgttgtt                                                                             <  2:3596129/150‑1 (MQ=255)
                                                                              tttAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTAATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCAt                                                                         <  2:2139961/150‑1 (MQ=255)
                                                                                tAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTAATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCTTTGTTGCATAg                                                                       <  2:3181697/150‑1 (MQ=255)
                                                                                tAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTAATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAg                                                                       <  2:2670268/150‑1 (MQ=255)
                                                                                tAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTAATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAg                                                                       >  1:856202/1‑150 (MQ=255)
                                                                                   cccACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTAt                                                                    <  1:43805/150‑1 (MQ=255)
                                                                                     cacaGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTAATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATcc                                                                  <  2:1857445/150‑1 (MQ=255)
                                                                                      acaGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATccc                                                                 >  1:894680/1‑150 (MQ=255)
                                                                                      acaGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATccc                                                                 >  1:1401556/1‑150 (MQ=255)
                                                                                        aGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATAGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCTGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATccccc                                                               >  2:331354/1‑150 (MQ=255)
                                                                                         gAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCa                                                              <  2:258516/150‑1 (MQ=255)
                                                                                         gAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCa                                                              <  1:4560168/150‑1 (MQ=255)
                                                                                          aaCGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCAt                                                             <  2:4633532/150‑1 (MQ=255)
                                                                                          aaCGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCAt                                                             <  2:1729214/150‑1 (MQ=255)
                                                                                          aaCGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTAATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCAt                                                             <  2:1783654/150‑1 (MQ=255)
                                                                                             gCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTAATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTcc                                                          <  1:1929325/150‑1 (MQ=255)
                                                                                               aaTCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTTTCCCCCTTTCCTg                                                        <  1:210225/150‑1 (MQ=255)
                                                                                                aTCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTAATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGa                                                       >  2:1929508/1‑150 (MQ=255)
                                                                                                 tCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGAc                                                      <  2:1668136/150‑1 (MQ=255)
                                                                                                 tCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGAc                                                      <  1:3333137/150‑1 (MQ=255)
                                                                                                    aCTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGtt                                                   >  1:1283192/1‑150 (MQ=255)
                                                                                                    aCTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTAATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGtt                                                   <  1:635104/150‑1 (MQ=255)
                                                                                                      tGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGtttt                                                 >  1:3021910/1‑150 (MQ=255)
                                                                                                       ggCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTc                                                <  2:4083044/150‑1 (MQ=255)
                                                                                                       ggCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTAATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTc                                                <  1:2409930/150‑1 (MQ=255)
                                                                                                        gCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTAATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTcc                                               >  2:1482401/1‑150 (MQ=255)
                                                                                                           ctTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTAATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCttt                                            >  2:2975681/1‑150 (MQ=255)
                                                                                                            ttCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCtttt                                           <  2:2015198/150‑1 (MQ=255)
                                                                                                             tCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCTTTTa                                          <  1:410816/150‑1 (MQ=255)
                                                                                                               aCGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCTTTTATc                                        <  1:306439/150‑1 (MQ=255)
                                                                                                               aCGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTAATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCTTTTATc                                        >  1:3861677/1‑150 (MQ=255)
                                                                                                                cGAAAGTGGCTAACCTGATTTATCGCCTTGCGTATTTATTTAATCGCTTTGATGCGATCATCTTACTGAAAGCTGGCGAGAAATTGCGTTAAATCCTGAGTTAATAGTTGTTGCATAGTATACCCCATTCCTGACGTTTTCCTTTTATcc                                       <  2:2073432/150‑1 (MQ=255)
                                                                                                                  aaaGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCTTTTATCcgt                                     <  2:4200157/150‑1 (MQ=255)
                                                                                                                  aaaGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCTTTTATCcgt                                     >  1:1904505/1‑150 (MQ=255)
                                                                                                                  aaaGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCTTTTATCcgt                                     <  2:2618123/150‑1 (MQ=255)
                                                                                                                   aaGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCTTTTATCcgtt                                    >  2:3917429/1‑150 (MQ=255)
                                                                                                                     ggggCTAACCTGATGGATCGCCATGCGGATTTAATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCTTTTATCcgttcg                                  <  2:2475775/150‑1 (MQ=255)
                                                                                                                      gggCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCTTTTATCcgttcgt                                 <  2:1523298/150‑1 (MQ=255)
                                                                                                                       ggCTAACCTGATGGATCGCCATGCGGATTTAATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCTTTTATCcgttcgtt                                >  2:4535172/1‑150 (MQ=255)
                                                                                                                       ggCTAACCTGATGGATCGCCATGCGGATTTAATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCTTTTATCcgttcgtt                                <  1:1817476/150‑1 (MQ=255)
                                                                                                                        gCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCTTTTATCCGTTCGTta                               <  1:1032038/150‑1 (MQ=255)
                                                                                                                           aaCCTGATGGATCGCCATGCGGATTTAATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCTTTTATCCGTTCGTTATaa                            <  2:4377763/150‑1 (MQ=255)
                                                                                                                            aCCTGATGGATCGCCATGCGGATTTAATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCTTTTATCCGTTCGTTATAAc                           >  1:3607443/1‑150 (MQ=255)
                                                                                                                                gATGGATCGCCATGCGGATTTAATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCTTTTATCCGTTCGTTATAACGCTg                       >  1:2104141/1‑150 (MQ=255)
                                                                                                                                  tGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCTTTTATCCGTTCGTTATAACGCTGAt                     <  1:2256355/150‑1 (MQ=255)
                                                                                                                                  tGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCTTTTATCCGTTCGTTATAACGCTGAt                     <  1:285650/150‑1 (MQ=255)
                                                                                                                                    gATCGCCATGCGGATTTAATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCTTTTATCCGTTCGTTATAACGCTGATAg                   <  1:883044/150‑1 (MQ=255)
                                                                                                                                      tCGCCATGCGGATTTAATTAATCTCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTCTCCCCCATTCCTGACGTTTTCCTTTTATCCGTGCGTTATAACGCTGATAggt                 >  1:3479766/1‑148 (MQ=255)
                                                                                                                                        gCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCTTTTATCCGTTCGTTATAACGCTGATAGttct               <  1:3455496/150‑1 (MQ=255)
                                                                                                                                          cATGCGGATTTCATTAATCGCGTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCTTTTATCCGTTCGTTATAACGCTGATAGttcttc             <  2:2505705/150‑1 (MQ=255)
                                                                                                                                           aTGCGGATTTAATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCTTTTATCCGTTCGTTATAACGCTGATAGTTCTTCg            <  1:3984513/150‑1 (MQ=255)
                                                                                                                                                gATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCTTTTATCCGTTCGTTATAACGCTGATAGTTCTTCGAtttt       >  2:307546/1‑150 (MQ=255)
                                                                                                                                                     cATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGTCGTTTTCCTTTTATCCGTTCGTTATAACGCTGATAGTTCTTCGATTTTGTggg  <  1:4446828/150‑1 (MQ=255)
                                                                                                                                                     cATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCTTTTATCCGTTCGTTATAACGCTGATAGTTCTTCGATTTTGTggg  <  2:4233076/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AATGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCTTTTATCCGTTCGTTATAACGCTGATAGTTCTTCGATTTTGTGGG  >  CP009273/630593‑630891

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: