Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,559,170 C→A 7.0% 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 7.0% 473.7 / 2.6 171D258Y (GAT→TAT) glgPglycogen phosphorylase
Reads supporting (aligned to +/- strand):  ref base C (59/90);  new base A (17/2);  total (79/92)
Fisher's exact test for biased strand distribution p-value = 3.57e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 175 total aligned reads displayed.

TTTTATGCAACTGATAATGGCGGCTTAAAATGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTC  >  CP009273/3559023‑3559317
                                                                                                                                                   |                                                                                                                                                   
ttttATGCAACTGATAATGGCGGCTTAAAATGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATc                                                                                                                                                     <  2:2843201/148‑1 (MQ=255)
 tttATGCAACTGATAATGGCGGCTTAAAATGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCt                                                                                                                                                    <  2:2436605/148‑1 (MQ=255)
 tttATGCAACTGATAATGGCGGCTTAAAATGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCt                                                                                                                                                    <  1:2359654/148‑1 (MQ=255)
   tATGCAACTGATAATGGCGGCTTAAAATGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTc                                                                                                                                                  <  1:2197916/148‑1 (MQ=255)
   tATGCAACTGATAATGGCGGCTTAAAATGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTc                                                                                                                                                  <  2:602674/148‑1 (MQ=255)
      gCAACTGATAATGGCGGCTTAAAATGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATATTCCaa                                                                                                                                               >  1:2443560/1‑147 (MQ=255)
              tAATGGCGGCTTAAAATGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTTTGTTACACAGCcgcgc                                                                                                                                       >  1:1466087/1‑147 (MQ=255)
              tAATGGCGGCTTAAAATGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTGTTTATATTACACTTCCGCGa                                                                                                                                       >  1:1466599/1‑148 (MQ=255)
              tAATGGCGGCTTAAAATGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTGTTTATATTAAAATGCcgcgc                                                                                                                                       >  1:1993649/1‑147 (MQ=255)
               aaTGGCGGCTTAAAATGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTGTTTATGTTCCACTGCcccgca                                                                                                                                      >  1:1218589/1‑146 (MQ=255)
                     ggcTTAAAATGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTc                                                                                                                                <  1:1945995/148‑1 (MQ=255)
                         tAAAATGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCAccc                                                                                                                            <  2:711490/148‑1 (MQ=255)
                         tAAAATGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCAccc                                                                                                                            <  1:661817/148‑1 (MQ=255)
                         tAAAATGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCAccc                                                                                                                            <  1:2826851/148‑1 (MQ=255)
                          aaaaTGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCt                                                                                                                           <  1:547986/148‑1 (MQ=255)
                           aaaTGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTg                                                                                                                          <  1:2356230/148‑1 (MQ=255)
                           aaaTGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATATTACACTGCCGCGAAGTAGTCAcccag                                                                                                                          >  1:2654254/1‑146 (MQ=255)
                             aTGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGt                                                                                                                        <  1:2507835/148‑1 (MQ=255)
                               gTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTg                                                                                                                      >  2:3158525/1‑148 (MQ=255)
                                 ccTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTCGaa                                                                                                                    <  2:1806829/148‑1 (MQ=255)
                                   tGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAAtt                                                                                                                  >  2:2839876/1‑148 (MQ=255)
                                   tGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAAtt                                                                                                                  <  2:1160515/148‑1 (MQ=255)
                                   tGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTGTTTATCTTGCACTGCCGCGCAGTAGTAACCCGGGttttatt                                                                                                                  >  1:2018806/1‑143 (MQ=255)
                                   tGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTGTTTATATTCCACTGCCGCGAAGTAGTCACCCGGGTTGAAtt                                                                                                                  >  1:354117/1‑148 (MQ=255)
                                     aaTGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCAAGCCCGGAGTAGCTGGTGAGATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACAGCCGCGAAGTAGTCACCCTGGTTGAATTTa                                                                                                                >  1:3249581/1‑148 (MQ=255)
                                      aTGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTAc                                                                                                               >  2:1558612/1‑148 (MQ=255)
                                        ggTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACANTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCg                                                                                                             <  1:1681474/148‑1 (MQ=255)
                                        ggTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCg                                                                                                             <  1:2626341/148‑1 (MQ=255)
                                        ggTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCg                                                                                                             <  1:1452285/148‑1 (MQ=255)
                                         gTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGa                                                                                                            >  2:2349795/1‑148 (MQ=255)
                                         gTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGa                                                                                                            >  2:2350226/1‑148 (MQ=255)
                                         gTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGa                                                                                                            >  2:412597/1‑148 (MQ=255)
                                         gTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGa                                                                                                            <  1:2679112/148‑1 (MQ=255)
                                         gTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGa                                                                                                            >  2:819877/1‑148 (MQ=255)
                                         gTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATATTCCACAGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGa                                                                                                            >  1:1045181/1‑148 (MQ=255)
                                           cGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAgg                                                                                                          <  2:440023/148‑1 (MQ=255)
                                           cGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAgg                                                                                                          <  1:2851228/148‑1 (MQ=255)
                                             aGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGtt                                                                                                        <  2:1443469/148‑1 (MQ=255)
                                             aGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTGTTTATCTTCCACAGCCGCGCAGTATTCACCCTGGTTGAATTTACCGAGGtt                                                                                                        >  1:2699344/1‑148 (MQ=255)
                                              ggAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACAGCCGCGCAGGAGGCACCCGGGGTGAATTTACCGAGGTTa                                                                                                       >  1:1056545/1‑148 (MQ=255)
                                               gAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTaa                                                                                                      <  2:871772/148‑1 (MQ=255)
                                                aaaCCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAAt                                                                                                     <  1:2519704/148‑1 (MQ=255)
                                                aaaCCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTGTTTATCTTCCACAGCCGCGACGAAGTCACCCGGGTTGAATTTACCGAGGtttag                                                                                                     >  1:2068038/1‑145 (MQ=255)
                                                 aaCCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAAtt                                                                                                    <  2:643316/148‑1 (MQ=255)
                                                  aCCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTTTATTACACAGCCGCGAAGTAGTCACCCGGGTTGAATTTAACGAGGGTTAttt                                                                                                   >  1:730843/1‑148 (MQ=255)
                                                  aCCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTGTTTTTCTTCCACAGCCGCGAAGTAGTCAACCGGGTTGAATTTACCGAGGGTAAttt                                                                                                   >  1:3124176/1‑148 (MQ=255)
                                                    cAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTACACTGCCGCGCAGTAGTCACCCAGGTTGAATTTAACGAGGTTAATTTCg                                                                                                 >  1:1840030/1‑148 (MQ=255)
                                                      ggAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCt                                                                                               >  2:2578026/1‑148 (MQ=255)
                                                       gAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATATTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTa                                                                                              <  1:722512/148‑1 (MQ=255)
                                                        aaGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTAc                                                                                             >  2:772345/1‑148 (MQ=255)
                                                        aaGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTAc                                                                                             >  2:772438/1‑148 (MQ=255)
                                                            aTTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGc                                                                                         <  2:3118182/148‑1 (MQ=255)
                                                             ttCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCt                                                                                        <  2:1355057/148‑1 (MQ=255)
                                                               ccTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTg                                                                                      <  2:1364807/148‑1 (MQ=255)
                                                                   aCGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCa                                                                                  <  1:3105810/148‑1 (MQ=255)
                                                                        ggCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATGTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTAACGACGTTAATTTAGCTACTTGCTTGCGCACTCca                                                                             >  1:2617014/1‑148 (MQ=255)
                                                                        ggCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACAGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTAGCTACTGGCTTGCGCACACca                                                                             >  1:3034537/1‑148 (MQ=255)
                                                                              gCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACTTTCTCGGAGTGGTTTTTTTATTACACAGCCCCGAACTAGTCACCCTGGTTGGATTTACCGAGGGTAATTTCGCTACATGCTTGCGCCCACCACcaaag                                                                       >  1:1259289/1‑146 (MQ=255)
                                                                              gCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTGTTTATCTTCCACTGCCGCGAAGTAGTCACCCGGGTTGAATTTACCGAGGTTTATTTCGCTACGGGCTTGCGCACTCCACAAACg                                                                       >  2:3092683/1‑148 (MQ=255)
                                                                               cTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCAATCCACAAACGc                                                                      >  2:1918194/1‑148 (MQ=255)
                                                                               cTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGGAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCAATCCACAAACGc                                                                      >  2:1918135/1‑148 (MQ=255)
                                                                               cTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTGTTTATCTTCCACTGCCGCGAACTAGTCACCCGGGTTGAATTTACCGAGGTTAATTTACCTACATGCATGCGCAATACAAAAACGc                                                                      >  2:2667509/1‑148 (MQ=255)
                                                                               cTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGAGACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGc                                                                      <  2:2705763/148‑1 (MQ=255)
                                                                                  aCGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTTTCTTGCACAGCCGCCAAGGAGTAACCCGGGTTGAATTTTACGAGGGTAATTTGGCTAAATGGTTGCGCACTCCACAAACCCAGc                                                                   >  1:999588/1‑148 (MQ=255)
                                                                                  aCGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGc                                                                   >  2:2147260/1‑148 (MQ=255)
                                                                                  aCGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCAGGTTGAAATTACCGAGGTTAATTTCGCTACTGGCTTGCGCCCTCCACCAAagcagc                                                                   >  1:2204516/1‑148 (MQ=255)
                                                                                  aCGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTGTTTATCTTGCACAGCCCCGAAGGAGTCACCCGGGGTGAATTTACAGAGGGTAATTTCGCTACAGGCTTGCGCCCACCACcaaaccacc                                                                   >  1:999848/1‑142 (MQ=255)
                                                                                     cccGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTGTTTTTATTCCACAGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTAGCTACTGGCTTTCGCACGCCACAAACGCAGCggg                                                                >  1:1301009/1‑146 (MQ=255)
                                                                                       cGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTTTATTACACAACCCCGAACTAGTCACCCGGGTTTTAATTAACCAGGTTAATTTAGCGACTGGCTTTGGCACTCCACACAACCAGCGGGtt                                                              >  1:3063160/1‑148 (MQ=255)
                                                                                        ggAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTg                                                             <  1:432579/148‑1 (MQ=255)
                                                                                              ggTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACAGCCGCGAAGTAGTCACCCAGGTTGAATTTACCGAGGTTAAATTAGCTACTGGCTTGCGCACTCCACAAAAGCAGCGTGTTGGTcgcg                                                       >  1:1354592/1‑148 (MQ=255)
                                                                                                tGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTACACAGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTc                                                     >  2:2050673/1‑148 (MQ=255)
                                                                                                    gTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGgtg                                                 <  1:2349795/148‑1 (MQ=255)
                                                                                                    gTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTACACAGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTAGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTTGTCGCGTCGgtg                                                 >  2:1776564/1‑148 (MQ=255)
                                                                                                        tCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGCAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTACACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGt                                             >  2:2537157/1‑148 (MQ=255)
                                                                                                         ccGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTa                                            <  1:760697/148‑1 (MQ=255)
                                                                                                          cGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTaa                                           <  2:1993887/148‑1 (MQ=255)
                                                                                                            gATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATATTACACAGCCGCGAAGTAGTCACCCGGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGGGTTGGTCGCGTCGGTGTCGTAAcc                                         >  2:2862994/1‑148 (MQ=255)
                                                                                                                cAGTACGCGAGATTCGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTGAATTTCTCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAggg                                     <  1:1755439/148‑1 (MQ=255)
                                                                                                                cAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAggg                                     <  2:1731172/148‑1 (MQ=255)
                                                                                                                cAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAggg                                     <  2:1729681/148‑1 (MQ=255)
                                                                                                                 aGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGa                                    <  1:2478146/148‑1 (MQ=255)
                                                                                                                 aGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGa                                    <  2:1944926/148‑1 (MQ=255)
                                                                                                                 aGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGa                                    <  1:1558612/148‑1 (MQ=255)
                                                                                                                        gagaTACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTg                             <  1:2578026/148‑1 (MQ=255)
                                                                                                                         agaTACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGa                            <  2:56261/148‑1 (MQ=255)
                                                                                                                         agaTACGTTCTCGGAGTGGTTTTTATCTTCCACAGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGGGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGa                            >  2:502598/1‑148 (MQ=255)
                                                                                                                          gaTACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGAt                           <  1:486627/148‑1 (MQ=255)
                                                                                                                             aCGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGt                        <  1:3119362/148‑1 (MQ=255)
                                                                                                                              cGTTCTCGGAGTGGTTTTTATATTCCACAGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGGGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTa                       >  2:2883460/1‑148 (MQ=255)
                                                                                                                               gTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTaa                      <  1:772438/148‑1 (MQ=255)
                                                                                                                                ttCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAg                     <  1:1320134/148‑1 (MQ=255)
                                                                                                                                    cGGTGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTAACCCTGGTTTAATTTACCGATGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGGCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGAc                 <  2:2108887/148‑1 (MQ=255)
                                                                                                                                      gAGTGGTTTTTATCTTCCACAGCCGCGAAGTAGTCACCCTGGTTGAATTTAACGAGGTTAATTTCGCTACTGGCTTGCGCACACCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAACCGACTc               >  2:1048129/1‑148 (MQ=255)
                                                                                                                                      gAGTGGTGTTTATCTTCCACAGCCGCGAAGTAGTCACCCTGGTTGAATTTAACCAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGGGTTGGGCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTc               >  2:2400303/1‑148 (MQ=255)
                                                                                                                                              tttATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTTAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCAAAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAAtc       <  2:3239518/148‑1 (MQ=255)
                                                                                                                                               ttATCTTCCACGGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAAtct      <  2:1983271/148‑1 (MQ=255)
                                                                                                                                                tATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAAtctc     >  1:3002054/1‑148 (MQ=255)
                                                                                                                                                  tCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCtt   >  1:2003229/1‑148 (MQ=255)
                                                                                                                                                   cTTCAACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACAGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGACTCTCTTc  >  2:3043071/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TTTTATGCAACTGATAATGGCGGCTTAAAATGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTC  >  CP009273/3559023‑3559317

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: