Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,246,190 G→A 22.9% V35V (GTC→GTT dhaK ← dihydroxyacetone kinase, PTS‑dependent, dihydroxyacetone‑binding subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,246,1900GA22.9% 559.6 / 165.4 257V35V (GTC→GTTdhaKdihydroxyacetone kinase, PTS‑dependent, dihydroxyacetone‑binding subunit
Reads supporting (aligned to +/- strand):  ref base G (98/100);  new base A (25/34);  total (123/134)
Fisher's exact test for biased strand distribution p-value = 3.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 261 total aligned reads displayed.

TTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATGCAGTGTCAGCGATGGATGCGCTTTCGCCAGTCCTGCCAGTTGTTCGTCCAGTACGTCTTGCACATCATTGATCAATTTTTTCATTATTTTGCTCCAGCAATTACGGTAGGGCATGGATGATGTTCAAC  >  CP009273/1246041‑1246338
                                                                                                                                                     |                                                                                                                                                    
ttATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTg                                                                                                                                                      <  1:1553732/150‑1 (MQ=255)
ttATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTg                                                                                                                                                      <  2:3102541/150‑1 (MQ=255)
   tCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACa                                                                                                                                                   <  2:775368/150‑1 (MQ=255)
    cGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACAt                                                                                                                                                  >  1:1450881/1‑150 (MQ=255)
    cGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACAt                                                                                                                                                  >  2:3322710/1‑150 (MQ=255)
      ggCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATac                                                                                                                                                >  2:1439914/1‑150 (MQ=255)
         gTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACAcc                                                                                                                                             >  1:1461394/1‑150 (MQ=255)
               gaggtgaAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGCCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATcc                                                                                                                                       <  1:815052/150‑1 (MQ=255)
                 ggtgaAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTg                                                                                                                                     >  1:1662419/1‑150 (MQ=255)
                 ggtgaAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTg                                                                                                                                     >  1:2858300/1‑150 (MQ=255)
                 ggtgaAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTg                                                                                                                                     >  2:4675175/1‑150 (MQ=255)
                    gaAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATg                                                                                                                                  >  1:1076260/1‑150 (MQ=255)
                            cGCCCGGACAGGCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATGCAGTGTCa                                                                                                                          >  2:4634437/1‑150 (MQ=255)
                            cGCCCGGACAGGCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTAACATACACCGGATCCTGATGCAGTGTCa                                                                                                                          >  2:1899070/1‑150 (MQ=255)
                                cGGACAGGCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTAACATACACCGGATCCTGATGCAGTGTCAGCga                                                                                                                      <  2:1338972/150‑1 (MQ=255)
                                 ggACAGGCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATGCAGTGTCAGCgat                                                                                                                     >  2:1389153/1‑150 (MQ=255)
                                     aGGCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTAACATACACCGGATCCTGATGCAGTGTCAGCgatggat                                                                                                                 <  1:4183564/150‑1 (MQ=255)
                                      ggCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATGCAGTGTCAGCgatggatg                                                                                                                >  2:4139366/1‑150 (MQ=255)
                                        cccccGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATGCAGTGTCAGCGATGGATgcg                                                                                                              >  2:4775252/1‑150 (MQ=255)
                                           ccGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATGCAGTGTCAGCGATGGATGCGCtt                                                                                                           >  1:1201278/1‑150 (MQ=255)
                                           ccGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATGCAGTGTCAGCGATGGATGCGCtt                                                                                                           <  1:4031133/150‑1 (MQ=255)
                                           ccGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATGCAGTGTCAGCGATGGATGCGCtt                                                                                                           <  1:3826039/150‑1 (MQ=255)
                                           ccGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATGCAGTGTCAGCGATGGATGCGCtt                                                                                                           >  1:3662710/1‑150 (MQ=255)
                                             gAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATGCAGTGTCAGCGATGGATGCGCTTTc                                                                                                         <  2:4595785/150‑1 (MQ=255)
                                                aGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATGCAGTGTCAGCGATGGATGCGCTTTCGcc                                                                                                      >  2:147818/1‑150 (MQ=255)
                                                 gCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATGCAGTGTCAGCGATGGATGCGCTTTCGCCa                                                                                                     >  2:3675960/1‑150 (MQ=255)
                                                   aTCCCCTGACCGATATAACAACAGTGCATCGGCTAGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCTGGTGAACTACACCGAATCATGATGAAAAGTCAGAGATGGATTATCTTACGCCAGt                                                                                                   >  1:3447901/1‑150 (MQ=255)
                                                    tCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTAACATACACCGGATCCTGATGCAGTGTCAGCGATGGATGCGCTTTCGCCAGTc                                                                                                  >  2:1288599/1‑150 (MQ=255)
                                                     ccccTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTAACATACACCGGATCCTGATGCAGTGTCAGCGATGGATGCGCTTTCGCCAGTcc                                                                                                 <  2:4543989/150‑1 (MQ=255)
                                                      cccTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATGCAGTGTCAGCGATGGATGCGCTTTCGCCAGTCCt                                                                                                <  1:407508/150‑1 (MQ=255)
                                                       ccTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATGCAGTGTCAGCGATGGATGCGCTTTCGCCAGTCCTg                                                                                               >  2:290658/1‑150 (MQ=255)
                                                        cTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATGCAGTGTCAGCGATGGATGCGCTTTCGCCAGTCCTGc                                                                                              >  1:2716068/1‑150 (MQ=255)
                                                          gACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTAACATACACCGGATCCTGATGCAGTGTCAGCGATGGATGCGCTTTCGCCAGTCCTGCCa                                                                                            <  1:1349974/150‑1 (MQ=255)
                                                          gACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTAACATACACCGGATCCTGATGCAGTGTCAGCGATGGATGCGCTTTCGCCAGTCCTGCCa                                                                                            <  2:1935524/150‑1 (MQ=255)
                                                            ccGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATGCAGTGTCAGCGATGGATGCGCTTTCGCCAGTCCTGCCAgt                                                                                          <  1:494142/150‑1 (MQ=255)
                                                             cGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATGCAGTGTCAGCGATGGATGCGCTTTCGCCAGTCCTGCCAgtt                                                                                         >  2:2435902/1‑150 (MQ=255)
                                                                tataACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATGCAGTGTCAGCGATGGATGCGCTTTCGCCAGTCCTGCCAgttgtt                                                                                      >  2:3017718/1‑150 (MQ=255)
                                                                tataACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATGCAGTGTCAGCGATGGATGCGCTTTCGCCAGTCCTGCCAgttgtt                                                                                      <  1:1234221/150‑1 (MQ=255)
                                                                 ataACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATGCAGTGTCAGCGATGGATGCGCTTTCGCCAGTCCTGCCAGTTGTTc                                                                                     <  1:3420105/150‑1 (MQ=255)
                                                                 ataACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATGCAGTGTCAGCGATGGATGCGCTTTCGCCAGTCCTGCCAGTTGTTc                                                                                     >  2:4488/1‑150 (MQ=255)
                                                                 ataACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGGGACAAACACCCGATCCTGATGAAGTGTAAGAGATGGATGCTCTGTCGCCAGTCCTGCCAGTTTTTc                                                                                     >  1:452974/1‑150 (MQ=255)
                                                                  taACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTAACATACACCGGATCCTGATGCAGTGTCAGCGATGGATGCGCTTTCGCCAGTCCTGCCAGTTGTTCg                                                                                    <  2:4050475/150‑1 (MQ=255)
                                                                      cacaGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTAACATACACCGGATCCTGATGCAGTGTCAGCGATGGATGCGCTTTCGCCAGTCCTGCCAGTTGTTCGTCCa                                                                                >  1:2572280/1‑150 (MQ=255)
                                                                       acaGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATGCAGTGTCAGCGATGGATGCGCTTTCGCCAGTCCTGCCAGTTGTTCGTCCAg                                                                               >  1:950723/1‑150 (MQ=255)
                                                                         aGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATGCAGTGTCAGCGATGGATGCGCTTTCGCCAGTCCTGCCAGTTGTTCGTCCGGTa                                                                             <  2:3947889/150‑1 (MQ=255)
                                                                         aGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATGCAGTGTCAGCGATGGATGCGCTTTCGCCAGTCCTGCCAGTTGTTCGTCCAGTa                                                                             >  1:1718038/1‑150 (MQ=255)
                                                                         aGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATGCAGTGTCAGCGATGGATGCGCTTTCGCCAGTCCTGCCAGTTGTTCGTCCAGTa                                                                             >  1:1843077/1‑150 (MQ=255)
                                                                               tCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATGCAGTGTCAGCGATGGATGCGCTTTCGCCAGTCCTGCCAGTTGTTCGTCCAGTACGTCtt                                                                       <  2:651553/150‑1 (MQ=255)
                                                                               tCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATGCAGTGTCAGCGATGGATGCGCTTTCGCCAGTCCTGCCAGTTGTTCGTCCAGTACGTCtt                                                                       <  1:4713318/150‑1 (MQ=255)
                                                                                cGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATGCAGTGTCAGCGATGGATGCGCTTTCGCCAGTCCTGCCAGTTGTTCGTCCAGTACGTCTTg                                                                      >  1:4655245/1‑150 (MQ=255)
                                                                                cGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATGCAGTGTCAGCGATGGATGCGCTTTCGCCAGTCCTGCCAGTTGTTCGTCCAGTACGTCTTg                                                                      <  1:4175211/150‑1 (MQ=255)
                                                                                 ggCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATGCAGTGTCAGCGATGGATGCGCTTTCGCCAGTCCTGCCAGTTGTTCGTCCAGTACGTCTTGc                                                                     >  1:4447616/1‑150 (MQ=255)
                                                                                 ggCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTAACATACACCGGATCCTGATGCAGTGTCAGCGATGGATGCGCTTTCGCCAGTCCTGCCAGTTGTTCGTCCAGTACGTCTTGc                                                                     >  2:2783756/1‑150 (MQ=255)
                                                                                      gtgtCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATGCAGTGTCAGCGATGGATGCGCTTTCGCCAGTCCTGCCAGTTGTTCGTCCAGTACGTCTTGCAcatc                                                                >  2:254627/1‑150 (MQ=255)
                                                                                      gtgtCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATGCAGTGTCAGCGATGGATGCGCTTTCGCCAGTCCTGCCAGTTGTTCGTCCAGTACGTCTTGCAcatc                                                                <  1:1000591/150‑1 (MQ=255)
                                                                                        gtCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATGCAGTGTCAGCGATGGATGCGCTTTCGCCAGTCCTGCCAGTTGTTCGTCCAGTACGTCTTGCAcatcat                                                              <  1:3276501/150‑1 (MQ=255)
                                                                                          ccGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATGCAGTGTCAGCGATGGATGCGCTTTCGCCAGTCCTGCCAGTTGTTCGTCCAGTACGTCTTGCACATCATTg                                                            <  1:4449814/150‑1 (MQ=255)
                                                                                          ccGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATGCAGTGTCAGCGATGGATGCGCTTTCGCCAGTCCTGCCAGTTGTTCGTCCAGTACGTCTTGCACATCATTg                                                            <  1:2844915/150‑1 (MQ=255)
                                                                                             cTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATGCAGTGTCAGCGATGGATGCGCTTTCGCCAGTCCTGCCAGTTGTTCGTCCAGTACGTCTTGCACATCATTGATc                                                         <  1:3637789/150‑1 (MQ=255)
                                                                                              tGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATGCAGTGTCAGCGATGGATGCGCTTTCGCCAGTCCTGCCAGTTGTTCGTCCAGTACGTCTTGCACATCATTGATCa                                                        >  2:2733417/1‑150 (MQ=255)
                                                                                              tGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATGCAGTGTCAGCGATGGATGCGCTTTCGCCAGTCCTGCCAGTTGTTCGTCCAGTACGTCTTGCACATCATTGATCa                                                        <  2:696351/150‑1 (MQ=255)
                                                                                              tGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATGCAGTGTCAGCGATGGATGCGCTTTCGCCAGTCCTGCCAGTTGTTCGTCCAGTACGTCTTGCACATCATTGATCa                                                        <  2:792537/150‑1 (MQ=255)
                                                                                                  gccACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATGCAGTGTCAGCGATGGATGCGCTTTCGCCAGTCCTGCCAGTTGTTCGTCCAGTACGTCTTGCACATCATTGATCAAttt                                                    >  2:1066634/1‑150 (MQ=255)
                                                                                                    cACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTAACATACACCGGATCCTGATGCAGTGTCAGCGATGGATGCGCTTTCGCCAGTCCTGCCAGTTGTTCGTCCAGTACGTCTTGCACATCATTGATCAAttttt                                                  <  2:3429919/150‑1 (MQ=255)
                                                                                                     aCCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATGCAGTGTCAGCGATGGATGCGCTTTCGCCAGTCCTGCCAGTTGTTCGTCCAGTACGTCTTGCACATCATTGATCAAtttttt                                                 <  1:2160353/150‑1 (MQ=255)
                                                                                                       ccGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATGCAGTGTCAGCGATGGATGCGCTTTCGCCAGTCCTGCCAGTTGTTCGTCCAGTACGTCTTGCACATCATTGATCAATTTTTTCa                                               >  1:2128303/1‑150 (MQ=255)
                                                                                                           cagcagGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTAACATACACCGGATCCTGATGCAGTGTCAGCGATGGATGCGCTTTCGCCAGTCCTGCCAGTTGTTCGTCCAGTACGTCTTGCACATCATTGATCAATTTTTTCattat                                           <  2:15823/150‑1 (MQ=255)
                                                                                                            agcagGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTAACATACACCGGATCCTGATGCAGTGTCAGCGATGGATGCGCTTTCGCCAGTCCTGCCAGTTGTTCGTCCAGTACGTCTTGCACATCATTGATCAATTTTTTCattatt                                          <  1:646674/150‑1 (MQ=255)
                                                                                                              cagGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATGCAGTGTCAGCGATGGATGCGCTTTCGCCAGTCCTGCCAGTTGTTCGTCCAGTACGTCTTGCACATCATTGATCAATTTTTTCATTAtttt                                        >  1:3708293/1‑150 (MQ=255)
                                                                                                              cagGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATGCAGTGTCAGCGATGGATGCGCTTTCGCCAGTCCTGCCAGTTGTTCGTCCAGTACGTCTTGCACATCATTGATCAATTTTTTCATTAtttt                                        <  1:3761667/150‑1 (MQ=255)
                                                                                                               agGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATGCAGTGTCAGCGATGGATGCGCTTTCGCCAGTCCTGCCAGTTGTTCGTCCAGTACGTCTTGCACATCATTGATCAATTTTTTCATTATTTTg                                       <  2:596730/150‑1 (MQ=255)
                                                                                                               agGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTAACATACACCGGATCCTGATGCAGTGTCAGCGATGGATGCGCTTTCGCCAGTCCTGCCAGTTGTTCGTCCAGTACGTCTTGCACATCATTGATCAATTTTTTCATTATTTTg                                       >  1:3938082/1‑150 (MQ=255)
                                                                                                                 ggCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATGCAGTGTCAGCGATGGATGCGCTTTCGCCAGTCCTGCCAGTTGTTCGTCCAGTACGTCTTGCACATCATTGATCAATTTTTTCATTATTTTGCt                                     >  2:3375909/1‑150 (MQ=255)
                                                                                                                 ggCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATGCAGTGTCAGCGATGGATGCGCTTTCGCCAGTCCTGCCAGTTGTTCGTCCAGTACGTCTTGCACATCATTGATCAATTTTTTCATTATTTTGCt                                     <  2:3221047/150‑1 (MQ=255)
                                                                                                                  gCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTAACATACACCGGATCCTGATGCAGTGTCAGCGATGGATGCGCTTTCGCCAGTCCTGCCAGTTGTTCGTCCAGTACGTCTTGCACATCATTGATCAATTTTTTCATTATTTTGCTc                                    <  1:1289285/150‑1 (MQ=255)
                                                                                                                    gACTTTTCCTGCAACAGGGGCATCAGCTCGGGTAACATACACCGGATCATGATGCAGTGTCAGCGATGGATGCGCTTTCGCCAGACCTGCCAGTTGTTCGTCCAGTACGACTTGCACATCATTGATCAATTTTTTCATTATTTAGCTCCa                                  >  2:2617734/1‑150 (MQ=255)
                                                                                                                       ttttCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATGCAGTGTCAGCGATGGATGCGCTTTCGCCAGTCCTGCCAGTTGTTCGTCCAGTACGTCTTGCACATCATTGATCAATTTTTTCATTATTTTGCTCCAGCa                               <  1:3044793/150‑1 (MQ=255)
                                                                                                                       ttttCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATGCAGTGTCAGCGATGGATGCGCTTTCGCCAGTCCTGCCAGTTGTTCGTCCAGTACGTCTTGCACATCATTGATCAATTTTTTCATTATTTTGCTCCAGCa                               >  2:2106475/1‑150 (MQ=255)
                                                                                                                        tttCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATGCAGTGTCAGCGATGGATGCGCTTTCGCCAGTCCTGCCAGTTGTTCGTCCAGTACGTCTTGCACATCATTGATCAATTTTTTCATTATTTTGCTCCAGCaa                              >  1:3396900/1‑150 (MQ=255)
                                                                                                                        tttCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATGCAGTGTCAGCGATGGATGCGCTTTCGCCAGTCCTGCCAGTTGTTCGTCCAGTAAGTCTTGCACATCATTGATCAATTTTTTCATTATTTTGCTCCAGCac                              >  1:1818553/1‑149 (MQ=255)
                                                                                                                         ttCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATGCAGTGTCAGCGATGGATGCGCTTTCGCCAGTCCTGCCAGTTGTTCGTCCAGTACGTCTTGCACATCATTGATCAATTTTTTCATTATTTTGCTCCAGCAAt                             <  2:1292994/150‑1 (MQ=255)
                                                                                                                          tCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATGCAGTGTCAGCGATGGATGCGCTTTCGCCAGTCCTGCCAGTTGTTCGTCCAGTACGTCTTGCACATCATTGATCAATTTTTTCATTATTTTGCTCCAGCAAtt                            <  2:1619509/150‑1 (MQ=255)
                                                                                                                           ccTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATGCAGTGTCAGCGATGGATGCGCTTTCGCCAGTCCTGCCAGTTGTTCGTCCAGTACGTCTTGCACATCATTGATCAATTTTTTCATTATTTTGCTCCAGCAATTa                           <  1:4424901/150‑1 (MQ=255)
                                                                                                                           ccTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATGCAGTGTCAGCGATGGATGCGCTTTCGCCAGTCCTGCCAGTTGTTCGTCCAGTACGTCTTGCACATCATTGATCAATTTTTTCATTATTTTGCTCCAGCAATTa                           <  1:2335344/150‑1 (MQ=255)
                                                                                                                              gCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATGCAGTGTCAGCGATGGATGCGCTTTCGCCAGTCCTGCCAGTTGTTCGTCCAGTACGTCTTGCACATCATTGATCAATTTTTTCATTATTTTGCTCCAGCAATTACgg                        >  2:4335459/1‑150 (MQ=255)
                                                                                                                               cAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATGCAGTGTCAGCGATGGATGCGCTTTCGCCAGTCCTGCCAGTTGTTCGTCCAGTACGTCTTGCACATCCTTGATCAATTTTTTCATTATTTTGCTCCAGCAATTACGGt                       >  2:590270/1‑150 (MQ=255)
                                                                                                                                aaCAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATGCAGTGTCAGCGATGGATGCGCTTTCGCCAGTCCTGCCAGTTGTTCGTCCAGTACGTCTTGCACATCATTGATCAATTTTTTCATTATTTTGCTCCAGCAATTACGGTa                      >  1:985260/1‑150 (MQ=255)
                                                                                                                                   aGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATGCAGTGTCAGCGATGGATGCGCTTTCGCCAGTCCTGCCAGTTGTTCGTCCAGTACGTCTTGCACATCATTGATCAATTTTTTCATTATTTTGCTCCAGCAATTACGGTAggg                   >  1:2612983/1‑150 (MQ=255)
                                                                                                                                     gggCATCAGCTCGGGTGACATACACCGGATCCTGATGCAGTGTCAGCGATGGATGCGCTTTCGCCAGTCCTGCCAGTTGTTCGTCCAGTACGTCTTGCACATCATTGATCAATTTTTTCATTATTTTGCTCCAGCAATTACGGTAGGGCa                 <  1:2432052/150‑1 (MQ=255)
                                                                                                                                     gggCATCAGCTCGGGTGACATACACCGGATCCTGATGCAGTGTCAGCGATGGATGCGCTTTCGCCAGTCCTGCCAGTTGTTCGTCCAGTACGTCTTGCACATCATTGATCAATTTTTTCATTATTTTGCTCCAGCAATTACGGTAGGGCa                 <  1:2223098/150‑1 (MQ=255)
                                                                                                                                      ggCATCAGCTCGGGTAACATACACCGGATCCTGATGCAGTGTCAGCGATGGATGCGCTTTCGCCAGTCCTGCCAGTTGTTCGTCCAGTACGTCTTGCACATCATTGATCAATTTTTTCATTATTTTGCTCCAGCAATTACGGTAGGGCAt                <  1:3067908/150‑1 (MQ=255)
                                                                                                                                       gCATCAGCTCGGGTGACATACACCGGATCCTGATGCAGTGTCAGCGATGGATGCGCTTTCGCCAGTCCTGCCAGTTGTTCGTCCAGTACGTCTTGCACATCATTGATCAATTTTTTCATTATTTTGCTCCAGCAATTACGGTAGGGCATg               <  1:3221018/150‑1 (MQ=255)
                                                                                                                                            aGCTCGGGTGACATACACCGGATCCTGATGCAGTGTCAGCGATGGATGCGCTTTCGCCAGTCCTGCCAGTTGTTCGTCCAGTATGTCTTGCACATCATTGATCAATTTTTTCATTATTTTGCTCCAGCAATTACGGTAGGGCATGgatga          >  1:2318770/1‑150 (MQ=255)
                                                                                                                                            aGCTCGGGTGACATACACCGGATCCTGATGCAGTGTCAGCGATGGATGCGCTTTCGCCAGTCCTGCCAGTTGTTCGTCCAGTACGTCTTGCACATCATTGATCAATTTTTTCATTATTTTGCTCCAGCAATTACGGTAGGGCATGgatga          >  2:1247269/1‑150 (MQ=255)
                                                                                                                                            aGCTCGGGTGACATACACCGGATCCTGATGCAGTGTCAGCGATGGATGCGCTTTCGCCAGTCCTGCCAGTTGTTCGTCCAGTACGTCTAGCACATCATTGATCAATTTTTTCATTATTTTGCTCCAGCAATTACGGTAGGGCATGgatga          <  2:3825101/150‑1 (MQ=255)
                                                                                                                                            aGCTCGGGTAACATACACCGGATCCTGATGCAGTGTCAGCGATGGATGCGCTTTCGCCAGTCCTGCCAGTTGTTCGTCCAGTACGTCTTGCACATCATTGATCAATTTTTTCATTATTTTGCTCCAGCAATTACGGTAGGGCATGgatga          <  1:884846/150‑1 (MQ=255)
                                                                                                                                             gCTCGGGTAACATACACCGGATCCTGATGCAGTGTCAGCGATGGATGCGCTTTCGCCAGTCCTGCCAGTTGTTCGTCCAGTACGTCTTGCACATCATTGATCAATTTTTTCATTATTTTGCTCCAGCAATTACGGTAGGGCATGgatgat         <  1:1868884/150‑1 (MQ=255)
                                                                                                                                              cTCGGGTGACATACACCGGATCCTGATGCAGTGTCAGCGATGGATGCGCTTTCGCCAGTCCTGCCAGTTGTTCGTCCAGTACGTCTTGCACATCATTGATCAATTTTTTCATTATTTTGCTCCAGCAATTACGGTAGGGCATGGatgatg        <  1:3573011/150‑1 (MQ=255)
                                                                                                                                               tCGGGTGACATACACCGGATCCTGATGCAGTGTCAGCGATGGATGCGCTTTCGCCAGTCCTGCCAGTTGTTCGTCCAGTACGTCTTGCACATCATTGATCAATTTTTTCATTATTTTGCTCCAGCAATTACGGTAGGGCATGGATGATGt       <  2:4809893/150‑1 (MQ=255)
                                                                                                                                                cGGGTAACATACACCGGATCCTGATGCAGTGTCAGCGATGGATGCGCTTTCGCCAGTCCTGCCAGTTGTTCGTCCAGTACGTCTTGCACATCATTGATCAATTTTTTCATTATTTTGCTCCAGCACTTACGGTCGGGCATGGATGATGtt      >  2:1218716/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATACACCGGATCCTGATGCAGTGTCAGCGATGGATGCGCTTTCGCCAGTCCTGCCAGTTGTTCGTCCAGTACGTCTTGCACATCATTGATCAATTTTTTCATTATTTTGCTCCAGCAATTACGGTAGGGCATGGATGATGTTCAAC  >  CP009273/1246041‑1246338

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: