Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 968,415 G→T 100% S148* (TCG→TAG)  ycbC ← DUF218 superfamily protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273968,4150GT97.0% 1003.7 / 18.7 302S148* (TCG→TAG) ycbCDUF218 superfamily protein
Reads supporting (aligned to +/- strand):  ref base G (7/2);  new base T (145/147);  total (153/149)
Fisher's exact test for biased strand distribution p-value = 1.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 302 total aligned reads displayed.

AAAAAGATCATCGCCCGCGGCAGATGTGATGCGGAGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCGATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCTTCATTCAGGCGTGGCAGACTGTTATTGATTAAATTA  >  CP009273/968268‑968561
                                                                                                                                                   |                                                                                                                                                  
 aaaaGATCATCGCCCGCGGCAGATGTGATGCGGAGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTAtt                                                                                                                                                 >  2:1212575/1‑150 (MQ=255)
  aaaGATCATCGCCCGCGGCAGATGTGATGCGGAGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTg                                                                                                                                                >  2:4577700/1‑150 (MQ=255)
       tcatcGCCCGCGGCAGATGTGATGCGGAGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCaa                                                                                                                                           >  2:227242/1‑150 (MQ=255)
        catcGCCCGCGGCAGATGTGATGCGGAGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAAc                                                                                                                                          >  1:1383277/1‑150 (MQ=255)
         atcGCCCGCGGCAGATGTGATGCGGAGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAAct                                                                                                                                         >  1:3659369/1‑150 (MQ=255)
           cGCCCGCGGCAGATGTGATGCGGAGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCGATTGCGCAActct                                                                                                                                       >  1:3416403/1‑150 (MQ=255)
            gCCCGCGGCAGATGTGATGCGGAGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTg                                                                                                                                      <  1:165703/150‑1 (MQ=255)
             cccGCGGCAGATGTGATGCGGAGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTgg                                                                                                                                     >  2:208637/1‑150 (MQ=255)
              ccGCGGCAGATGTGATGCGGAGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGc                                                                                                                                    <  1:84192/150‑1 (MQ=255)
               cgcgGCAGATGTGATGCGGAGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCa                                                                                                                                   <  1:3481321/150‑1 (MQ=255)
               cgcgGCAGATGTGATGCGGAGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCa                                                                                                                                   <  2:3879990/150‑1 (MQ=255)
                  ggCAGATGTGATGCGGAGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCt                                                                                                                                <  1:4298144/150‑1 (MQ=255)
                       aTGTGATGCGGAGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTc                                                                                                                           <  1:1518832/150‑1 (MQ=255)
                        tgtgATGCGGAGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTcc                                                                                                                          <  1:2019062/150‑1 (MQ=255)
                            aTGCGGAGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGTTATCTTTGGAAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGt                                                                                                                     >  2:1864900/1‑150 (MQ=255)
                            aTGCGGAGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTg                                                                                                                      >  2:4491762/1‑150 (MQ=255)
                            aTGCGGAGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTg                                                                                                                      <  2:1380498/150‑1 (MQ=255)
                              gCGGAGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTa                                                                                                                    >  1:1166054/1‑150 (MQ=255)
                              gCGGAGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCGATTGCGCAACTCTGGCACCTACTTCCGCTGTa                                                                                                                    <  1:2403576/150‑1 (MQ=255)
                               cGGAGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTAc                                                                                                                   <  1:2705744/150‑1 (MQ=255)
                                  aGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCa                                                                                                                >  2:1907373/1‑150 (MQ=255)
                                   ggTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCAc                                                                                                               >  1:3285067/1‑150 (MQ=255)
                                   ggTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCAc                                                                                                               <  2:1271480/150‑1 (MQ=255)
                                   ggTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCAc                                                                                                               >  1:142320/1‑150 (MQ=255)
                                      cACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGt                                                                                                            >  1:2779155/1‑150 (MQ=255)
                                       aCCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTa                                                                                                           >  1:3297917/1‑150 (MQ=255)
                                        ccAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTAt                                                                                                          <  2:525044/150‑1 (MQ=255)
                                         cagcagGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTAtt                                                                                                         >  1:2819558/1‑150 (MQ=255)
                                            cagGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGt                                                                                                      <  2:715388/150‑1 (MQ=255)
                                             agGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGtt                                                                                                     >  1:3163124/1‑150 (MQ=255)
                                               gAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCATCTTCTTCGGTATCTTTAGGCAAATACAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTAGTGCGCAACTCAGGCACCTACTTCCGCTGTACTCACCGTATTGGtttt                                                                                                   >  1:73060/1‑150 (MQ=255)
                                                     ggCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTAc                                                                                             <  1:2918642/150‑1 (MQ=255)
                                                       cGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGc                                                                                           >  1:3382400/1‑150 (MQ=255)
                                                         tCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCt                                                                                         >  1:1082409/1‑150 (MQ=255)
                                                         tCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCt                                                                                         <  2:1313960/150‑1 (MQ=255)
                                                          cACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTc                                                                                        <  2:2195860/150‑1 (MQ=255)
                                                           aCCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTcc                                                                                       <  2:397875/150‑1 (MQ=255)
                                                           aCCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTcc                                                                                       >  2:1491872/1‑150 (MQ=255)
                                                           aCCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTcc                                                                                       >  2:1135157/1‑150 (MQ=255)
                                                            ccAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTccc                                                                                      >  2:4312707/1‑150 (MQ=255)
                                                               aTCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTg                                                                                   <  2:472273/150‑1 (MQ=255)
                                                                tCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGa                                                                                  <  2:111888/150‑1 (MQ=255)
                                                                  gCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAg                                                                                >  1:4556013/1‑150 (MQ=255)
                                                                   ccTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGa                                                                               >  2:3741991/1‑150 (MQ=255)
                                                                     tGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGTTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATc                                                                             <  2:910594/150‑1 (MQ=255)
                                                                      gTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCa                                                                            <  1:185611/150‑1 (MQ=255)
                                                                         tCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGtt                                                                         <  2:2540755/150‑1 (MQ=255)
                                                                           aCTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGtttt                                                                       <  1:3944903/150‑1 (MQ=255)
                                                                           aCTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGtttt                                                                       <  2:2825666/150‑1 (MQ=255)
                                                                           aCTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGtttt                                                                       >  2:179687/1‑150 (MQ=255)
                                                                            ctgctgCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTg                                                                      <  1:3838747/150‑1 (MQ=255)
                                                                                tgcAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATcc                                                                  >  1:4447782/1‑150 (MQ=255)
                                                                                   aGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCgg                                                               >  1:1499025/1‑150 (MQ=255)
                                                                                   aGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCgg                                                               <  2:1145753/150‑1 (MQ=255)
                                                                                   aGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCgg                                                               <  1:2315327/150‑1 (MQ=255)
                                                                                    gCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCACCTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGTTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGa                                                              <  2:3295012/150‑1 (MQ=255)
                                                                                    gCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGa                                                              <  2:3776700/150‑1 (MQ=255)
                                                                                       tcttcttcGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGAttt                                                           <  1:2448987/150‑1 (MQ=255)
                                                                                        cttcttcGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTCTTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATGGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGAtttt                                                          >  2:1532917/1‑150 (MQ=255)
                                                                                        cttcttcGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGAtttt                                                          <  2:2442676/150‑1 (MQ=255)
                                                                                         ttcttcGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTc                                                         <  1:4177357/150‑1 (MQ=255)
                                                                                          tcttcGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTcg                                                        <  2:3623228/150‑1 (MQ=255)
                                                                                          tcttcGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTcg                                                        <  1:2260205/150‑1 (MQ=255)
                                                                                           cttcGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTcgc                                                       >  2:1278033/1‑150 (MQ=255)
                                                                                               ggTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCa                                                   <  1:3118073/150‑1 (MQ=255)
                                                                                               ggTAACTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGAGGCACGCCCAGCTATTTCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTAGCTACGACTCCCGAGAAGATAAGTTGTGATCCCGGATATACGCGCCa                                                   >  1:1684848/1‑150 (MQ=255)
                                                                                                gTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCAt                                                  <  1:528900/150‑1 (MQ=255)
                                                                                                 tATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATa                                                 >  2:4061486/1‑150 (MQ=255)
                                                                                                 tATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGCCACGCCCCGCAAATGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATAGGTTTTTGCTACGCCTCCCGTCAAGATCACTTTTGATCCCGGATTATCGCGCCaaa                                                 >  1:2876336/1‑148 (MQ=255)
                                                                                                    cTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTGTTGATCCCGGATTTTCGCGCCATAAGc                                              >  1:3018512/1‑150 (MQ=255)
                                                                                                     ttttGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCg                                             <  2:1593912/150‑1 (MQ=255)
                                                                                                       ttGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGaa                                           >  1:4013738/1‑150 (MQ=255)
                                                                                                        tGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAAt                                          >  2:2392884/1‑150 (MQ=255)
                                                                                                        tGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGATATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAAt                                          <  2:2732483/150‑1 (MQ=255)
                                                                                                           cAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATAcc                                       >  2:2672945/1‑150 (MQ=255)
                                                                                                               tCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCACCTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCttca                                   <  2:2855630/150‑1 (MQ=255)
                                                                                                               tCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCttca                                   >  2:461330/1‑150 (MQ=255)
                                                                                                                 cAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCttcatt                                 <  2:4398505/150‑1 (MQ=255)
                                                                                                                 cAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCttcatt                                 >  1:3181408/1‑150 (MQ=255)
                                                                                                                     gTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCTTCATTCAgg                             >  2:3414860/1‑150 (MQ=255)
                                                                                                                       gATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCTTCATTCAGGCg                           <  1:2556160/150‑1 (MQ=255)
                                                                                                                         tAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCTTCATTCAGGCGTg                         >  1:3753544/1‑150 (MQ=255)
                                                                                                                         tAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCTTCATTCAGGCGTg                         <  2:1697506/150‑1 (MQ=255)
                                                                                                                          aaTTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCTTCATTCAGGCGTgg                        >  1:2637522/1‑150 (MQ=255)
                                                                                                                          aaTTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCTTCATTCAGGCGTgg                        <  1:1011545/150‑1 (MQ=255)
                                                                                                                          aaTTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCTTCATTCAGGCGTgg                        >  1:2197553/1‑150 (MQ=255)
                                                                                                                          aaTTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCTTCATTCAGGCGTgg                        >  1:2229645/1‑150 (MQ=255)
                                                                                                                            tttGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCTTCATTCAGGCGTGGCa                      <  1:1589152/150‑1 (MQ=255)
                                                                                                                             ttGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCTTCATTCAGGCGTGGCAg                     >  2:22522/1‑150 (MQ=255)
                                                                                                                               gCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCTTCATTCAGGCGTGGCAGAc                   <  2:2819447/150‑1 (MQ=255)
                                                                                                                               gCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCTTCATTCAGGCGTGGCAGAc                   <  2:3217828/150‑1 (MQ=255)
                                                                                                                                 tCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCTTCATTCAGGCGTGGCAGACTg                 <  1:1297189/150‑1 (MQ=255)
                                                                                                                                     gcgGCACGCCCAGCGATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCTTCATTCAGGCGTGGCAGACTGTTAt             >  1:1654134/1‑150 (MQ=255)
                                                                                                                                      cgGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCTTCATTCAGGCGTGGCAGACTGTTAtt            <  1:1612247/150‑1 (MQ=255)
                                                                                                                                          aCGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCTTCATTCAGGCGTGGCAGACTGTTATTGAtt        <  2:2075831/150‑1 (MQ=255)
                                                                                                                                            gCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCTTCATTCAGGCGTGGCAGACTGTTATTGATTaa      <  2:4119282/150‑1 (MQ=255)
                                                                                                                                            gCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCTTCATTCAGGCGTGGCAGACTGTTATTGATTaa      >  2:2901820/1‑150 (MQ=255)
                                                                                                                                            gCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCTTCATTCAGGCGTGGCAGACTGTTATTGATTaa      >  2:3048253/1‑150 (MQ=255)
                                                                                                                                               cAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCTTCATTCAGGCGTGGCAGACTGTTATTGATTAAAtt   <  1:4634583/150‑1 (MQ=255)
                                                                                                                                                aGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCTTCATTCAGGCGTGGCAGACTGTTATTGATTAAATTa  <  2:775946/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AAAAAGATCATCGCCCGCGGCAGATGTGATGCGGAGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCGATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCTTCATTCAGGCGTGGCAGACTGTTATTGATTAAATTA  >  CP009273/968268‑968561

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: