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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273968,4150GT85.0% 753.2 / 128.9 291S148* (TCG→TAG) ycbCDUF218 superfamily protein
Reads supporting (aligned to +/- strand):  ref base G (25/19);  new base T (123/124);  total (148/143)
Fisher's exact test for biased strand distribution p-value = 4.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 292 total aligned reads displayed.

GAAAAAAGATCATCGCCCGCGGCAGATGTGATGCGGAGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCGATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCTTCATTCAGGCGTGGCAGACTGTTATTGATTAAATTA  >  CP009273/968266‑968561
                                                                                                                                                     |                                                                                                                                                  
   aaaaGATCATCGCCCGCGGCAGATGTGATGCGGAGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTAtt                                                                                                                                                 >  1:3478097/1‑150 (MQ=255)
   aaaaGATCATCGCCCGCGGCAGATGTGATGCGGAGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTAtt                                                                                                                                                 >  1:3884840/1‑150 (MQ=255)
   aaaaGATCATCGCCCGCGGCAGATGTGATGCGGAGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTAtt                                                                                                                                                 >  1:1607596/1‑150 (MQ=255)
   aaaaGATCATCGCCCGCGGCAGATGTGATGCGGAGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCGAtt                                                                                                                                                 <  1:2465586/150‑1 (MQ=255)
     aaGATCATCGCCCGCGGCAGATGTGATGCGGAGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTCTTgc                                                                                                                                               >  1:1083726/1‑150 (MQ=255)
     aaGATCATCGCCCGCGGCAGATGTGATGCGGAGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTgc                                                                                                                                               <  1:2315460/150‑1 (MQ=255)
       gATCATCGCCCGCGGCAGATGTGATGCGGAGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCGATTgcgc                                                                                                                                             >  1:3918659/1‑150 (MQ=255)
        atcatcGCCCGCGGCAGATGTGATGCGGAGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCGATTGCGCa                                                                                                                                            >  1:2359887/1‑150 (MQ=255)
          catcGCCCGCGGCAGATGTGATGCGGAGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAAc                                                                                                                                          >  1:693472/1‑150 (MQ=255)
          catcGCCCGCGGCAGATGTGATGCGGAGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAAc                                                                                                                                          >  1:2174328/1‑150 (MQ=255)
           atcGCCCGCGGCAGATGTGATGCGGAGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAAct                                                                                                                                         <  1:4343210/150‑1 (MQ=255)
             cGCCCGCGGCAGATGTGATGCGGAGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAActct                                                                                                                                       >  2:2927788/1‑150 (MQ=255)
               cccGCGGCAGATGTGATGCGGAGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTCTTGCGCAACTCTgg                                                                                                                                     >  1:3845034/1‑150 (MQ=255)
                   cgGCAGATGTGATGCGGAGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCAcc                                                                                                                                 <  1:2985411/150‑1 (MQ=255)
                    ggCAGATGTGATGCGGAGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTCTTGCGCAACTCTGGCACCt                                                                                                                                >  1:2734962/1‑150 (MQ=255)
                      cAGATGTGATGCGGAGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTAc                                                                                                                              >  2:1692106/1‑150 (MQ=255)
                       aGATGTGATGCGGAGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACt                                                                                                                             >  2:269502/1‑150 (MQ=255)
                       aGATGTGATGCGGAGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACt                                                                                                                             <  1:793015/150‑1 (MQ=255)
                          tgtgATGCGGAGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTcc                                                                                                                          >  1:3264008/1‑150 (MQ=255)
                          tgtgATGCGGAGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTcc                                                                                                                          >  2:3060245/1‑150 (MQ=255)
                             gATGCGGAGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCt                                                                                                                       <  1:4427179/150‑1 (MQ=255)
                             gATGCGGAGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCGTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCt                                                                                                                       <  1:1871909/150‑1 (MQ=255)
                               tGCGGAGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCGATTGCGCAACTCTGGCACCTACTTCCGCTGt                                                                                                                     <  1:1504794/150‑1 (MQ=255)
                                 cGGAGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTAc                                                                                                                   <  2:4187022/150‑1 (MQ=255)
                                   gAGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTc                                                                                                                 >  2:1919762/1‑150 (MQ=255)
                                    aGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCa                                                                                                                <  1:4062477/150‑1 (MQ=255)
                                    aGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCa                                                                                                                <  2:27553/150‑1 (MQ=255)
                                    aGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCGATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCa                                                                                                                <  2:787366/150‑1 (MQ=255)
                                     ggTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCAc                                                                                                               >  1:3143735/1‑150 (MQ=255)
                                      gTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCAcc                                                                                                              <  2:4405997/150‑1 (MQ=255)
                                         aCCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCGATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTa                                                                                                           <  1:2612443/150‑1 (MQ=255)
                                           cagcagGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTAtt                                                                                                         <  2:2356553/150‑1 (MQ=255)
                                            agcagGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTg                                                                                                        <  1:856850/150‑1 (MQ=255)
                                            agcagGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTg                                                                                                        <  1:2910473/150‑1 (MQ=255)
                                            agcagGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTg                                                                                                        <  2:3179790/150‑1 (MQ=255)
                                            agcagGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCGATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTg                                                                                                        >  1:3176350/1‑150 (MQ=255)
                                              cagGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGt                                                                                                      >  2:1620378/1‑150 (MQ=255)
                                                ggAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGttt                                                                                                    <  1:949833/150‑1 (MQ=255)
                                                  aaaGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCGATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGttttt                                                                                                  >  1:3718678/1‑150 (MQ=255)
                                                         cGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGc                                                                                           >  2:3770336/1‑150 (MQ=255)
                                                          gTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGcc                                                                                          >  1:1355714/1‑150 (MQ=255)
                                                             aCCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTcc                                                                                       <  2:2581016/150‑1 (MQ=255)
                                                               cAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCGATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCg                                                                                     >  1:3636387/1‑150 (MQ=255)
                                                               cAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCGATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCg                                                                                     >  1:1352662/1‑150 (MQ=255)
                                                                  tCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGa                                                                                  >  2:3163063/1‑150 (MQ=255)
                                                                  tCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGa                                                                                  <  2:623620/150‑1 (MQ=255)
                                                                     ccTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGa                                                                               >  1:2374007/1‑150 (MQ=255)
                                                                     ccTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGa                                                                               >  1:3330694/1‑150 (MQ=255)
                                                                      cTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCGATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGag                                                                              >  1:805083/1‑149 (MQ=255)
                                                                      cTGTTTCACTGCTGCAGCTTCGTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCGATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGAt                                                                              <  2:3795923/150‑1 (MQ=255)
                                                                        gTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCATATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCa                                                                            <  2:1076042/150‑1 (MQ=255)
                                                                        gTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCa                                                                            <  2:3362674/150‑1 (MQ=255)
                                                                        gTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCa                                                                            <  2:120969/150‑1 (MQ=255)
                                                                        gTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCa                                                                            <  2:4460852/150‑1 (MQ=255)
                                                                           tCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGtt                                                                         <  1:2276441/150‑1 (MQ=255)
                                                                            cACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGttt                                                                        <  1:2199708/150‑1 (MQ=255)
                                                                            cACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGttt                                                                        >  1:3750460/1‑150 (MQ=255)
                                                                             aCTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGtttt                                                                       <  2:2144130/150‑1 (MQ=255)
                                                                             aCTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCGATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGtttt                                                                       >  1:3520706/1‑150 (MQ=255)
                                                                             aCTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCGATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGtttt                                                                       >  1:1235237/1‑150 (MQ=255)
                                                                                   gcAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATccc                                                                 <  2:4129510/150‑1 (MQ=255)
                                                                                     aGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCgg                                                               >  1:377461/1‑150 (MQ=255)
                                                                                      gCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGa                                                              <  2:3478097/150‑1 (MQ=255)
                                                                                        ttcttcttcGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCGATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGAtt                                                            >  1:3086252/1‑150 (MQ=255)
                                                                                           ttcttcGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTc                                                         >  2:4055190/1‑150 (MQ=255)
                                                                                           ttcttcGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTc                                                         <  1:1788495/150‑1 (MQ=255)
                                                                                           ttcttcGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTc                                                         <  1:3994058/150‑1 (MQ=255)
                                                                                           ttcttcGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTc                                                         <  1:4115506/150‑1 (MQ=255)
                                                                                                  gTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCAt                                                  >  2:864104/1‑150 (MQ=255)
                                                                                                   tATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATa                                                 >  2:2346222/1‑150 (MQ=255)
                                                                                                   tATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATa                                                 >  2:2876048/1‑150 (MQ=255)
                                                                                                   tATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATa                                                 >  1:635187/1‑150 (MQ=255)
                                                                                                          tGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAAt                                          <  2:1508751/150‑1 (MQ=255)
                                                                                                           ggCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATa                                         >  2:1953651/1‑150 (MQ=255)
                                                                                                              aaaTCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCt                                      <  2:2104699/150‑1 (MQ=255)
                                                                                                              aaaTCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAAGACCt                                      >  2:3269004/1‑150 (MQ=255)
                                                                                                              aaaTCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCAACTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTACCGTGAAGATCAGTTTTGAGCCCGGATTTTCGCGCCATAAGCGAATACCt                                      <  2:3041404/150‑1 (MQ=255)
                                                                                                                aTCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCttc                                    <  2:3394045/150‑1 (MQ=255)
                                                                                                                  ccAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCttcat                                  <  1:1045959/150‑1 (MQ=255)
                                                                                                                  ccAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCttcat                                  >  2:3042060/1‑150 (MQ=255)
                                                                                                                   cAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCttcatt                                 >  1:2215906/1‑150 (MQ=255)
                                                                                                                    aGGGTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCttcattc                                <  1:2015532/150‑1 (MQ=255)
                                                                                                                     gggTGATAATTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCttcattca                               <  1:2767232/150‑1 (MQ=255)
                                                                                                                     gggTGATAATTTGCTCGCGCGGCACGCCCAGCGATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCttcattca                               <  1:1668894/150‑1 (MQ=255)
                                                                                                                            aaTTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACGTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCTTCATTCAGGCGTgg                        <  2:3432917/150‑1 (MQ=255)
                                                                                                                             aTTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCTTCATTCAGGCGTGGc                       >  2:1726550/1‑150 (MQ=255)
                                                                                                                             aTTTGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCTTCATTCAGGCGTGGc                       >  2:1883707/1‑150 (MQ=255)
                                                                                                                               ttGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCTTCATTCAGGCGTGGCAg                     <  1:3393764/150‑1 (MQ=255)
                                                                                                                               ttGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCTTCATTCAGGCGTGGCAg                     <  2:217549/150‑1 (MQ=255)
                                                                                                                                tGCTCGCGCGGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCTTCATTCAGGCGTGGCAGa                    <  2:298426/150‑1 (MQ=255)
                                                                                                                                    cgcgcgGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCTTCATTCAGGCGTGGCAGACTGt                <  2:478356/150‑1 (MQ=255)
                                                                                                                                      cgcgGCACGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCTTCATTCAGGCGTGGCAGACTGttc              >  1:1741097/1‑149 (MQ=255)
                                                                                                                                            aCGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCTTCATTCAGGCGTGGCAGACTGTTATTGAtt        <  1:3326938/150‑1 (MQ=255)
                                                                                                                                             cGCCCAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCTTCATTCAGGCGTGGCAGACTGTTATTGATTa       <  2:3691439/150‑1 (MQ=255)
                                                                                                                                              gCCCAGCGATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCTTCATTCAGGCGTGGCAGACTGTTATTGATTaa      <  2:4379407/150‑1 (MQ=255)
                                                                                                                                                ccAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCTTCATTCAGGCGTGGCAGACTGTTATTGATTAAAt    <  2:1287462/150‑1 (MQ=255)
                                                                                                                                                 cAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCTTCATTCAGGCGTGGCAGACTGTTATTGATTAAAtt   <  2:3568458/150‑1 (MQ=255)
                                                                                                                                                 cAGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCTTCATTCAGGCGTGGCAGACTGTTATTGATTAAAtt   <  1:4534212/150‑1 (MQ=255)
                                                                                                                                                  aGCTATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCTTCATTCAGGCGTGGCAGACTGTTATTGATTAAATTa  >  1:3562734/1‑150 (MQ=255)
                                                                                                                                                  aGCGATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCTTCATTCAGGCGTGGCAGACTGTTATTGATTAAATTa  <  2:3918659/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GAAAAAAGATCATCGCCCGCGGCAGATGTGATGCGGAGGTCACCAGCAGGAAAGGGGCGTCACCAATCGCCTGTTTCACTGCTGCAGCTTCTTCTTCGGTATCTTTTGGCAAATCCAGGGTGATAATTTGCTCGCGCGGCACGCCCAGCGATTGCGCAACTCTGGCACCTACTTCCGCTGTACTCACCGTATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCTTCATTCAGGCGTGGCAGACTGTTATTGATTAAATTA  >  CP009273/968266‑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: