Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,405,946 T→G 100% 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,9460TG96.9% 1010.9 / 8.7 257Q43P (CAG→CCG) yfbUUPF0304 family protein
Reads supporting (aligned to +/- strand):  ref base T (5/3);  new base G (112/137);  total (117/140)
Fisher's exact test for biased strand distribution p-value = 4.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.14e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 257 total aligned reads displayed.

GGTGACACGACGTTCATCGATGGATTGCTGGTCCTGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCTGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACT  >  CP009273/2405801‑2406092
                                                                                                                                                 |                                                                                                                                                  
ggTGACACGACGTTCATCGATGGATTGCTGGTCCTGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTaa                                                                                                                                                >  1:1417213/1‑150 (MQ=255)
  tGACACGACGTTCATCGATGGATTGCTGGTCCTGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATc                                                                                                                                              <  2:3720/150‑1 (MQ=255)
  tGACACGACGTTCATCGATGGATTGCTGGTCCTGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATc                                                                                                                                              >  2:408368/1‑150 (MQ=255)
   gACACGACGTTCATCGATGGATTGCTGGTCCTGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATcc                                                                                                                                             <  2:1785955/150‑1 (MQ=255)
   gACACGACGTTCATCGATGGATTGCTGGTCCTGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATcc                                                                                                                                             <  2:676584/150‑1 (MQ=255)
   gACACGACGTTCATCGATGGATTGCTGGTCCTGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATcc                                                                                                                                             <  1:81949/150‑1 (MQ=255)
      acgacgTTCATCGATGGATTGCTGGTCCTGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTa                                                                                                                                          <  1:1546630/150‑1 (MQ=255)
      acgacgTTCATCGATGGATTGCTGGTCCTGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTa                                                                                                                                          <  2:1728077/150‑1 (MQ=255)
        gacgTTCATCGATGGATTGCTGGTCCTGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAAc                                                                                                                                        >  1:2375709/1‑150 (MQ=255)
        gacgTTCATCGATGGATTGCTGGTCCTGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAAc                                                                                                                                        <  1:3940451/150‑1 (MQ=255)
           gTTCATCGATGGATTGCTGGTCCTGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCAc                                                                                                                                     >  2:4220002/1‑150 (MQ=255)
            ttCATCGATGGATTGCTGGTCCTGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACg                                                                                                                                    <  1:1178695/150‑1 (MQ=255)
            ttCATCGATGGATTGCTGGTCCTGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACg                                                                                                                                    <  2:2485244/150‑1 (MQ=255)
               aTCGATGGATTGCTGGTCCTGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTc                                                                                                                                 >  1:3443102/1‑150 (MQ=255)
                tCGATGGATTGCTGGTCCTGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCa                                                                                                                                <  1:3628015/150‑1 (MQ=255)
                tCGATGGATTGCTGGTCCTGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCa                                                                                                                                >  1:2518699/1‑150 (MQ=255)
                   aTGGATTGCTGGTCCTGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCaata                                                                                                                             >  2:425164/1‑150 (MQ=255)
                   aTGGATTGCTGGTCCTGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCaata                                                                                                                             >  1:3773445/1‑150 (MQ=255)
                   aTGGATTGCTGGTCCTGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCaata                                                                                                                             >  1:4426146/1‑150 (MQ=255)
                    tGGATTGCTGGTCCTGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCaataa                                                                                                                            >  2:149812/1‑150 (MQ=255)
                    tGGATTGCTGGTCCTGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCaataa                                                                                                                            >  1:85302/1‑150 (MQ=255)
                       aTTGCTGGTCCTGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTg                                                                                                                         >  2:3323436/1‑150 (MQ=255)
                        ttGCTGGTCCTGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGt                                                                                                                        <  2:1596481/150‑1 (MQ=255)
                          gCTGGTCCTGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGttt                                                                                                                      <  2:1382551/150‑1 (MQ=255)
                            tGGTCCTGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTGTGc                                                                                                                    >  1:1002677/1‑150 (MQ=255)
                            tGGTCCTGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACACGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGc                                                                                                                    >  1:1757284/1‑150 (MQ=255)
                              gTCCTGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAg                                                                                                                  <  1:1022998/150‑1 (MQ=255)
                                     aaaTTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCTGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACgg                                                                                                           >  2:3710388/1‑150 (MQ=255)
                                     aaaTTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGAGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACgg                                                                                                           >  2:4531367/1‑150 (MQ=255)
                                      aaTTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGt                                                                                                          >  1:2582062/1‑150 (MQ=255)
                                         tAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAAc                                                                                                       <  1:1493301/150‑1 (MQ=255)
                                         tAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAAc                                                                                                       <  1:2237872/150‑1 (MQ=255)
                                         tAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAAc                                                                                                       >  2:181453/1‑150 (MQ=255)
                                          aGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACg                                                                                                      <  2:1135263/150‑1 (MQ=255)
                                          aGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACg                                                                                                      <  1:3212669/150‑1 (MQ=255)
                                          aGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACg                                                                                                      <  1:894975/150‑1 (MQ=255)
                                          aGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACg                                                                                                      <  1:1393138/150‑1 (MQ=255)
                                          aGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACg                                                                                                      >  1:468287/1‑150 (MQ=255)
                                              cAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCa                                                                                                  >  2:213763/1‑150 (MQ=255)
                                               aGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAg                                                                                                 >  2:809922/1‑150 (MQ=255)
                                               aGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAg                                                                                                 >  2:3127306/1‑150 (MQ=255)
                                                 gAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCa                                                                                               <  2:3958613/150‑1 (MQ=255)
                                                  aaaCATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCAt                                                                                              <  1:955930/150‑1 (MQ=255)
                                                    aCATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCAttt                                                                                            <  2:4473578/150‑1 (MQ=255)
                                                     cATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTg                                                                                           <  2:2734247/150‑1 (MQ=255)
                                                     cATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTg                                                                                           <  2:1673981/150‑1 (MQ=255)
                                                        gCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCTGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCg                                                                                        <  1:2036354/150‑1 (MQ=255)
                                                         cAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCTGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCgg                                                                                       <  1:824122/150‑1 (MQ=255)
                                                           aCGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGAt                                                                                     >  1:2310886/1‑150 (MQ=255)
                                                               aTGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAg                                                                                 <  2:2934759/150‑1 (MQ=255)
                                                                      aTCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGtc                                                                          >  1:2739579/1‑150 (MQ=255)
                                                                        ctcCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGtcat                                                                        >  2:2253144/1‑150 (MQ=255)
                                                                         tcCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGtcatc                                                                       <  1:3310744/150‑1 (MQ=255)
                                                                         tcCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGtcatc                                                                       >  2:4551876/1‑150 (MQ=255)
                                                                             tAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCt                                                                   <  2:2572170/150‑1 (MQ=255)
                                                                              aaTGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCtt                                                                  <  1:2699134/150‑1 (MQ=255)
                                                                              aaTGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCtt                                                                  <  2:836123/150‑1 (MQ=255)
                                                                              aaTGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTAATCATCtt                                                                  <  2:1793952/150‑1 (MQ=255)
                                                                               aTGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTg                                                                 >  1:1692599/1‑150 (MQ=255)
                                                                                tGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGt                                                                <  1:370843/150‑1 (MQ=255)
                                                                                  tCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTAc                                                              >  1:1169294/1‑150 (MQ=255)
                                                                                   cGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACt                                                             <  1:3953594/150‑1 (MQ=255)
                                                                                    gatgatAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTg                                                            <  1:986162/150‑1 (MQ=255)
                                                                                    gatgatAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTg                                                            <  2:467388/150‑1 (MQ=255)
                                                                                    gatgatAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTg                                                            <  2:630795/150‑1 (MQ=255)
                                                                                      tgatAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGAt                                                          >  2:2121959/1‑150 (MQ=255)
                                                                                           gTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGAc                                                     <  2:1032498/150‑1 (MQ=255)
                                                                                              cGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCTGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACaaa                                                  >  2:4429105/1‑150 (MQ=255)
                                                                                               ggCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACaaaa                                                 <  1:2894051/150‑1 (MQ=255)
                                                                                                  aGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATaa                                              <  1:3881623/150‑1 (MQ=255)
                                                                                                   ggTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAg                                             >  2:1926713/1‑150 (MQ=255)
                                                                                                   ggTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAg                                             <  2:3629515/150‑1 (MQ=255)
                                                                                                   ggTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAg                                             <  2:3153392/150‑1 (MQ=255)
                                                                                                    gTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGa                                            <  2:3177244/150‑1 (MQ=255)
                                                                                                     tttCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGAc                                           <  2:4088926/150‑1 (MQ=255)
                                                                                                         ttcttTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTg                                       >  2:1188265/1‑150 (MQ=255)
                                                                                                               cAGCTCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGtt                                 <  2:2963349/150‑1 (MQ=255)
                                                                                                                   tCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTc                             <  1:248553/150‑1 (MQ=255)
                                                                                                                   tCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTc                             >  1:3017876/1‑150 (MQ=255)
                                                                                                                   tCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTAAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGAAGTTGCGCGTTGGTc                             >  1:995901/1‑150 (MQ=255)
                                                                                                                   tCGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAAACACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCTTTGGTc                             >  1:3274296/1‑150 (MQ=255)
                                                                                                                    cGCCAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCa                            >  2:3858148/1‑150 (MQ=255)
                                                                                                                      ccAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCAtt                          <  2:2678118/150‑1 (MQ=255)
                                                                                                                       cAAACTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCAttt                         <  2:871103/150‑1 (MQ=255)
                                                                                                                         aaCTCGCGATCCAATTCACGCATCGGTACTCCGTAACCACGCTCAATAATTGTTTGCAGTCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGAAGTTGCGCGTTGGTCATTTcc                       <  2:3402100/150‑1 (MQ=255)
                                                                                                                           cTCGCGATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCAt                     >  1:1403568/1‑150 (MQ=255)
                                                                                                                           cTCGCGATCCAATTCACGCATCGCTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCAt                     >  1:2066603/1‑150 (MQ=255)
                                                                                                                               cgATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCAttgtt                 <  2:2501198/150‑1 (MQ=255)
                                                                                                                                gATCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCAttgttg                >  1:2730088/1‑150 (MQ=255)
                                                                                                                                 aTCCAATTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGa               >  1:4298663/1‑150 (MQ=255)
                                                                                                                                    aaaTTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAAAGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTc            >  1:3487537/2‑150 (MQ=255)
                                                                                                                                      aTTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCt          <  2:3164503/150‑1 (MQ=255)
                                                                                                                                      aTTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCt          <  2:1858830/150‑1 (MQ=255)
                                                                                                                                      aTTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCt          >  1:1873730/1‑150 (MQ=255)
                                                                                                                                      aTTCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCt          <  2:2229636/150‑1 (MQ=255)
                                                                                                                                       ttCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTg         <  1:3801060/150‑1 (MQ=255)
                                                                                                                                       ttCACGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTg         <  1:685996/150‑1 (MQ=255)
                                                                                                                                           cGCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCTTATTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATc     >  1:1989125/1‑150 (MQ=255)
                                                                                                                                            gCATCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATNa    <  2:403620/150‑1 (MQ=255)
                                                                                                                                              aTCGGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCAct  >  1:2627483/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
GGTGACACGACGTTCATCGATGGATTGCTGGTCCTGCAAATTAGACCAGGAAACATGCAACGCATGATACATCTCCATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCTGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACT  >  CP009273/2405801‑2406092

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