Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,574,820 C→A 86.6% 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,8200CA86.6% 421.1 / 71.7 187A609E (GCG→GAG) tktBtransketolase 2, thiamine triphosphate‑binding
Reads supporting (aligned to +/- strand):  ref base C (8/17);  new base A (70/92);  total (78/109)
Fisher's exact test for biased strand distribution p-value = 3.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.57e-01

Only 100 of 187 total aligned reads displayed.

AAATTACCCTGCAAGCGGCAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGGCAAA  >  CP009273/2574676‑2574964
                                                                                                                                                |                                                                                                                                                
    tACCCTGCAAGCGGCAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGaa                                                                                                                                           <  2:160453/148‑1 (MQ=255)
       ccTGCAAGCGGCAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCa                                                                                                                                        <  1:766522/148‑1 (MQ=255)
        cTGCAAGCGGCAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAg                                                                                                                                       >  2:695729/1‑148 (MQ=255)
        cTGCAAGCGGCAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAg                                                                                                                                       >  1:1960753/1‑148 (MQ=255)
        cTGCAAGCGGCAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAg                                                                                                                                       >  1:1713621/1‑148 (MQ=255)
         tGCAAGCGGCAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAgg                                                                                                                                      >  1:2253230/1‑148 (MQ=255)
          gCAAGCGGCAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGt                                                                                                                                     <  1:1394858/148‑1 (MQ=255)
              gCGGCAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTg                                                                                                                                 <  1:3036013/148‑1 (MQ=255)
              gCGGCAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTg                                                                                                                                 <  1:3035961/148‑1 (MQ=255)
               cGGCAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGc                                                                                                                                >  1:1345405/1‑148 (MQ=255)
                 gCAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCg                                                                                                                              <  1:2433732/148‑1 (MQ=255)
                  cAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGa                                                                                                                             <  1:3092247/148‑1 (MQ=255)
                  cAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGa                                                                                                                             <  1:3092242/148‑1 (MQ=255)
                   agagAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATTTTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGAt                                                                                                                            <  2:1181519/148‑1 (MQ=255)
                   agagAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATAGTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGACGCAGGTATTGCCGAt                                                                                                                            >  1:64417/1‑148 (MQ=255)
                     agaAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTa                                                                                                                          >  1:2980630/1‑148 (MQ=255)
                      gaAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTAc                                                                                                                         <  2:1395647/148‑1 (MQ=255)
                       aaaTTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACt                                                                                                                        >  2:2177127/1‑148 (MQ=255)
                           tAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTa                                                                                                                    >  2:3250626/1‑148 (MQ=255)
                             gCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACa                                                                                                                  >  1:2188453/1‑148 (MQ=255)
                                ggAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGt                                                                                                               >  1:2732930/1‑148 (MQ=255)
                                 gagaAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTa                                                                                                              <  1:1041400/148‑1 (MQ=255)
                                 gagaAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTa                                                                                                              <  2:1787139/148‑1 (MQ=255)
                                    aaGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGt                                                                                                           <  1:1540231/148‑1 (MQ=255)
                                       gTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTgg                                                                                                        <  1:1842427/148‑1 (MQ=255)
                                       gTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTgg                                                                                                        <  1:365918/148‑1 (MQ=255)
                                        tCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGt                                                                                                       >  2:654028/1‑148 (MQ=255)
                                         cGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTc                                                                                                      >  1:155805/1‑148 (MQ=255)
                                               gTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAg                                                                                                <  1:2553546/148‑1 (MQ=255)
                                                tACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAgg                                                                                               >  2:1294743/1‑148 (MQ=255)
                                                tACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAgg                                                                                               <  2:826811/148‑1 (MQ=255)
                                                    cgTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCa                                                                                           >  2:1969196/1‑148 (MQ=255)
                                                      tAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAAt                                                                                         >  2:2610637/1‑148 (MQ=255)
                                                      tAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAAt                                                                                         <  1:40052/148‑1 (MQ=255)
                                                       aGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAAtt                                                                                        >  2:933086/1‑148 (MQ=255)
                                                            cccTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCgg                                                                                   >  1:98526/1‑148 (MQ=255)
                                                              cTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGa                                                                                 <  1:2384346/148‑1 (MQ=255)
                                                                gCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATg                                                                               <  1:3064959/148‑1 (MQ=255)
                                                                gCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATg                                                                               <  2:1618301/148‑1 (MQ=255)
                                                                   cTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACg                                                                            >  2:1848951/1‑148 (MQ=255)
                                                                   cTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACg                                                                            <  2:2843588/148‑1 (MQ=255)
                                                                    tCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACgg                                                                           <  2:1137082/148‑1 (MQ=255)
                                                                     cGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACggg                                                                          <  2:279586/148‑1 (MQ=255)
                                                                     cGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACggg                                                                          <  1:3178553/148‑1 (MQ=255)
                                                                       aCCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGtt                                                                        >  1:27331/1‑148 (MQ=255)
                                                                       aCCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGtt                                                                        >  2:2089436/1‑148 (MQ=255)
                                                                         cGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTAc                                                                      >  1:1373655/1‑148 (MQ=255)
                                                                         cGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTAc                                                                      >  1:1373710/1‑148 (MQ=255)
                                                                         cGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTAc                                                                      >  1:909877/1‑148 (MQ=255)
                                                                          gATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACg                                                                     <  2:1129235/148‑1 (MQ=255)
                                                                          gATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACg                                                                     <  2:1129549/148‑1 (MQ=255)
                                                                          gATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACg                                                                     <  2:1129649/148‑1 (MQ=255)
                                                                          gATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACg                                                                     <  2:855376/148‑1 (MQ=255)
                                                                               tttCGACGCCCAGGATGAGGAAGATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGaa                                                                <  2:867743/148‑1 (MQ=255)
                                                                                ttCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAAt                                                               <  1:591397/148‑1 (MQ=255)
                                                                                ttCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAAt                                                               <  1:560618/148‑1 (MQ=255)
                                                                                ttCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAAt                                                               <  2:2588499/148‑1 (MQ=255)
                                                                                     cGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCt                                                          >  2:396150/1‑148 (MQ=255)
                                                                                      gCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTc                                                         <  2:1464864/148‑1 (MQ=255)
                                                                                      gCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTc                                                         <  1:1472950/148‑1 (MQ=255)
                                                                                      gCCCAGGATGAGGAATATCGGGAGTCGGTGGTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTc                                                         <  2:3106320/148‑1 (MQ=255)
                                                                                       cccAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTcc                                                        >  2:1923428/1‑148 (MQ=255)
                                                                                        ccAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCcg                                                       >  2:1745936/1‑148 (MQ=255)
                                                                                         cAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCcgg                                                      >  2:2574090/1‑148 (MQ=255)
                                                                                            gATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCcggcgg                                                   >  1:3074473/1‑148 (MQ=255)
                                                                                            gATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCcggcgg                                                   >  2:898635/1‑148 (MQ=255)
                                                                                               gAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATa                                                >  2:2038263/1‑148 (MQ=255)
                                                                                                       tCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTc                                        <  1:1537251/148‑1 (MQ=255)
                                                                                                             gTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTc                                  <  2:487775/148‑1 (MQ=255)
                                                                                                             gTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTc                                  <  1:2089436/148‑1 (MQ=255)
                                                                                                              tCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCt                                 <  2:2257296/148‑1 (MQ=255)
                                                                                                              tCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCt                                 <  1:2230216/148‑1 (MQ=255)
                                                                                                              tCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCt                                 <  2:451952/148‑1 (MQ=255)
                                                                                                              tCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCt                                 >  2:2750333/1‑148 (MQ=255)
                                                                                                                ggTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCttt                               <  1:2660438/148‑1 (MQ=255)
                                                                                                                ggTGTTGCCTTCTAACGGTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCttt                               <  2:3100955/148‑1 (MQ=255)
                                                                                                                 gtgtTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTg                              <  1:847637/148‑1 (MQ=255)
                                                                                                                 gtgtTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTg                              <  1:847706/148‑1 (MQ=255)
                                                                                                                  tgtTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTgg                             >  1:921273/1‑148 (MQ=255)
                                                                                                                     tGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCtt                          <  1:1294743/148‑1 (MQ=255)
                                                                                                                      gCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCttt                         >  1:45077/1‑148 (MQ=255)
                                                                                                                       ccTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTa                        <  2:1080008/148‑1 (MQ=255)
                                                                                                                       ccTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTa                        <  2:1554107/148‑1 (MQ=255)
                                                                                                                        cTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTAc                       >  2:1694578/1‑148 (MQ=255)
                                                                                                                             aaCGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTAccgccg                  >  2:1983048/1‑148 (MQ=255)
                                                                                                                             aaCGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTAccgccg                  >  2:834758/1‑148 (MQ=255)
                                                                                                                               cGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCgag                >  2:1671472/1‑148 (MQ=255)
                                                                                                                                gTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCgaga               <  1:90652/148‑1 (MQ=255)
                                                                                                                                   gCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAata            >  1:1316217/1‑148 (MQ=255)
                                                                                                                                   gCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAata            <  1:2045912/148‑1 (MQ=255)
                                                                                                                                    cGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAatat           >  2:1834354/1‑148 (MQ=255)
                                                                                                                                    cGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAatat           >  2:2292219/1‑148 (MQ=255)
                                                                                                                                    cGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAatat           <  1:2175784/148‑1 (MQ=255)
                                                                                                                                      gCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATtg         >  2:557613/1‑148 (MQ=255)
                                                                                                                                       cTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATtgt        <  1:401141/148‑1 (MQ=255)
                                                                                                                                        tCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATtgtg       <  1:2423360/148‑1 (MQ=255)
                                                                                                                                         cgcgTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCACGTTCTTTGGCTTTACAGACGAGAATATTGTgg      >  2:2585303/1‑148 (MQ=255)
                                                                                                                                           cgTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGGCa    >  2:53139/1‑148 (MQ=255)
                                                                                                                                             tGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGGCaaa  <  1:1032179/148‑1 (MQ=255)
                                                                                                                                             tGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGGCaaa  <  1:1031084/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                
AAATTACCCTGCAAGCGGCAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGGCAAA  >  CP009273/2574676‑2574964

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: