Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,405,946 T→G 94.4% Q43P (CAG→CCG)  yfbU ← UPF0304 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,405,9460TG94.4% 962.3 / 24.8 248Q43P (CAG→CCG) yfbUUPF0304 family protein
Reads supporting (aligned to +/- strand):  ref base T (8/6);  new base G (112/122);  total (120/128)
Fisher's exact test for biased strand distribution p-value = 5.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

Only 100 of 248 total aligned reads displayed.

GGTGACACGACGTTCATCGATGGATTGCTGGTCCTGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCTGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCAC  >  CP009273/2405801‑2406091
                                                                                                                                                 |                                                                                                                                                 
  tGACACGACGTTCATCGATGGATTGCTGGTCCTGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATc                                                                                                                                             <  2:2045170/150‑1 (MQ=255)
    acacGACGTTCATCGATGGATTGCTGGTCCTGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCg                                                                                                                                           <  2:974494/150‑1 (MQ=255)
     cacGACGTTCATCGATGGATTGCTGGTCCTGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGt                                                                                                                                          <  2:3990026/150‑1 (MQ=255)
       cgacgTTCATCGATGGATTGCTGGTCCTGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTaa                                                                                                                                        >  1:1169141/1‑150 (MQ=255)
        gacgTTCATCGATGGATTGCTGGTCCTGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAAc                                                                                                                                       >  2:2033534/1‑150 (MQ=255)
           gTTCATCTATGGTTTGCTGGTACTGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCAc                                                                                                                                    <  2:1865958/150‑1 (MQ=255)
           gTTCATCGATGGAGTGCTGGTCCTGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCAc                                                                                                                                    <  2:4270636/150‑1 (MQ=255)
              cATCGATGGATTGCTGGTCCTGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCt                                                                                                                                 <  1:383376/150‑1 (MQ=255)
                  gATGGATTGCTGGTCCTGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCaat                                                                                                                             <  1:3612969/150‑1 (MQ=255)
                  gATGGATTGCTGGTCCTGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCaat                                                                                                                             >  2:1370919/1‑150 (MQ=255)
                      gATTGCTGGTCCTGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAAtt                                                                                                                         <  2:2453357/150‑1 (MQ=255)
                      gATTGCTGGTCCTGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAAtt                                                                                                                         >  1:158829/1‑150 (MQ=255)
                      gATTGCTGGTCCTGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAAtt                                                                                                                         <  2:1758867/150‑1 (MQ=255)
                       aTTGCTGGTCCTGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTg                                                                                                                        <  1:147222/150‑1 (MQ=255)
                        ttGCTGGTCCTGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGt                                                                                                                       <  2:2820509/150‑1 (MQ=255)
                        ttGCTGGTCCTGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGt                                                                                                                       <  1:1327400/150‑1 (MQ=255)
                        ttGCTGGTCCTGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGt                                                                                                                       <  1:3738539/150‑1 (MQ=255)
                           cTGGTCCTGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGAGGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTg                                                                                                                    <  2:1361398/150‑1 (MQ=255)
                            tGGTCCTGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGc                                                                                                                   >  2:3414530/1‑150 (MQ=255)
                              gTCCTGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAg                                                                                                                 <  1:813298/150‑1 (MQ=255)
                                 cTGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTTTTTGCAGGCg                                                                                                              <  1:503095/150‑1 (MQ=255)
                                  tGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGa                                                                                                             <  2:3799129/150‑1 (MQ=255)
                                     aaaTTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACgg                                                                                                          >  2:222444/1‑150 (MQ=255)
                                       aTTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTa                                                                                                        <  1:4359823/150‑1 (MQ=255)
                                        ttAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTaa                                                                                                       <  2:4472826/150‑1 (MQ=255)
                                        ttAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTaa                                                                                                       <  2:1156976/150‑1 (MQ=255)
                                         tAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAAc                                                                                                      <  1:1541608/150‑1 (MQ=255)
                                         tAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAAc                                                                                                      <  2:843183/150‑1 (MQ=255)
                                          aGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACg                                                                                                     >  1:4413032/1‑150 (MQ=255)
                                              cAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCTGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCa                                                                                                 <  1:236346/150‑1 (MQ=255)
                                               aGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAg                                                                                                <  2:77818/150‑1 (MQ=255)
                                                ggAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGc                                                                                               <  1:3021525/150‑1 (MQ=255)
                                                ggAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGc                                                                                               <  1:845973/150‑1 (MQ=255)
                                                ggAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGc                                                                                               >  2:982331/1‑150 (MQ=255)
                                                 gAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCa                                                                                              <  2:413440/150‑1 (MQ=255)
                                                 gAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCa                                                                                              >  2:594642/1‑150 (MQ=255)
                                                  aaaCATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCAt                                                                                             >  1:2747516/1‑150 (MQ=255)
                                                  aaaCATGCAACACATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCAt                                                                                             >  2:4081483/1‑150 (MQ=255)
                                                   aaCATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCAtt                                                                                            <  2:4438167/150‑1 (MQ=255)
                                                     cATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTg                                                                                          >  2:2382203/1‑150 (MQ=255)
                                                       tGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGcc                                                                                        >  2:1251772/1‑150 (MQ=255)
                                                         cAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCATGCGACGGTAACGTTAAGCATTTGCCgg                                                                                      >  2:3893792/1‑150 (MQ=255)
                                                         cAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTAACGCATCGGTAATACGTAACCACGCTCAATAATTGTTTGCATGCGACGGTAACGTTCAGCATTTGACgg                                                                                      >  2:3893979/1‑150 (MQ=255)
                                                           aCGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGAt                                                                                    <  1:1122870/150‑1 (MQ=255)
                                                            cGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATc                                                                                   <  2:1868250/150‑1 (MQ=255)
                                                            cGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATc                                                                                   >  1:3376346/1‑150 (MQ=255)
                                                                tGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGc                                                                               >  2:3614294/1‑150 (MQ=255)
                                                                 gATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCa                                                                              >  2:2354465/1‑150 (MQ=255)
                                                                 gATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCa                                                                              <  2:285968/150‑1 (MQ=255)
                                                                 gATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCa                                                                              >  2:4755106/1‑150 (MQ=255)
                                                                       tctcCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGtca                                                                        >  1:1646407/1‑150 (MQ=255)
                                                                        ctcCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGtcat                                                                       <  2:3529410/150‑1 (MQ=255)
                                                                        ctcCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGtcat                                                                       >  1:4530252/1‑150 (MQ=255)
                                                                         tcCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGtcatc                                                                      >  2:3063148/1‑150 (MQ=255)
                                                                          ccATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGtcatca                                                                     <  2:3830708/150‑1 (MQ=255)
                                                                            aTAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCatcatc                                                                   <  2:4564462/150‑1 (MQ=255)
                                                                              aaTGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCtt                                                                 >  2:2989377/1‑150 (MQ=255)
                                                                              aaTGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCtt                                                                 <  1:4496908/150‑1 (MQ=255)
                                                                               aTGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTg                                                                <  1:1441389/150‑1 (MQ=255)
                                                                                 gTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTa                                                              <  2:4690948/150‑1 (MQ=255)
                                                                                 gTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTa                                                              <  1:196215/150‑1 (MQ=255)
                                                                                  tCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCTGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTAc                                                             <  2:2797503/150‑1 (MQ=255)
                                                                                   cGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACt                                                            >  1:2825239/1‑150 (MQ=255)
                                                                                   cGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACt                                                            >  1:4113209/1‑150 (MQ=255)
                                                                                      tgatAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCTGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGAt                                                         >  1:3195544/1‑150 (MQ=255)
                                                                                          aGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGa                                                     >  1:3112570/1‑150 (MQ=255)
                                                                                          aGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGa                                                     >  1:1680451/1‑150 (MQ=255)
                                                                                           gTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGAc                                                    >  1:4279167/1‑150 (MQ=255)
                                                                                             gCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACaa                                                  <  1:2417871/150‑1 (MQ=255)
                                                                                             gCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACaa                                                  <  1:22442/150‑1 (MQ=255)
                                                                                             gCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACaa                                                  <  1:2083531/150‑1 (MQ=255)
                                                                                              cGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACaaa                                                 <  1:4655948/150‑1 (MQ=255)
                                                                                                gCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAAt                                               <  2:3145517/150‑1 (MQ=255)
                                                                                                  aGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATaa                                             >  2:3282901/1‑150 (MQ=255)
                                                                                                   ggTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAg                                            >  1:2941596/1‑150 (MQ=255)
                                                                                                        cttcttTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGtt                                       >  2:2432650/1‑150 (MQ=255)
                                                                                                        cttcttTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGtt                                       <  1:1760045/150‑1 (MQ=255)
                                                                                                         ttcttTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTg                                      <  1:2195432/150‑1 (MQ=255)
                                                                                                         ttcttTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTg                                      >  1:65186/1‑150 (MQ=255)
                                                                                                          tcttTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTgc                                     >  2:2080851/1‑150 (MQ=255)
                                                                                                              tCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGt                                 >  1:3192976/1‑150 (MQ=255)
                                                                                                              tCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGt                                 >  1:3979146/1‑150 (MQ=255)
                                                                                                                aGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTg                               <  1:2033534/150‑1 (MQ=255)
                                                                                                                  cTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGt                             <  1:3975865/150‑1 (MQ=255)
                                                                                                                    cGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCa                           <  2:878523/150‑1 (MQ=255)
                                                                                                                     gCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCAt                          <  1:2952337/150‑1 (MQ=255)
                                                                                                                      ccAAACTCGCGATCCAATTCACGCATCTGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCAtt                         >  2:658633/1‑150 (MQ=255)
                                                                                                                      ccAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCAtt                         >  1:3939621/1‑150 (MQ=255)
                                                                                                                       cAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCAttt                        >  2:4128307/1‑150 (MQ=255)
                                                                                                                       cAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCAttt                        >  2:1152782/1‑150 (MQ=255)
                                                                                                                        aaaCTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTc                       >  2:470246/1‑150 (MQ=255)
                                                                                                                           cTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCAt                    >  2:2749115/1‑150 (MQ=255)
                                                                                                                             cgcgATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCAttg                  <  2:1330289/150‑1 (MQ=255)
                                                                                                                             cgcgATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCAttg                  <  1:2605027/150‑1 (MQ=255)
                                                                                                                             cgcgATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCAttg                  <  1:2606333/150‑1 (MQ=255)
                                                                                                                                gATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCAttgttg               >  2:2480467/1‑150 (MQ=255)
                                                                                                                                       ttCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTg        <  2:175318/150‑1 (MQ=255)
                                                                                                                                       ttCACGCATCGGTAAGCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTg        <  2:2579023/150‑1 (MQ=255)
                                                                                                                                         cACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTa      <  2:1432511/150‑1 (MQ=255)
                                                                                                                                             cATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCAc  >  1:4708730/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
GGTGACACGACGTTCATCGATGGATTGCTGGTCCTGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCTGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCAC  >  CP009273/2405801‑2406091

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: