Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,574,820 C→A 100% A609E (GCG→GAG)  tktB → transketolase 2, thiamine triphosphate‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,574,8200CA99.6% 729.2 / ‑5.7 225A609E (GCG→GAG) tktBtransketolase 2, thiamine triphosphate‑binding
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (113/111);  total (114/111)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 225 total aligned reads displayed.

GAAATTACCCTGCAAGCGGCAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGG  >  CP009273/2574675‑2574960
                                                                                                                                                 |                                                                                                                                            
gAAATTACCCTGCAAGCGGCAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAgg                                                                                                                                            <  2:489809/148‑1 (MQ=255)
 aaaTTACCCTGCAAGCGGCAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGt                                                                                                                                           >  1:1644134/1‑148 (MQ=255)
   aTTACCCTGCAAGCGGCAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGACTATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTgg                                                                                                                                         >  2:1951126/1‑148 (MQ=255)
   aTTACCCTGCAAGCGGCAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTgg                                                                                                                                         >  2:3500810/1‑148 (MQ=255)
     tACCCTGCAAGCGGCAGAGAAATTAGCAGGAGAAGGTCGCACTGTACGCGTAGTTTCCCTGCCCTCGACCGATACTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGaa                                                                                                                                       >  2:1501032/1‑148 (MQ=255)
     tACCCTGCAAGCGGCAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGaa                                                                                                                                       >  1:1855827/1‑148 (MQ=255)
     tACCCTGCAAGCGGCAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGaa                                                                                                                                       >  2:3496606/1‑148 (MQ=255)
     tACCCTGCAAGCGGCAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGGGGAGGTGGaa                                                                                                                                       >  2:2983736/1‑148 (MQ=255)
       cccTGCAAGCGGCAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGc                                                                                                                                     <  1:3342320/148‑1 (MQ=255)
           gCAAGCGGCAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGt                                                                                                                                 >  2:1162241/1‑148 (MQ=255)
             aaGCGGCAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTAt                                                                                                                               <  1:1957583/148‑1 (MQ=255)
               gCGGCAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTg                                                                                                                             >  2:615347/1‑148 (MQ=255)
               gCGGCAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTg                                                                                                                             <  1:2900750/148‑1 (MQ=255)
               gCGGCAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTg                                                                                                                             <  2:3062700/148‑1 (MQ=255)
                 ggCAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGcc                                                                                                                           <  2:1569736/148‑1 (MQ=255)
                  gCAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCg                                                                                                                          <  2:2336365/148‑1 (MQ=255)
                  gCAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCg                                                                                                                          <  2:2336594/148‑1 (MQ=255)
                    agagAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGAt                                                                                                                        <  2:3796235/148‑1 (MQ=255)
                     gagaAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGAtt                                                                                                                       >  1:612695/1‑148 (MQ=255)
                        aaaTTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACt                                                                                                                    <  1:972665/148‑1 (MQ=255)
                         aaTTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTg                                                                                                                   <  2:2374861/148‑1 (MQ=255)
                          aTTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCGGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTgg                                                                                                                  >  1:458154/1‑148 (MQ=255)
                          aTTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTgg                                                                                                                  <  1:2512780/148‑1 (MQ=255)
                           ttAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGt                                                                                                                 >  1:2737171/1‑148 (MQ=255)
                            tAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTa                                                                                                                >  2:1635219/1‑148 (MQ=255)
                                 ggAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGt                                                                                                           >  1:2222660/1‑148 (MQ=255)
                                  gagaAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTa                                                                                                          <  2:410222/148‑1 (MQ=255)
                                     aaGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGt                                                                                                       >  1:2117208/1‑148 (MQ=255)
                                       ggTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTg                                                                                                     >  1:3650949/1‑148 (MQ=255)
                                        gTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTgg                                                                                                    >  1:1583076/1‑148 (MQ=255)
                                          cGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTc                                                                                                  <  1:1151446/148‑1 (MQ=255)
                                          cGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTc                                                                                                  <  1:2087331/148‑1 (MQ=255)
                                            cAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTg                                                                                                >  2:2827308/1‑148 (MQ=255)
                                              aTGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGaa                                                                                              <  1:771600/148‑1 (MQ=255)
                                               tGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGaaa                                                                                             <  1:900083/148‑1 (MQ=255)
                                               tGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGaaa                                                                                             <  1:738800/148‑1 (MQ=255)
                                                gTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAg                                                                                            >  2:235516/1‑148 (MQ=255)
                                                 tACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAgg                                                                                           >  2:1133485/1‑148 (MQ=255)
                                                  aCGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAggg                                                                                          >  1:26103/1‑148 (MQ=255)
                                                     cgTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGAATGTTGGTCTGAAAGGGGCa                                                                                       >  1:1559667/1‑148 (MQ=255)
                                                      gTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCaa                                                                                      <  1:2302735/148‑1 (MQ=255)
                                                       tAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAAt                                                                                     <  2:3086337/148‑1 (MQ=255)
                                                       tAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTAGGTTGGTCTGAAAGGGGCAAt                                                                                     >  2:591003/1‑148 (MQ=255)
                                                         gTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTg                                                                                   >  2:2160270/1‑148 (MQ=255)
                                                          tttCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGt                                                                                  >  1:444489/1‑148 (MQ=255)
                                                              ccTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCggg                                                                              <  2:1855827/148‑1 (MQ=255)
                                                                  cccTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGa                                                                          <  1:3503012/148‑1 (MQ=255)
                                                                  cccTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGa                                                                          <  2:2555030/148‑1 (MQ=255)
                                                                    cTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACg                                                                        <  1:560013/148‑1 (MQ=255)
                                                                      cGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACggg                                                                      >  2:2660333/1‑148 (MQ=255)
                                                                      cGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACggg                                                                      <  1:2323831/148‑1 (MQ=255)
                                                                           gATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACg                                                                 <  1:110334/148‑1 (MQ=255)
                                                                            atatTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACgg                                                                <  2:1077800/148‑1 (MQ=255)
                                                                             tatTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATAGTCGGGATGACGGGTTACggg                                                               >  2:3113083/1‑148 (MQ=255)
                                                                              atTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACgggg                                                              >  1:2095140/1‑148 (MQ=255)
                                                                               ttttCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGa                                                             >  2:1728920/1‑148 (MQ=255)
                                                                                 ttCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAAt                                                           <  1:3280513/148‑1 (MQ=255)
                                                                                   cGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCt                                                         >  2:859241/1‑148 (MQ=255)
                                                                                     aCGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGc                                                       >  2:3494711/1‑148 (MQ=255)
                                                                                     aCGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGc                                                       >  2:3207475/1‑148 (MQ=255)
                                                                                     aCGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGc                                                       >  2:3361174/1‑148 (MQ=255)
                                                                                      cGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCt                                                      <  2:3383908/148‑1 (MQ=255)
                                                                                          cAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCcgg                                                  >  1:2303351/1‑148 (MQ=255)
                                                                                           aGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCcggc                                                 >  1:3305682/1‑148 (MQ=255)
                                                                                            ggATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCcggcg                                                <  1:2798745/148‑1 (MQ=255)
                                                                                                  ggAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAg                                          >  2:2892275/1‑148 (MQ=255)
                                                                                                    aaTATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCt                                        <  2:1798198/148‑1 (MQ=255)
                                                                                                     atatCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTg                                       <  2:2307391/148‑1 (MQ=255)
                                                                                                      tatCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGt                                      >  1:2547451/1‑148 (MQ=255)
                                                                                                       atCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTAGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGtt                                     <  1:2559982/148‑1 (MQ=255)
                                                                                                        tCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTc                                    <  1:1887799/148‑1 (MQ=255)
                                                                                                        tCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTc                                    <  1:1887634/148‑1 (MQ=255)
                                                                                                        tCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTc                                    >  2:1706562/1‑148 (MQ=255)
                                                                                                           ggAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCg                                 >  2:833735/1‑148 (MQ=255)
                                                                                                              gTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTc                              <  2:299118/148‑1 (MQ=255)
                                                                                                               tCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCt                             >  2:373395/1‑148 (MQ=255)
                                                                                                                   tgtTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTgg                         >  2:1786065/1‑148 (MQ=255)
                                                                                                                     ttGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCt                       >  2:639992/1‑148 (MQ=255)
                                                                                                                          ttCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTAcc                  >  1:1534025/1‑148 (MQ=255)
                                                                                                                          ttCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTAcc                  >  2:2058854/1‑148 (MQ=255)
                                                                                                                          ttCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTAcc                  >  1:3433024/1‑148 (MQ=255)
                                                                                                                           tCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTAccg                 >  2:2709762/1‑148 (MQ=255)
                                                                                                                           tCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTAccg                 >  2:2709563/1‑148 (MQ=255)
                                                                                                                              aaCGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTAccgccg              <  2:2899295/148‑1 (MQ=255)
                                                                                                                              aaCGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTAccgccg              >  2:2100574/1‑148 (MQ=255)
                                                                                                                               aCGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCga             >  2:3247911/1‑148 (MQ=255)
                                                                                                                                 gTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCgaga           <  1:956633/148‑1 (MQ=255)
                                                                                                                                 gTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCgaga           <  2:1288936/148‑1 (MQ=255)
                                                                                                                                   tGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAat         >  2:2640510/1‑148 (MQ=255)
                                                                                                                                   tGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAat         <  2:3555139/148‑1 (MQ=255)
                                                                                                                                   tGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAat         <  1:18044/148‑1 (MQ=255)
                                                                                                                                   tGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAat         <  2:3152077/148‑1 (MQ=255)
                                                                                                                                    gCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAata        <  1:3013412/148‑1 (MQ=255)
                                                                                                                                     cGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAatat       <  1:1290116/148‑1 (MQ=255)
                                                                                                                                       gCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATtg     <  2:1836571/148‑1 (MQ=255)
                                                                                                                                       gCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATtg     <  2:913762/148‑1 (MQ=255)
                                                                                                                                       gCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATtg     <  2:747051/148‑1 (MQ=255)
                                                                                                                                       gCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATtg     <  2:3727549/148‑1 (MQ=255)
                                                                                                                                        cTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATtgt    <  1:1021386/148‑1 (MQ=255)
                                                                                                                                          cgcgTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTgg  <  1:681842/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                            
GAAATTACCCTGCAAGCGGCAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGG  >  CP009273/2574675‑2574960

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: