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,8200CA97.7% 536.1 / 3.5 174A609E (GCG→GAG) tktBtransketolase 2, thiamine triphosphate‑binding
Reads supporting (aligned to +/- strand):  ref base C (4/0);  new base A (85/84);  total (89/85)
Fisher's exact test for biased strand distribution p-value = 1.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01
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 174 total aligned reads displayed.

GAAATTACCCTGCAAGCGGCAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGGCA  >  CP009273/2574675‑2574962
                                                                                                                                                 |                                                                                                                                              
gAAATTACCCTGCAAGCGGCAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAgg                                                                                                                                              >  2:3409834/1‑148 (MQ=255)
 aaaTTACCCTGCAAGCGGCAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGt                                                                                                                                             >  2:233568/1‑148 (MQ=255)
 aaaTTACCCTGCAAGCGGCAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGt                                                                                                                                             >  1:2047915/1‑148 (MQ=255)
   aTTACCCTGCAAGCGGCAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTgg                                                                                                                                           >  1:2900351/1‑148 (MQ=255)
    gtACCCTGCAAGCGGCAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGa                                                                                                                                          <  2:2419958/147‑1 (MQ=255)
         cTGCAAGCGGCAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAg                                                                                                                                     <  1:670614/148‑1 (MQ=255)
         cTGCAAGCGGCAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAg                                                                                                                                     <  1:670689/148‑1 (MQ=255)
         cTGCAAGCGGCAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAg                                                                                                                                     >  1:947960/1‑148 (MQ=255)
         cTGCAAGCGGCAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAg                                                                                                                                     <  1:1821572/148‑1 (MQ=255)
           gCACGCGGCAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGt                                                                                                                                   >  2:3505896/1‑148 (MQ=255)
           gCAAGCGGCAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGt                                                                                                                                   <  1:1030388/148‑1 (MQ=255)
               gCGGCAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTg                                                                                                                               >  2:3400266/1‑148 (MQ=255)
               gCGGCAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTg                                                                                                                               >  1:1471353/1‑148 (MQ=255)
               gCGGCAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTg                                                                                                                               >  2:1793187/1‑148 (MQ=255)
               gCGGCAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTg                                                                                                                               >  1:1471349/1‑148 (MQ=255)
               gCGGCAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTg                                                                                                                               <  2:1704979/148‑1 (MQ=255)
                cGGCAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGc                                                                                                                              <  2:630835/148‑1 (MQ=255)
                  gCAGATAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGACCTCGACCGATATTTTCGACGCCCAGGATGAGGAAAATCGGGAGACGGTGTTGCCTAATAACGTTGCGGCTCGCGTGGAGCTGGAATCAGGTATTGCCg                                                                                                                            >  2:217287/1‑148 (MQ=255)
                       gaAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTAc                                                                                                                       <  1:1001034/148‑1 (MQ=255)
                       gaAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTAc                                                                                                                       <  1:1000933/148‑1 (MQ=255)
                          aTTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTgg                                                                                                                    >  2:3535195/1‑148 (MQ=255)
                              gCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACa                                                                                                                <  1:1927018/148‑1 (MQ=255)
                              gCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACa                                                                                                                <  2:2192140/148‑1 (MQ=255)
                               cAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGACGATTACTGGTACaa                                                                                                               >  2:198844/1‑148 (MQ=255)
                                 ggAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGt                                                                                                             >  2:1417756/1‑148 (MQ=255)
                                 ggAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGt                                                                                                             >  2:268245/1‑148 (MQ=255)
                                 ggAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGt                                                                                                             <  2:529061/148‑1 (MQ=255)
                                 ggAGAAGGTCGCAATGTACACGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGt                                                                                                             >  2:1417868/1‑148 (MQ=255)
                                  gagaAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTa                                                                                                            >  1:1209322/1‑148 (MQ=255)
                                     aaGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGt                                                                                                         >  1:3011240/1‑148 (MQ=255)
                                     aaGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGt                                                                                                         >  1:3108242/1‑148 (MQ=255)
                                       ggTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTg                                                                                                       >  1:1984452/1‑148 (MQ=255)
                                       ggTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTg                                                                                                       >  1:106130/1‑148 (MQ=255)
                                          cGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTc                                                                                                    <  1:1616301/148‑1 (MQ=255)
                                            cAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTg                                                                                                  >  1:2158342/1‑148 (MQ=255)
                                             aaTGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGa                                                                                                 <  1:1949084/148‑1 (MQ=255)
                                               tGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGaaa                                                                                               >  2:1879338/1‑148 (MQ=255)
                                               tGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGaaa                                                                                               >  2:3320966/1‑148 (MQ=255)
                                                gTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAg                                                                                              >  2:389151/1‑148 (MQ=255)
                                                 tACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAgg                                                                                             <  1:717191/148‑1 (MQ=255)
                                                     cgTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCa                                                                                         <  1:2448943/148‑1 (MQ=255)
                                                       tAGTGTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAAt                                                                                       <  1:1986671/148‑1 (MQ=255)
                                                          tttCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGTGAGTCGGTGTTGCCTTCTATCGTTGCGGCTCGCGTGGAGGTTGTAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGAAATTGt                                                                                    <  1:668068/148‑1 (MQ=255)
                                                            tCCCCGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACGGGTACAAGTCTGTTGGTCTGACAGGGGCAATTGTCg                                                                                  >  1:1544410/1‑148 (MQ=255)
                                                              ccTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCggg                                                                                <  1:909871/148‑1 (MQ=255)
                                                               cTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGa                                                                               <  1:2187168/148‑1 (MQ=255)
                                                                tGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGAt                                                                              >  1:126343/1‑148 (MQ=255)
                                                                 gCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATg                                                                             >  1:985608/1‑148 (MQ=255)
                                                                      cGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACggg                                                                        <  2:1356940/148‑1 (MQ=255)
                                                                        aCCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGtt                                                                      <  2:2950480/148‑1 (MQ=255)
                                                                         ccGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTa                                                                     >  2:347303/1‑148 (MQ=255)
                                                                          cGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTAc                                                                    <  1:2739178/148‑1 (MQ=255)
                                                                          cGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTAc                                                                    <  1:963677/148‑1 (MQ=255)
                                                                           gATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACg                                                                   <  2:1961128/148‑1 (MQ=255)
                                                                            atatTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGACGGGATGACGTGTTACgg                                                                  >  2:1974784/1‑148 (MQ=255)
                                                                              atTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACgggg                                                                >  1:111207/1‑148 (MQ=255)
                                                                              atTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACgggg                                                                >  1:226098/1‑148 (MQ=255)
                                                                              atTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACgggg                                                                <  2:106130/148‑1 (MQ=255)
                                                                              atTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACgggg                                                                >  2:2470298/1‑148 (MQ=255)
                                                                                tttCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGaa                                                              >  2:3550149/1‑148 (MQ=255)
                                                                                tttCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGaa                                                              >  2:3273583/1‑148 (MQ=255)
                                                                                 ttCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAAt                                                             >  1:587299/1‑148 (MQ=255)
                                                                                  tCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATc                                                            <  2:1126351/148‑1 (MQ=255)
                                                                                     aCGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGc                                                         >  2:302452/1‑148 (MQ=255)
                                                                                       gCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTc                                                       <  1:2670624/148‑1 (MQ=255)
                                                                                                 aGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATaa                                             >  1:3205912/1‑148 (MQ=255)
                                                                                                  ggAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAg                                            <  2:2928353/148‑1 (MQ=255)
                                                                                                   gAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGc                                           <  2:1853885/148‑1 (MQ=255)
                                                                                                     atatCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTg                                         >  2:2160599/1‑148 (MQ=255)
                                                                                                      tatCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGt                                        <  2:2036242/148‑1 (MQ=255)
                                                                                                        tCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTc                                      >  2:1765281/1‑148 (MQ=255)
                                                                                                        tCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTc                                      <  1:806392/148‑1 (MQ=255)
                                                                                                         cGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTcc                                     >  1:2614047/1‑148 (MQ=255)
                                                                                                         cGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTcc                                     <  1:2634473/148‑1 (MQ=255)
                                                                                                           ggAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCg                                   >  2:2257249/1‑148 (MQ=255)
                                                                                                           ggAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCg                                   >  2:749872/1‑148 (MQ=255)
                                                                                                               tCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCt                               <  2:622028/148‑1 (MQ=255)
                                                                                                               tCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCt                               <  1:1668665/148‑1 (MQ=255)
                                                                                                               tCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCt                               <  1:3428663/148‑1 (MQ=255)
                                                                                                               tCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCt                               >  2:881669/1‑148 (MQ=255)
                                                                                                               tCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCt                               >  2:191197/1‑148 (MQ=255)
                                                                                                                      tGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCtt                        >  1:882822/1‑148 (MQ=255)
                                                                                                                      tGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCtt                        <  2:1546910/148‑1 (MQ=255)
                                                                                                                      tGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCtt                        >  1:882833/1‑148 (MQ=255)
                                                                                                                       gCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACTAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCttt                       <  2:2321352/148‑1 (MQ=255)
                                                                                                                       gCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCttt                       <  2:3108242/148‑1 (MQ=255)
                                                                                                                        ccTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTa                      >  2:3138693/1‑148 (MQ=255)
                                                                                                                        ccTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTa                      <  1:469061/148‑1 (MQ=255)
                                                                                                                        ccTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTa                      <  1:2763921/148‑1 (MQ=255)
                                                                                                                            cTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTAccgc                  <  1:389151/148‑1 (MQ=255)
                                                                                                                             tAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTAccgcc                 >  2:3129907/1‑148 (MQ=255)
                                                                                                                              aaCGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTAccgccg                <  2:2202527/148‑1 (MQ=255)
                                                                                                                               aCGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCga               >  2:1303829/1‑148 (MQ=255)
                                                                                                                                    gCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAata          <  1:1688409/148‑1 (MQ=255)
                                                                                                                                    gCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAata          <  1:3535195/148‑1 (MQ=255)
                                                                                                                                       gCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATtg       >  2:1687607/1‑148 (MQ=255)
                                                                                                                                       gCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATtg       <  2:2895624/148‑1 (MQ=255)
                                                                                                                                         tCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATtgtg     <  2:1181770/148‑1 (MQ=255)
                                                                                                                                          cgcgTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTgg    >  1:911430/1‑148 (MQ=255)
                                                                                                                                            cgTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGGCa  >  2:1668168/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
GAAATTACCCTGCAAGCGGCAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGGCA  >  CP009273/2574675‑2574962

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: