Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,297,846 C→A 92.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,8460CA92.0% 936.7 / 70.5 311M150I (ATG→ATTalsEallulose‑6‑phosphate 3‑epimerase
Reads supporting (aligned to +/- strand):  ref base C (14/11);  new base A (134/152);  total (148/163)
Fisher's exact test for biased strand distribution p-value = 4.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.66e-01

Only 100 of 312 total aligned reads displayed.

AACAGGCCGGAAGTGCCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTT  >  CP009273/4297697‑4297994
                                                                                                                                                     |                                                                                                                                                    
  cAGGCCGGAAGTGCCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAAt                                                                                                                                                    <  2:571837/150‑1 (MQ=255)
   aGGCCGGAAGTGCCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAAtt                                                                                                                                                   <  2:1023263/150‑1 (MQ=255)
   aGGCCGGAAGTGCCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAAtt                                                                                                                                                   <  1:4404140/150‑1 (MQ=255)
   aGGCCGGAAGTGCCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAAtt                                                                                                                                                   >  1:2707042/1‑150 (MQ=255)
       cGGAAGTGCCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAg                                                                                                                                               >  2:4722738/1‑150 (MQ=255)
       cGGAAGTGCCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAg                                                                                                                                               >  1:2860884/1‑150 (MQ=255)
        ggAAGTGCCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAgg                                                                                                                                              >  1:1006118/1‑150 (MQ=255)
         gAAGTGCCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGa                                                                                                                                             >  2:2569546/1‑150 (MQ=255)
          aaGTGCCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGaa                                                                                                                                            <  1:523969/150‑1 (MQ=255)
             tGCCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGa                                                                                                                                         <  2:1550743/150‑1 (MQ=255)
                    gaTAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTg                                                                                                                                  >  2:1629719/1‑150 (MQ=255)
                    gaTAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTg                                                                                                                                  >  1:619794/1‑150 (MQ=255)
                    gaTAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTg                                                                                                                                  <  1:1166212/150‑1 (MQ=255)
                      tAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTc                                                                                                                                <  2:476172/150‑1 (MQ=255)
                       aaaGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTcc                                                                                                                               <  2:2060872/150‑1 (MQ=255)
                       aaaGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTcc                                                                                                                               <  2:2141556/150‑1 (MQ=255)
                        aaGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCg                                                                                                                              >  1:1772753/1‑150 (MQ=255)
                              tCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCGGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAg                                                                                                                        >  1:1947469/1‑150 (MQ=255)
                                cGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGcc                                                                                                                      <  2:1731179/150‑1 (MQ=255)
                                 gCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCg                                                                                                                     >  1:3633303/1‑150 (MQ=255)
                                  ccccTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCgg                                                                                                                    >  1:2230560/1‑150 (MQ=255)
                                   cccTGCCGCCATGAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAATTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCggg                                                                                                                   <  2:1293479/150‑1 (MQ=255)
                                    ccTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGa                                                                                                                  <  2:2978968/150‑1 (MQ=255)
                                     cTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGAt                                                                                                                 <  1:2992056/150‑1 (MQ=255)
                                      tGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATc                                                                                                                <  2:4495500/150‑1 (MQ=255)
                                      tGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATc                                                                                                                <  1:2134053/150‑1 (MQ=255)
                                      tGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATGTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATc                                                                                                                <  1:1425809/150‑1 (MQ=255)
                                       gccgccATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCg                                                                                                               >  1:3653796/1‑150 (MQ=255)
                                           ccATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAATTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAg                                                                                                           >  2:1854197/1‑150 (MQ=255)
                                            cATCAGTTTTTCGTATGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGAAATTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGt                                                                                                          <  1:1600415/150‑1 (MQ=255)
                                             aTCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTc                                                                                                         >  2:746568/1‑150 (MQ=255)
                                             aTCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTc                                                                                                         >  1:1296737/1‑150 (MQ=255)
                                              tCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCa                                                                                                        >  1:2779448/1‑150 (MQ=255)
                                               cAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCAt                                                                                                       >  2:2379120/1‑150 (MQ=255)
                                               cAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCAt                                                                                                       >  1:4097798/1‑150 (MQ=255)
                                               cAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCAt                                                                                                       <  1:4690642/150‑1 (MQ=255)
                                                aGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATg                                                                                                      <  2:3698490/150‑1 (MQ=255)
                                                 gTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGa                                                                                                     <  2:3496293/150‑1 (MQ=255)
                                                   ttttCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGAcc                                                                                                   <  1:50833/150‑1 (MQ=255)
                                                      tCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTa                                                                                                >  2:2176202/1‑150 (MQ=255)
                                                       cGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTaa                                                                                               <  1:4355963/150‑1 (MQ=255)
                                                       cGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTaa                                                                                               >  1:3814301/1‑150 (MQ=255)
                                                           aGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAAtttt                                                                                           >  2:3983701/1‑150 (MQ=255)
                                                              tGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATc                                                                                        <  1:3959847/150‑1 (MQ=255)
                                                                ccTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTACAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCgg                                                                                      >  1:3613262/1‑150 (MQ=255)
                                                                   ggTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCt                                                                                   <  1:4588226/150‑1 (MQ=255)
                                                                      tGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTAt                                                                                <  1:4117180/150‑1 (MQ=255)
                                                                         aGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGa                                                                             >  2:2920483/1‑150 (MQ=255)
                                                                         aGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGa                                                                             >  1:2155053/1‑150 (MQ=255)
                                                                           gAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGata                                                                           <  1:3035620/150‑1 (MQ=255)
                                                                           gAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGata                                                                           <  1:786562/150‑1 (MQ=255)
                                                                            aaCCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGatat                                                                          >  2:1814006/1‑150 (MQ=255)
                                                                            aaCCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGagat                                                                          >  2:3027173/1‑150 (MQ=255)
                                                                             aCCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGAtata                                                                         <  2:1896071/150‑1 (MQ=255)
                                                                              ccGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAg                                                                        >  2:2265413/1‑150 (MQ=255)
                                                                               cGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGt                                                                       <  1:4720040/150‑1 (MQ=255)
                                                                                gTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTa                                                                      <  2:2797942/150‑1 (MQ=255)
                                                                                gTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTa                                                                      >  2:3199896/1‑150 (MQ=255)
                                                                                  ccACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTAtt                                                                    >  2:5044587/1‑150 (MQ=255)
                                                                                   cACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATTCCTTAAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTAGCGGCCTTATGGATATAGTAttt                                                                   >  1:4880860/1‑150 (MQ=255)
                                                                                   cACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAATTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTAttt                                                                   <  2:5092691/150‑1 (MQ=255)
                                                                                   cACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTAttt                                                                   <  2:1083737/150‑1 (MQ=255)
                                                                                   cACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTAttt                                                                   >  1:1531326/1‑150 (MQ=255)
                                                                                       tCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATg                                                               <  1:712927/150‑1 (MQ=255)
                                                                                        cAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATgg                                                              <  2:4964114/150‑1 (MQ=255)
                                                                                         aaTTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGc                                                             >  1:4775098/1‑150 (MQ=255)
                                                                                          aTTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGAc                                                                                                    <  1:2310439/110‑1 (MQ=255)
                                                                                            ttCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTc                                                          <  2:1245162/150‑1 (MQ=255)
                                                                                             tCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCa                                                         <  2:38183/150‑1 (MQ=255)
                                                                                             tCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCa                                                         <  2:3433764/150‑1 (MQ=255)
                                                                                               gTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAAc                                                       >  2:5086643/1‑150 (MQ=255)
                                                                                                tACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGTAATTAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACt                                                      <  2:4313904/150‑1 (MQ=255)
                                                                                                tACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACt                                                      >  2:2853694/1‑150 (MQ=255)
                                                                                                    ccAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAATTTTATCCAGAATTTCAGGAAGGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCg                                                  <  2:4530925/150‑1 (MQ=255)
                                                                                                     cAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGt                                                 >  2:2729812/1‑150 (MQ=255)
                                                                                                     cAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGt                                                 >  2:4633383/1‑150 (MQ=255)
                                                                                                     cAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGt                                                 >  2:2309740/1‑150 (MQ=255)
                                                                                                           ttCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCgg                                           <  1:1693899/150‑1 (MQ=255)
                                                                                                              tcGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGtt                                        >  2:4329661/1‑150 (MQ=255)
                                                                                                                       ccATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAg                               >  1:3510838/1‑150 (MQ=255)
                                                                                                                       ccATGCCTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAg                               <  1:4515642/150‑1 (MQ=255)
                                                                                                                             cTTCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCAc                         <  2:4379576/150‑1 (MQ=255)
                                                                                                                              ttCAGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACt                        <  1:1556491/150‑1 (MQ=255)
                                                                                                                                 aGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTtc                     <  2:4069842/150‑1 (MQ=255)
                                                                                                                                 aGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTtc                     >  2:3379440/1‑150 (MQ=255)
                                                                                                                                 aGTTCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTtc                     <  2:1156521/150‑1 (MQ=255)
                                                                                                                                   ttCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTtcat                   >  2:2881948/1‑150 (MQ=255)
                                                                                                                                    tCGGCAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTtcatg                  >  1:2059992/1‑150 (MQ=255)
                                                                                                                                        cAAGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTtcatgtcat              >  2:2655254/1‑150 (MQ=255)
                                                                                                                                         aaGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTtcatgtcatg             <  1:3102864/150‑1 (MQ=255)
                                                                                                                                          aGTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGa            >  1:10719/1‑150 (MQ=255)
                                                                                                                                           gTTTATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGAc           <  2:3571077/150‑1 (MQ=255)
                                                                                                                                            tttATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGAcg          <  1:4174630/150‑1 (MQ=255)
                                                                                                                                            tttATCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGAcg          <  1:3682085/150‑1 (MQ=255)
                                                                                                                                               aTCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGAcggcg       <  2:619794/150‑1 (MQ=255)
                                                                                                                                                tCCAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGAcggcgg      <  1:3273678/150‑1 (MQ=255)
                                                                                                                                                 ccAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGa     >  1:2444441/1‑150 (MQ=255)
                                                                                                                                                  cAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGAt    <  2:2707042/150‑1 (MQ=255)
                                                                                                                                                  cAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGAt    <  1:499305/150‑1 (MQ=255)
                                                                                                                                                  cAGAATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGAt    >  2:84201/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AACAGGCCGGAAGTGCCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTT  >  CP009273/4297697‑4297994

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: