Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,557,815 T→A 5.8% K709N (AAA→AAT)  glgP ← glycogen phosphorylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,557,8150TA 5.8% 874.9 / 34.0 275K709N (AAA→AAT) glgPglycogen phosphorylase
Reads supporting (aligned to +/- strand):  ref base T (120/138);  new base A (10/6);  total (131/144)
Fisher's exact test for biased strand distribution p-value = 3.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.70e-01

Only 100 of 277 total aligned reads displayed.

CGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCATG  >  CP009273/3557668‑3557962
                                                                                                                                                   |                                                                                                                                                   
 gCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGttt                                                                                                                                                   <  2:2494337/149‑1 (MQ=255)
  ccAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTg                                                                                                                                                  >  1:3517907/1‑149 (MQ=255)
  ccAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTg                                                                                                                                                  <  1:2519839/149‑1 (MQ=255)
  ccAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTg                                                                                                                                                  <  1:4259725/149‑1 (MQ=255)
   cAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGt                                                                                                                                                 <  1:262512/149‑1 (MQ=255)
   cAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGt                                                                                                                                                 <  2:4217505/149‑1 (MQ=255)
    aGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTa                                                                                                                                                >  1:89483/1‑149 (MQ=255)
        ccTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGccc                                                                                                                                            >  1:3998739/1‑149 (MQ=255)
        ccTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGccc                                                                                                                                            >  1:418047/1‑149 (MQ=255)
           ggTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTga                                                                                                                                         <  1:3780790/149‑1 (MQ=255)
           ggTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTga                                                                                                                                         >  1:3878346/1‑149 (MQ=255)
              aGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTgacga                                                                                                                                      <  2:4104651/149‑1 (MQ=255)
               gTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTgacgac                                                                                                                                     <  2:4550028/149‑1 (MQ=255)
               gTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTgacgac                                                                                                                                     >  1:187682/1‑149 (MQ=255)
                tGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGacgacg                                                                                                                                    >  1:4022230/1‑149 (MQ=255)
                tGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGacgacg                                                                                                                                    <  1:4379224/149‑1 (MQ=255)
                tGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGacgacg                                                                                                                                    >  2:1728079/1‑149 (MQ=255)
                tGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGacgacg                                                                                                                                    >  1:2797300/1‑149 (MQ=255)
                 gATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGc                                                                                                                                   >  2:694213/1‑149 (MQ=255)
                 gATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGc                                                                                                                                   >  1:3651711/1‑149 (MQ=255)
                 gATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGATTGTAGCCCTGACGACGc                                                                                                                                   <  2:1444954/149‑1 (MQ=255)
                   tCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAg                                                                                                                                 <  1:115815/149‑1 (MQ=255)
                    cGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGt                                                                                                                                >  2:2739442/1‑149 (MQ=255)
                     gCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGATTGTAGCCCTGACGACGCAGtt                                                                                                                               >  2:609243/1‑149 (MQ=255)
                        gAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGATTGTAGCCCTGACGACGCAGttctt                                                                                                                            >  2:307121/1‑149 (MQ=255)
                         aaGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGttcttc                                                                                                                           >  2:1508781/1‑149 (MQ=255)
                          aGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTcc                                                                                                                          <  1:3913128/149‑1 (MQ=255)
                           gTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCa                                                                                                                         <  2:3421889/149‑1 (MQ=255)
                           gTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCa                                                                                                                         >  1:2557019/1‑149 (MQ=255)
                           gTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTGTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCa                                                                                                                         >  1:2393982/1‑149 (MQ=255)
                            ttGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCAc                                                                                                                        >  2:3342627/1‑149 (MQ=255)
                             tGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCAct                                                                                                                       <  1:605444/149‑1 (MQ=255)
                              gATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCActt                                                                                                                      >  2:191619/1‑149 (MQ=255)
                                tCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCActtct                                                                                                                    <  2:3503739/149‑1 (MQ=255)
                                 cAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCActtctt                                                                                                                   <  1:1648767/149‑1 (MQ=255)
                                  aGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACttcttc                                                                                                                  <  1:838467/149‑1 (MQ=255)
                                   gCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTcc                                                                                                                 <  2:3940991/149‑1 (MQ=255)
                                     gAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGc                                                                                                               >  1:2194720/1‑149 (MQ=255)
                                     gAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGc                                                                                                               <  2:3691168/149‑1 (MQ=255)
                                          aaCCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTgtgt                                                                                                          <  1:1703221/149‑1 (MQ=255)
                                            ccAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTa                                                                                                        <  2:2825715/149‑1 (MQ=255)
                                                 tCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCCGCTCCTCATCTTTCTCGTAGTATTCACGCGGTGTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGGGGTACCaac                                                                                                   >  2:3715782/1‑148 (MQ=255)
                                                   gcgATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCaaaaa                                                                                                 <  1:4209736/149‑1 (MQ=255)
                                                   gcgATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCaaaaa                                                                                                 <  2:1791426/149‑1 (MQ=255)
                                                      aTAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATaa                                                                                              <  2:4686777/149‑1 (MQ=255)
                                                      aTAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATaa                                                                                              >  2:2055289/1‑149 (MQ=255)
                                                        aGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAg                                                                                            <  2:2239904/149‑1 (MQ=255)
                                                          cGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGat                                                                                          <  1:2687114/149‑1 (MQ=255)
                                                           gACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGata                                                                                         <  1:730656/149‑1 (MQ=255)
                                                            aCCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGatat                                                                                        <  2:1066428/149‑1 (MQ=255)
                                                                ggATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTc                                                                                    >  1:797837/1‑149 (MQ=255)
                                                                  aTCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAg                                                                                  <  1:899704/149‑1 (MQ=255)
                                                                        cGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGa                                                                            <  2:2582112/149‑1 (MQ=255)
                                                                         ggACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGAc                                                                           >  1:3670568/1‑149 (MQ=255)
                                                                          gACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACa                                                                          >  2:3023952/1‑149 (MQ=255)
                                                                          gACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACa                                                                          <  1:1220336/149‑1 (MQ=255)
                                                                          gACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACa                                                                          <  1:2739442/149‑1 (MQ=255)
                                                                          gACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACa                                                                          <  1:4673464/149‑1 (MQ=255)
                                                                            cTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATg                                                                        >  1:3021222/1‑149 (MQ=255)
                                                                             tGACTACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGa                                                                       <  2:4129970/149‑1 (MQ=255)
                                                                             tGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGa                                                                       >  1:2222844/1‑149 (MQ=255)
                                                                              gAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGAt                                                                      >  1:3467491/1‑149 (MQ=255)
                                                                              gAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGAt                                                                      >  2:3232278/1‑149 (MQ=255)
                                                                               aaTACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATc                                                                     >  2:80227/1‑149 (MQ=255)
                                                                               aaTACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATc                                                                     >  2:2390556/1‑149 (MQ=255)
                                                                                  acacCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAg                                                                  >  2:1983720/1‑149 (MQ=255)
                                                                                  acacCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAg                                                                  >  2:3454652/1‑149 (MQ=255)
                                                                                        cTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCTGTTCTTCCACTTCTTCCGATGTGTTACCAAAAATAAAGATATTGTCAGCACCTACATGATCCAGCAtctc                                                            <  2:2771105/149‑1 (MQ=255)
                                                                                             gATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACAt                                                       <  1:214924/149‑1 (MQ=255)
                                                                                                     tCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCg                                               <  1:4363450/149‑1 (MQ=255)
                                                                                                      cAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGt                                              <  2:3601856/149‑1 (MQ=255)
                                                                                                      cAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGt                                              <  1:350779/149‑1 (MQ=255)
                                                                                                        gCACGTGATTCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGTCGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTcc                                            <  1:2081179/149‑1 (MQ=255)
                                                                                                        gCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTcc                                            >  2:1115996/1‑149 (MQ=255)
                                                                                                        gCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTcc                                            <  2:1358729/149‑1 (MQ=255)
                                                                                                           ccTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAAc                                         <  1:2797019/149‑1 (MQ=255)
                                                                                                                  cAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGa                                  <  2:3699099/149‑1 (MQ=255)
                                                                                                                  cAGCTCCTCATCTTTCTCGTAGTATTCACGCGGATTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGa                                  <  1:307121/149‑1 (MQ=255)
                                                                                                                   aGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGAt                                 >  2:2987433/1‑149 (MQ=255)
                                                                                                                   aGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGAt                                 >  1:4120744/1‑149 (MQ=255)
                                                                                                                   aGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGAt                                 <  1:1832743/149‑1 (MQ=255)
                                                                                                                       cctcATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTc                             >  1:4587052/1‑149 (MQ=255)
                                                                                                                         tcATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAg                           <  1:154471/149‑1 (MQ=255)
                                                                                                                           aTCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAgcg                         >  2:316874/1‑149 (MQ=255)
                                                                                                                           aTCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAgcg                         >  2:316703/1‑149 (MQ=255)
                                                                                                                            tCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAgcgc                        <  1:4690505/149‑1 (MQ=255)
                                                                                                                            tCTTTCTCGTAGTATTCACGCGGATTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAgcgc                        >  1:3716749/1‑149 (MQ=255)
                                                                                                                                  tcGTAGTATTCACGCGGATTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGtt                  >  1:3266632/1‑149 (MQ=255)
                                                                                                                                    gtagtaTTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTaa                <  2:2943064/149‑1 (MQ=255)
                                                                                                                                    gtagtaTTCACGCGGATTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTaa                >  2:2456681/1‑149 (MQ=255)
                                                                                                                                     tagtaTTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAg               <  2:3295964/149‑1 (MQ=255)
                                                                                                                                        taTTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAgcgc            >  1:932990/1‑149 (MQ=255)
                                                                                                                                         aTTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAGCGCa           <  1:4317508/149‑1 (MQ=255)
                                                                                                                                         aTTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAGCGCa           <  1:1389055/149‑1 (MQ=255)
                                                                                                                                           tCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAGCGCaaa         <  2:2194720/149‑1 (MQ=255)
                                                                                                                                            cACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAGCGCAAAc        >  1:63919/1‑149 (MQ=255)
                                                                                                                                            cACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAGCGCAAAc        >  2:1033930/1‑149 (MQ=255)
                                                                                                                                             aCGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAGCGCAAACt       <  2:12608/149‑1 (MQ=255)
                                                                                                                                              cgcgGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAGCGCAAACtt      <  2:4653269/149‑1 (MQ=255)
                                                                                                                                                cgGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAGCGCAAACTTcc    >  2:4327870/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCATG  >  CP009273/3557668‑3557962

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