Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,297,846 C→A 9.0% 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,8460CA 9.0% 1293.4 / 78.7 348M150I (ATG→ATTalsEallulose‑6‑phosphate 3‑epimerase
Reads supporting (aligned to +/- strand):  ref base C (154/162);  new base A (16/16);  total (170/178)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 350 total aligned reads displayed.

ACAGGCCGGAAGTGCCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTC  >  CP009273/4297698‑4297995
                                                                                                                                                    |                                                                                                                                                      
aCAGGCCGGAAGTGCCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCa                                                                                                                                                       >  1:3105391/1‑150 (MQ=255)
 cAGGCCGGAAGTGCCGACGATAATGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCAt                                                                                                                                                      <  2:1153856/150‑1 (MQ=255)
 cAGGCCGGAAGTGCCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCAt                                                                                                                                                      <  1:4921359/150‑1 (MQ=255)
 cAGGCCGGAAGTGCCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCAt                                                                                                                                                      <  1:2199985/150‑1 (MQ=255)
  aGGCCGGAAGTGCCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCAtt                                                                                                                                                     <  1:3887811/150‑1 (MQ=255)
   ggCCGGAAGTGCCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCAttt                                                                                                                                                    <  1:3524618/150‑1 (MQ=255)
   ggCCGGAAGTGCCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCAttt                                                                                                                                                    <  1:4653524/150‑1 (MQ=255)
    gCCGGAAGTGCCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTc                                                                                                                                                   >  1:4993700/1‑150 (MQ=255)
     ccGGAAGTGCCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCa                                                                                                                                                  <  1:4784451/150‑1 (MQ=255)
     ccGGAAGTGCCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCa                                                                                                                                                  <  1:4367690/150‑1 (MQ=255)
       ggAAGTGCCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAgg                                                                                                                                                <  2:4767109/150‑1 (MQ=255)
        gAAGTGCCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCCGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGa                                                                                                                                              >  2:4899537/1‑150 (MQ=255)
         aaGTGCCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTATCGTAAGTTGCCTGGTTGCAGGACCCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTAGGCAAGTTTCTCCAGCATTTCAGGaa                                                                                                                                             >  1:4682086/1‑150 (MQ=255)
            tGCCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGa                                                                                                                                          <  2:2837530/150‑1 (MQ=255)
                 acgaTAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGt                                                                                                                                     >  1:1461020/1‑150 (MQ=255)
                    aTAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGt                                                                                                                                  >  1:3367767/1‑150 (MQ=255)
                        aGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCgg                                                                                                                              >  1:1838529/1‑150 (MQ=255)
                         gACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGc                                                                                                                             >  2:4300122/1‑150 (MQ=255)
                          aCATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCa                                                                                                                            >  1:3702392/1‑150 (MQ=255)
                          aCATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCa                                                                                                                            <  2:195619/150‑1 (MQ=255)
                            aTCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCaaa                                                                                                                          >  1:5077493/1‑150 (MQ=255)
                                   ccTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGa                                                                                                                   <  2:4422277/150‑1 (MQ=255)
                                     tGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATc                                                                                                                 <  2:827789/150‑1 (MQ=255)
                                     tGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATc                                                                                                                 <  1:594009/150‑1 (MQ=255)
                                     tGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATc                                                                                                                 <  1:2650140/150‑1 (MQ=255)
                                     tGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATc                                                                                                                 <  2:3354133/150‑1 (MQ=255)
                                        cgccATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGAc                                                                                                              <  2:1270667/150‑1 (MQ=255)
                                           cATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGt                                                                                                           <  2:80666/150‑1 (MQ=255)
                                           cATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGt                                                                                                           >  1:3730560/1‑150 (MQ=255)
                                            aTCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTc                                                                                                          >  2:3596513/1‑150 (MQ=255)
                                                 tttttCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGAc                                                                                                     >  2:5209888/1‑150 (MQ=255)
                                                 tttttCGTAAGTGGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGAc                                                                                                     <  2:3003081/150‑1 (MQ=255)
                                                       gTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAAt                                                                                               <  2:1835872/150‑1 (MQ=255)
                                                       gTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAAt                                                                                               <  1:3907844/150‑1 (MQ=255)
                                                       gTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACACCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGTATCGACAGTCATGACCGTAAt                                                                                               <  1:3184019/150‑1 (MQ=255)
                                                          aGTTGCCTGGTTGCAGGAACCGTCCACCTCAATGTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAAtttt                                                                                            <  2:4269142/150‑1 (MQ=255)
                                                            ttGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCCGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTAt                                                                                          <  2:1041813/150‑1 (MQ=255)
                                                            ttGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTAt                                                                                          <  2:2323/150‑1 (MQ=255)
                                                             tGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATc                                                                                         <  2:3277200/150‑1 (MQ=255)
                                                             tGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATc                                                                                         <  1:5308988/150‑1 (MQ=255)
                                                              gCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCg                                                                                        <  1:1324798/150‑1 (MQ=255)
                                                               ccTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCgg                                                                                       <  2:1580850/150‑1 (MQ=255)
                                                               ccTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCgg                                                                                       >  1:962352/1‑150 (MQ=255)
                                                                cTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGc                                                                                      <  1:847231/150‑1 (MQ=255)
                                                                   gTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCtt                                                                                   >  2:1325756/1‑150 (MQ=255)
                                                                    ttGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTa                                                                                  >  2:2162291/1‑150 (MQ=255)
                                                                    ttGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTa                                                                                  <  1:1145671/150‑1 (MQ=255)
                                                                       cAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATgg                                                                               <  1:75922/150‑1 (MQ=255)
                                                                         ggAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGat                                                                             >  1:5152961/1‑150 (MQ=255)
                                                                         ggAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGat                                                                             <  1:2059478/150‑1 (MQ=255)
                                                                          gAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGcta                                                                            >  2:3810059/1‑150 (MQ=255)
                                                                          gAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGata                                                                            >  2:2705881/1‑150 (MQ=255)
                                                                            aCCGTCCACCTCAATTTCGTACGCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAATGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTAGGGAtata                                                                          <  2:3530725/150‑1 (MQ=255)
                                                                              cGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGt                                                                        <  1:5444209/150‑1 (MQ=255)
                                                                               gTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTa                                                                       >  2:3885675/1‑150 (MQ=255)
                                                                                tCCACCTCAAGTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTAt                                                                      <  1:2892220/150‑1 (MQ=255)
                                                                                 ccACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGAAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTAtt                                                                     <  2:2542811/150‑1 (MQ=255)
                                                                                   aCCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTc                                                                   >  2:2760200/1‑150 (MQ=255)
                                                                                     cTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCAt                                                                 <  2:2249604/150‑1 (MQ=255)
                                                                                      tCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCTTGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATg                                                                <  2:1896409/150‑1 (MQ=255)
                                                                                      tCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCAtt                                                                >  2:4187240/1‑149 (MQ=255)
                                                                                      tCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATg                                                                <  1:1958245/150‑1 (MQ=255)
                                                                                      tCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATg                                                                <  1:3421402/150‑1 (MQ=255)
                                                                                        aaTTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGc                                                              >  2:4600697/1‑150 (MQ=255)
                                                                                        aaTTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGc                                                              >  2:1097080/1‑150 (MQ=255)
                                                                                         aTTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGcc                                                             >  1:1072762/1‑150 (MQ=255)
                                                                                            tCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCa                                                          <  1:234396/150‑1 (MQ=255)
                                                                                            tCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCa                                                          >  1:4482995/1‑150 (MQ=255)
                                                                                            tCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCa                                                          <  2:164315/150‑1 (MQ=255)
                                                                                                  tCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGc                                                    <  1:4296928/150‑1 (MQ=255)
                                                                                                     aGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGtc                                                 <  2:3631180/150‑1 (MQ=255)
                                                                                                      gACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGtct                                                >  1:5273767/1‑150 (MQ=255)
                                                                                                       aCCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGtctc                                               >  2:2448087/1‑150 (MQ=255)
                                                                                                       aCCATCTCGTTCACGCCATGACTTCAGAAAGGAAAGTTTATCCAGCATTACAGGAATCAACGGTTGTCCGGCAAAGCCGGGATCGACATTCATGACCTTAATTATAAAGGACTTATGGATATAGAATATCATGGCATCAACTCGCGtctc                                               >  2:5162104/1‑150 (MQ=255)
                                                                                                        ccTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTcc                                              <  1:1407235/150‑1 (MQ=255)
                                                                                                        ccTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTcc                                              <  1:1563870/150‑1 (MQ=255)
                                                                                                         cTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCg                                             <  1:3591262/150‑1 (MQ=255)
                                                                                                           tctcGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCggg                                           >  1:624183/1‑150 (MQ=255)
                                                                                                               gTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTaa                                       >  1:2786708/1‑150 (MQ=255)
                                                                                                                ttCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAg                                      >  2:556083/1‑150 (MQ=255)
                                                                                                                    cGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATc                                  >  1:4778523/1‑150 (MQ=255)
                                                                                                                        aTGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGcc                              <  2:4510903/150‑1 (MQ=255)
                                                                                                                          gCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGcccc                            >  2:5033147/1‑150 (MQ=255)
                                                                                                                            cTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCAc                          >  2:619033/1‑150 (MQ=255)
                                                                                                                            cTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCAc                          >  2:2487552/1‑150 (MQ=255)
                                                                                                                              tCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACtt                        <  1:2426766/150‑1 (MQ=255)
                                                                                                                                aGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTtc                      <  2:1774468/150‑1 (MQ=255)
                                                                                                                                      gCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTtcatgtca                <  1:3299405/150‑1 (MQ=255)
                                                                                                                                      gCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTtcatgtca                >  1:2166873/1‑150 (MQ=255)
                                                                                                                                      gCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTtcatgtca                >  1:4135615/1‑150 (MQ=255)
                                                                                                                                        aaGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTtcatgtcatg              >  1:3762539/1‑150 (MQ=255)
                                                                                                                                        aaGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTtcatgtcatg              <  1:4038859/150‑1 (MQ=255)
                                                                                                                                          gTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGAc            >  2:5424419/1‑150 (MQ=255)
                                                                                                                                              aTCCAGCAGTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGAcggcg        <  1:3278871/150‑1 (MQ=255)
                                                                                                                                                ccATCATTTCATGAATGAACGGTTGTCCGGCAAAGCCGGGATAGACTGTCATGACCGTAATTTTATCGTCCTTATGGATTTATTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGa      <  1:1745122/150‑1 (MQ=255)
                                                                                                                                                ccAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGa      >  2:317123/1‑150 (MQ=255)
                                                                                                                                                  aGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGAtt    >  2:2024000/1‑150 (MQ=255)
                                                                                                                                                   gCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGAttt   <  1:2571690/150‑1 (MQ=255)
                                                                                                                                                   gCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGAttt   >  1:736988/1‑150 (MQ=255)
                                                                                                                                                    cATTTCAGGAAATGAACGGTTGTCCGGCAAAGCCGGGATAGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGAttt   <  2:3702521/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
ACAGGCCGGAAGTGCCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTC  >  CP009273/4297698‑4297995

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: