Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,559,170 C→A 8.4% D258Y (GAT→TAT)  glgP ← glycogen phosphorylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,559,1700CA 8.4% 768.0 / 5.9 263D258Y (GAT→TAT) glgPglycogen phosphorylase
Reads supporting (aligned to +/- strand):  ref base C (77/147);  new base A (35/2);  total (114/149)
Fisher's exact test for biased strand distribution p-value = 1.15e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 265 total aligned reads displayed.

GTTTTATGCAACTGATAATGGCGGCTTAAAATGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGG  >  CP009273/3559022‑3559319
                                                                                                                                                    |                                                                                                                                                     
 ttttATGCAACTGATAATGGCGGCTTAAAATGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCtt                                                                                                                                                     <  1:46780/150‑1 (MQ=255)
    tATGCAACTGATAATGGCGGCTTAAAATGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCa                                                                                                                                                  <  2:579702/150‑1 (MQ=255)
     aTGCAACTGATAATGGCGGCTTAAAATGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCAc                                                                                                                                                 <  2:2567349/150‑1 (MQ=255)
         aaCTGATAATGGCGGCTTAAAATGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCATGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATATTCCACAGcc                                                                                                                                             >  1:284373/1‑150 (MQ=255)
         aaCTGATAATGGCGGCTTAAAATGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGcc                                                                                                                                             <  1:3729586/150‑1 (MQ=255)
         aaCTGATAATGGCGGCTTAAAATGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGcc                                                                                                                                             <  2:3339503/150‑1 (MQ=255)
           cTGATAATGGCGGCTTAAAATGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATATTACACAGCcgc                                                                                                                                           >  1:1899696/1‑150 (MQ=255)
            tGATAATGGCGGCTTAAAATGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCcgcg                                                                                                                                          <  1:782026/150‑1 (MQ=255)
            tGATAATGGCGGCTTAAAATGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATATTACACTGCcgcg                                                                                                                                          >  1:2982270/1‑150 (MQ=255)
             gATAATGGCGGCTTAAAATGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGa                                                                                                                                         >  1:3401719/1‑150 (MQ=255)
                  tGGCGGCTTAAAATGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAagtag                                                                                                                                    <  1:2778137/150‑1 (MQ=255)
                   ggcggcTTAAAATGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAagtagt                                                                                                                                   <  2:1813373/150‑1 (MQ=255)
                   ggcggcTTAAAATGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAagtagt                                                                                                                                   <  2:3314422/150‑1 (MQ=255)
                        cTTAAAATGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTGTTTATCTTACACAGCCCCGAAGTATTAAccc                                                                                                                              >  1:438228/1‑150 (MQ=255)
                         ttACAATGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCt                                                                                                                             <  1:2484628/150‑1 (MQ=255)
                          tAAAATGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTg                                                                                                                            <  2:1461406/150‑1 (MQ=255)
                          tAAAATGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTACACTGCCGCGAAGTAGTCAcccgg                                                                                                                            >  2:3235704/1‑148 (MQ=255)
                           aaaaTGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTgg                                                                                                                           <  2:115830/150‑1 (MQ=255)
                           aaaaTGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACAGCCGCGAAGTAGTCACCCggg                                                                                                                           >  1:3515823/1‑150 (MQ=255)
                            aaaTGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATATTCCACTTCCGCGAAGTAGTCACCCAGGt                                                                                                                          >  1:3020752/1‑150 (MQ=255)
                             aaTGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGtt                                                                                                                         <  2:2692483/150‑1 (MQ=255)
                              aTGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTGTTTATATTACACAGCCGCGCAGTAGTAACCCAGGTTg                                                                                                                        >  1:3474969/1‑150 (MQ=255)
                                gTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGaa                                                                                                                      >  2:3591474/1‑150 (MQ=255)
                                gTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTGTTTATATTCCACAGCCGCGAAGTAGTCACCCTGGTTGaa                                                                                                                      >  1:1636499/1‑150 (MQ=255)
                                 tCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAAt                                                                                                                     <  2:4295362/150‑1 (MQ=255)
                                 tCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAAt                                                                                                                     <  1:2423790/150‑1 (MQ=255)
                                  ccTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAAtt                                                                                                                    >  1:269434/1‑150 (MQ=255)
                                    tGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTGTTTATATTACACAGCCGCGAAGTATTACACCAGGTTGAATTTa                                                                                                                  >  1:1750560/1‑150 (MQ=255)
                                     gAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTAc                                                                                                                 <  1:4757235/150‑1 (MQ=255)
                                      aaTGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCAGGTTGAATTTAcc                                                                                                                >  1:747184/1‑150 (MQ=255)
                                        tGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTGTTTATCTTCCACAGCCGCGAAGTAGTCACCCGGGTTGAATTTAACGa                                                                                                              >  1:1093341/1‑150 (MQ=255)
                                            cGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGtt                                                                                                          <  1:3337969/150‑1 (MQ=255)
                                            cGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGtt                                                                                                          <  1:1021662/150‑1 (MQ=255)
                                            cGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGtt                                                                                                          <  1:212611/150‑1 (MQ=255)
                                              aGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTACACTGCCGCGAAGTAGTCACCCTGGTTTAATTTTACGAGGTTaa                                                                                                        >  2:2571019/1‑150 (MQ=255)
                                               ggAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGCAGTAGTCACCCAGGTTGAATTTACCGAGGTTCAt                                                                                                       >  2:3694410/1‑150 (MQ=255)
                                                 aaaCCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAAttt                                                                                                     <  1:1486922/150‑1 (MQ=255)
                                                 aaaCCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATATTCCACAGCCGCGAAGTAGTCACCCGGGTTGAATTTACCGAGGTTAAttt                                                                                                     >  2:505320/1‑150 (MQ=255)
                                                  aaCCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTc                                                                                                    <  2:1068331/150‑1 (MQ=255)
                                                  aaCCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTc                                                                                                    <  2:4704043/150‑1 (MQ=255)
                                                  aaCCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTc                                                                                                    <  2:3471621/150‑1 (MQ=255)
                                                   aCCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCg                                                                                                   <  1:4519336/150‑1 (MQ=255)
                                                     cAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCt                                                                                                 <  1:3318149/150‑1 (MQ=255)
                                                       ggAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTAc                                                                                               <  1:1155876/150‑1 (MQ=255)
                                                       ggAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTAc                                                                                               <  2:722880/150‑1 (MQ=255)
                                                          aGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTgg                                                                                            <  1:4579555/150‑1 (MQ=255)
                                                           gTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTTTATTACACAGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTAGCTaaaggc                                                                                           >  2:2837545/1‑145 (MQ=255)
                                                            tATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCt                                                                                          <  1:311884/150‑1 (MQ=255)
                                                              ttCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTg                                                                                        >  2:3089719/1‑150 (MQ=255)
                                                               tCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTgc                                                                                       <  2:739414/150‑1 (MQ=255)
                                                                ccTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTgcg                                                                                      <  1:3467213/150‑1 (MQ=255)
                                                                ccTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGAGACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTgcg                                                                                      <  2:4580763/150‑1 (MQ=255)
                                                                       cAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCca                                                                               <  1:1988399/150‑1 (MQ=255)
                                                                       cAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTGTTTATATTCCCCAGCCCCGCAGTAGTCACCCAGGGTGAATTTAACGAGGTTAATTTCGCTACTGGCTTGCGCACTCca                                                                               >  1:4233224/1‑150 (MQ=255)
                                                                            gcAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACAGCCGCGAAGGAGTCACCCGGGTTGAATTTACCGAGGTTAATTTCGCTACAGGCTTGCGCACTCCACCAAc                                                                          >  2:423226/1‑150 (MQ=255)
                                                                             cAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTACACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACg                                                                         >  2:3945642/1‑150 (MQ=255)
                                                                              aGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATATTCCACAGCCGCGAAGTAGTCACCCGGGTTGAATTTAACGAGGTTAATTTAGCTACTGGCTTGCGCACTCCACAAACGc                                                                        >  2:907144/1‑150 (MQ=255)
                                                                                 tCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGc                                                                     <  1:985864/150‑1 (MQ=255)
                                                                                 tCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTGTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTGATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGc                                                                     >  2:4143789/1‑150 (MQ=255)
                                                                                     gCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCgtgt                                                                 <  1:266714/150‑1 (MQ=255)
                                                                                       ccGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTg                                                               <  2:1720338/150‑1 (MQ=255)
                                                                                         ggAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATAATACACAGCCGCGAAGTAGTAACCCTGGTTGAATTTACCGAGGTTAATTTAGCTACTGGCTTGCGCACTCCACAAAAGCAGCGGGTTggg                                                             >  2:529951/1‑149 (MQ=255)
                                                                                          gAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTc                                                            <  2:3416570/150‑1 (MQ=255)
                                                                                           aGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTGTTTATCTTCCACAGCCGCGAAGTAGTCACCCAGGTTGAATTTACCGAGGTTAATTTCGCTACAGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTcg                                                           >  2:3788291/1‑150 (MQ=255)
                                                                                               ggTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTc                                                       <  1:2237925/150‑1 (MQ=255)
                                                                                                  ggAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGgt                                                    <  2:1899696/150‑1 (MQ=255)
                                                                                                  ggAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGgt                                                    <  1:3796523/150‑1 (MQ=255)
                                                                                                  ggAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGgt                                                    >  2:4713876/1‑150 (MQ=255)
                                                                                                  ggAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACAGCCGCGAAGTAGTCACCCGGGTTGAATTTACCGAGGGTAATTTCGCTACTGGCTTGCGCACACCACAAACGCAGCGTGTTGGTCGCGTCGgt                                                    >  2:1625099/1‑150 (MQ=255)
                                                                                                  ggAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTACACAGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGgt                                                    >  2:4517926/1‑150 (MQ=255)
                                                                                                      tCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCg                                                <  2:3278338/150‑1 (MQ=255)
                                                                                                       cATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGt                                               <  1:1058722/150‑1 (MQ=255)
                                                                                                        aTCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATATTACACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCGTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTa                                              >  2:4672677/1‑150 (MQ=255)
                                                                                                         tCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACAGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTaa                                             >  2:4377623/1‑150 (MQ=255)
                                                                                                         tCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCATCCACAGCCCCGAAGTAGTCACCCGGGTTGAATTTACCGAGGTTAATTTGGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTaa                                             >  1:1964884/1‑150 (MQ=255)
                                                                                                            ggATACAGTACGCGAGATACGTTCTCGGAGTGGTGTTTATCTTCCACAGCCGCGAAGTAGTCACCCGGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCCCAAACGCAGCGTGGTGGTCGCGGCGGTGTCGTAACCa                                          >  2:1504899/1‑150 (MQ=255)
                                                                                                              aTACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAgg                                        <  1:3887904/150‑1 (MQ=255)
                                                                                                                aCAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGa                                      <  2:1401117/150‑1 (MQ=255)
                                                                                                                 cAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGAt                                     <  1:419219/150‑1 (MQ=255)
                                                                                                                 cAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCGTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTGAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGAt                                     <  1:498175/150‑1 (MQ=255)
                                                                                                                  aGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGAtt                                    <  2:2441265/150‑1 (MQ=255)
                                                                                                                  aGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGAtt                                    <  2:2442286/150‑1 (MQ=255)
                                                                                                                  aGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGGGTTGGTCGCGTCGGTGTCGTAACCCGGGAtt                                    >  2:3542618/1‑150 (MQ=255)
                                                                                                                    tACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTAt                                  <  1:2808940/150‑1 (MQ=255)
                                                                                                                    tACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCAGGTTGAATTTAACGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTAt                                  >  2:2257112/1‑150 (MQ=255)
                                                                                                                     aCGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATc                                 <  1:2747251/150‑1 (MQ=255)
                                                                                                                      cgcgAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCt                                <  2:618212/150‑1 (MQ=255)
                                                                                                                            aTACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGt                          <  1:3627888/150‑1 (MQ=255)
                                                                                                                             tACGTTCTCGGAGTGGTTTTTATCTTCCACAGCCGCGAAGTAGTCACCCGGGGTGAATTTACCGACGTTAATTTCGCTCCTGGCTTGCGCACTCCACAAACGCACCGGGGTGGTCGCGGCGGTGTCGTAACCCGGGAATATCTGATCGTa                         >  2:3511412/1‑150 (MQ=255)
                                                                                                                               cGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAg                       <  2:3474969/150‑1 (MQ=255)
                                                                                                                                gTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGc                      <  1:3085598/150‑1 (MQ=255)
                                                                                                                                 ttCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCg                     <  1:838173/150‑1 (MQ=255)
                                                                                                                                 ttCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCg                     <  1:4314942/150‑1 (MQ=255)
                                                                                                                                   ctcGGAGTGGTTTTTATCTTCCACAGCCGCGAAGTATTCACCCTGGTTGAATTTAACCAGGTTAATTTCGCGACTGGCTTGCGCACACCACAAACGCAGCGGGGTGGTCGCGTCGGTGTCGTAACCAGGGCTTATATGATCGTAAGCGAc                   >  2:2580016/1‑150 (MQ=255)
                                                                                                                                         gtgtTTTTTATCTTCCACTGCCGCTAAGTAGTCACCCTGTTTGATTTTACCGTGGTTATTTTCGCTAGTGGCTTGCGCACTCAACAAACGCAGGGTGTGGGTCGCGTGGGTGTCGTAACCAGGTTTTATCTGATCGTAAGCGACTCCCAg             <  2:3846873/150‑1 (MQ=255)
                                                                                                                                               tttATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAAtctc       >  1:455684/1‑150 (MQ=255)
                                                                                                                                                ttATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATctct      >  2:2971906/1‑150 (MQ=255)
                                                                                                                                                ttATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATctct      >  2:1349996/1‑150 (MQ=255)
                                                                                                                                                   tCTTCCACTGCCGCGATTTAGTCACCCTGGTTGATTTTACCGTGTTTAATTTCGCTACTGGCTTGCGCACCCCACAAACGCTGCGTTTTGGTCGCGTCGGTGTCTTAACCAGGGATTATCTTATCGTAAGCGACTCCCAGAATCTCTTCg   <  2:2875403/150‑1 (MQ=255)
                                                                                                                                                   tCTTCCACTGCCGCGAAGTAGTCACCCTGGTTTAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCg   <  2:1954362/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GTTTTATGCAACTGATAATGGCGGCTTAAAATGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGG  >  CP009273/3559022‑3559319

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: