Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,297,846 C→A 94.4% M150I (ATG→ATT alsE ← allulose‑6‑phosphate 3‑epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,297,8460CA94.4% 829.8 / 36.0 269M150I (ATG→ATTalsEallulose‑6‑phosphate 3‑epimerase
Reads supporting (aligned to +/- strand):  ref base C (5/10);  new base A (147/105);  total (153/116)
Fisher's exact test for biased strand distribution p-value = 6.53e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.16e-01

Only 100 of 269 total aligned reads displayed.

CAGGCCGGAAGTGCCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGAT  >  CP009273/4297699‑4297992
                                                                                                                                                   |                                                                                                                                                  
cAGGCCGGAAGTGCCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAAt                                                                                                                                                  >  1:1759209/1‑150 (MQ=255)
cAGGCCGGAAGTGCCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAAATTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAAt                                                                                                                                                  >  2:1590217/1‑150 (MQ=255)
 aGGCCGGAAGTGCCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAAtt                                                                                                                                                 >  2:4059928/1‑150 (MQ=255)
 aGGCCGGAAGTGCCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAAtt                                                                                                                                                 >  1:773210/1‑150 (MQ=255)
 aGGCCGGAAGTGCCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAAtt                                                                                                                                                 >  1:773225/1‑150 (MQ=255)
 aGGCCGGAAGTGCCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAAtt                                                                                                                                                 >  2:2909326/1‑150 (MQ=255)
    ccGGAAGTGCCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCa                                                                                                                                              <  2:2311884/150‑1 (MQ=255)
      ggAAGTGCCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAgg                                                                                                                                            >  1:4114224/1‑150 (MQ=255)
       gAAGTGCCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAATTTTATCCAGAATTTCAGGa                                                                                                                                           <  1:1630232/150‑1 (MQ=255)
       gAAGTGCCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGa                                                                                                                                           >  1:3955533/1‑150 (MQ=255)
         aGTGCCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAAt                                                                                                                                         >  1:1083834/1‑150 (MQ=255)
               gacgaTAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGGTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACgg                                                                                                                                   >  2:4458227/1‑150 (MQ=255)
                acgaTAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGt                                                                                                                                  <  2:3989220/150‑1 (MQ=255)
                acgaTAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGt                                                                                                                                  >  2:1177519/1‑150 (MQ=255)
                 cgaTAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGtt                                                                                                                                 <  1:2551064/150‑1 (MQ=255)
                  gaTAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTg                                                                                                                                >  2:2146568/1‑150 (MQ=255)
                   aTAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGt                                                                                                                               <  2:3124924/150‑1 (MQ=255)
                   aTAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGt                                                                                                                               <  2:844734/150‑1 (MQ=255)
                    tAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTc                                                                                                                              <  1:1464895/150‑1 (MQ=255)
                     aaaGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTcc                                                                                                                             >  2:1886821/1‑150 (MQ=255)
                      aaGACATCCGCCCCTGCCGCCATCAGTTTTTCGTATGTTGCCTGGTTGCAGGAACCCTACACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACTGTTGTCCg                                                                                                                            <  2:4406835/150‑1 (MQ=255)
                         aCATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCa                                                                                                                         >  1:1392210/1‑150 (MQ=255)
                           aTCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCaaa                                                                                                                       >  1:4350639/1‑150 (MQ=255)
                           aTCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCaaa                                                                                                                       >  1:4350550/1‑150 (MQ=255)
                           aTCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCaaa                                                                                                                       >  1:2380364/1‑150 (MQ=255)
                              cGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGcc                                                                                                                    >  1:1015259/1‑150 (MQ=255)
                              cGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGcc                                                                                                                    <  2:3751910/150‑1 (MQ=255)
                                 cccTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCggg                                                                                                                 >  1:3436352/1‑150 (MQ=255)
                                   cTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGAt                                                                                                               <  2:2376475/150‑1 (MQ=255)
                                   cTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGAt                                                                                                               <  2:2376585/150‑1 (MQ=255)
                                    tGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATc                                                                                                              <  2:269747/150‑1 (MQ=255)
                                    tGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATc                                                                                                              >  2:4199055/1‑150 (MQ=255)
                                        gccATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACa                                                                                                          <  2:1862996/150‑1 (MQ=255)
                                              aGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATg                                                                                                    <  2:670499/150‑1 (MQ=255)
                                               gTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGa                                                                                                   <  1:4036249/150‑1 (MQ=255)
                                                tttttCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGAc                                                                                                  >  1:3724515/1‑150 (MQ=255)
                                                tttttCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGAc                                                                                                  >  1:2822930/1‑150 (MQ=255)
                                                tttttCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGAc                                                                                                  >  1:538527/1‑150 (MQ=255)
                                                 ttttCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGAcc                                                                                                 <  2:1083828/150‑1 (MQ=255)
                                                    tCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTa                                                                                              >  1:2449166/1‑150 (MQ=255)
                                                     cGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTaa                                                                                             >  2:2918896/1‑150 (MQ=255)
                                                       tAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAATTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAAtt                                                                                           <  1:3069862/150‑1 (MQ=255)
                                                         aGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAAtttt                                                                                         >  2:1691523/1‑150 (MQ=255)
                                                         aGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAAtttt                                                                                         <  2:2598039/150‑1 (MQ=255)
                                                           ttGCCTGGTTGCATGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTAt                                                                                       >  1:1333064/1‑150 (MQ=255)
                                                           ttGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTAt                                                                                       <  2:337398/150‑1 (MQ=255)
                                                           ttGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTAt                                                                                       <  2:2633646/150‑1 (MQ=255)
                                                              ccTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCgg                                                                                    <  1:291937/150‑1 (MQ=255)
                                                                 ggTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCt                                                                                 <  1:4191755/150‑1 (MQ=255)
                                                                 ggTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCt                                                                                 >  2:2565666/1‑150 (MQ=255)
                                                                     gCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATg                                                                             >  1:4444245/1‑150 (MQ=255)
                                                                       aGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGa                                                                           >  1:1665335/1‑150 (MQ=255)
                                                                       aGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGa                                                                           >  1:2432253/1‑150 (MQ=255)
                                                                       aGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGa                                                                           <  1:1967303/150‑1 (MQ=255)
                                                                       aGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGa                                                                           >  2:1285509/1‑150 (MQ=255)
                                                                       aGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGa                                                                           >  2:1285721/1‑150 (MQ=255)
                                                                        ggAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGat                                                                          >  2:468910/1‑150 (MQ=255)
                                                                        ggAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCGTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGat                                                                          >  2:468689/1‑150 (MQ=255)
                                                                           aCCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGAtata                                                                       <  2:2339650/150‑1 (MQ=255)
                                                                             cGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAATTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGt                                                                     >  2:3589708/1‑150 (MQ=255)
                                                                               tCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTAt                                                                   >  1:598613/1‑150 (MQ=255)
                                                                                   ccTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCa                                                               <  2:2563951/150‑1 (MQ=255)
                                                                                    cTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCAt                                                              <  2:1727440/150‑1 (MQ=255)
                                                                                       aaTTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGc                                                           <  2:2978218/150‑1 (MQ=255)
                                                                                        aTTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGTTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGcc                                                          >  1:2932717/1‑150 (MQ=255)
                                                                                        aTTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGTATATAGTATTTCATGGcc                                                          >  1:2976231/1‑150 (MQ=255)
                                                                                          ttCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTc                                                        <  1:4403415/150‑1 (MQ=255)
                                                                                          ttCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTc                                                        <  1:3713648/150‑1 (MQ=255)
                                                                                           tCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCa                                                       >  1:2243963/1‑150 (MQ=255)
                                                                                            cGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCaa                                                      >  1:2747957/1‑150 (MQ=255)
                                                                                            cGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCaa                                                      <  2:773223/150‑1 (MQ=255)
                                                                                            cGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGACTCat                                                      <  2:773885/150‑2 (MQ=255)
                                                                                             gTACTCCAGACCTTCTCGTACACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGACGTGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAAc                                                     >  2:2867010/1‑150 (MQ=255)
                                                                                                cTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTgg                                                  >  2:1249293/1‑150 (MQ=255)
                                                                                                  ccAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCg                                                <  1:954322/150‑1 (MQ=255)
                                                                                                  ccAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGACATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCg                                                <  1:4458227/150‑1 (MQ=255)
                                                                                                     gACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGtct                                             <  1:705895/150‑1 (MQ=255)
                                                                                                      aCCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGtctc                                            <  2:4432958/150‑1 (MQ=255)
                                                                                                       ccTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTcc                                           <  2:3941747/150‑1 (MQ=255)
                                                                                                         ttCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCgg                                         <  2:3249251/150‑1 (MQ=255)
                                                                                                         ttCTCGGTCACGCCATGCCTTCAGTTCGGCAATTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCgg                                         <  1:4000536/150‑1 (MQ=255)
                                                                                                           ctcGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGt                                       <  1:3610416/150‑1 (MQ=255)
                                                                                                                   cGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATc                               <  1:2339600/150‑1 (MQ=255)
                                                                                                                     ccATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAg                             >  2:1549181/1‑150 (MQ=255)
                                                                                                                       aTGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGcc                           <  1:3200900/150‑1 (MQ=255)
                                                                                                                         gCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGcccc                         <  2:3574130/150‑1 (MQ=255)
                                                                                                                         gCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGcccc                         >  1:4167292/1‑150 (MQ=255)
                                                                                                                          ccTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCa                        >  1:2320506/1‑150 (MQ=255)
                                                                                                                            ttCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACt                      <  2:3219046/150‑1 (MQ=255)
                                                                                                                                   cGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTtcatgt               >  2:1889034/1‑150 (MQ=255)
                                                                                                                                    ggCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTtcatgtc              <  1:1890683/150‑1 (MQ=255)
                                                                                                                                     gCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTtcatgtca             >  1:4208576/1‑150 (MQ=255)
                                                                                                                                     gCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTtcatgtca             <  1:934862/150‑1 (MQ=255)
                                                                                                                                       aaGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTtcatgtcatg           <  1:2147278/150‑1 (MQ=255)
                                                                                                                                       aaGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTtcatgtcatg           >  1:4316931/1‑150 (MQ=255)
                                                                                                                                          tttATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTCCATGTCATGAcg        <  2:3888964/150‑1 (MQ=255)
                                                                                                                                            tATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGAcggc      >  1:557615/1‑150 (MQ=255)
                                                                                                                                            tATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGAcggc      <  2:2500933/150‑1 (MQ=255)
                                                                                                                                             aTCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGAcggcg     >  1:3910177/1‑150 (MQ=255)
                                                                                                                                                cAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGAt  <  1:1230105/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CAGGCCGGAAGTGCCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGAT  >  CP009273/4297699‑4297992

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: