Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,074,629 G→A 18.0% N122N (AAC→AAT tktA ← transketolase 1, thiamine triphosphate‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,074,6290GA18.0% 1021.3 / 164.3 430N122N (AAC→AATtktAtransketolase 1, thiamine triphosphate‑binding
Reads supporting (aligned to +/- strand):  ref base G (177/175);  new base A (40/38);  total (217/213)
Fisher's exact test for biased strand distribution p-value = 9.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.51e-05

Only 100 of 435 total aligned reads displayed.

GTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCGGCAGATCGTAACCGGTGAGGTGCAGCAGGCT  >  CP009273/3074480‑3074778
                                                                                                                                                     |                                                                                                                                                     
 tACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGt                                                                                                                                                      >  2:37170/1‑150 (MQ=255)
 tACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGt                                                                                                                                                      >  1:425362/1‑150 (MQ=255)
    ccGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTgg                                                                                                                                                   <  1:757132/150‑1 (MQ=255)
    ccGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTgg                                                                                                                                                   <  2:204033/150‑1 (MQ=255)
       ccaatgagCAAGCTTCATGCGAGATGCCTGCCATTAAGCAGCCGTCGCCCAGGAAGGCGTAGGTGTAGTGGTCGACGATTTCGTGGCCGGGGCGGTTGAATTGGCCCGCCAGCGCGCGTTCGGCTATCGCGAAACCCACGGCATTggtg                                                                                                                                                 <  1:210482/144‑3 (MQ=255)
       ccAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCaa                                                                                                                                                <  2:411669/150‑1 (MQ=255)
         agagagCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATa                                                                                                                                              <  1:622111/150‑1 (MQ=255)
         agagagCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATa                                                                                                                                              >  2:833302/1‑150 (MQ=255)
           agagCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATATCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATAcc                                                                                                                                            <  1:977927/150‑1 (MQ=255)
             agCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCt                                                                                                                                          <  1:637827/150‑1 (MQ=255)
                  gcTTCATGCGAGATGCCTTCCATTAAGCAGCCGTCGCCCAGGAAGGCGTAGGTGTAGTGGTCGACGATTTCGTGGCCGGGGCGGTTGAATTGGCCCGCCAGCGCGCGTTCGGCTATCGCGAAACCCACGGCATTGGTGATGCCTTGAcc                                                                                                                                      >  2:692078/2‑149 (MQ=255)
                  aCTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGAccc                                                                                                                                     <  1:440340/150‑1 (MQ=255)
                  aCTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGAcc                                                                                                                                      >  1:888428/1‑149 (MQ=255)
                       gTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCg                                                                                                                                 <  1:830832/149‑1 (MQ=255)
                        tGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGa                                                                                                                               >  2:867344/1‑150 (MQ=255)
                          cgagaTGCCTTCCATTAAGCAGCCGTCGCCCAGGAAGGCGTAGGTGTAGTGGTCGACGATTTCGTGGCCGGGGCGGTTGAATTGGCCCGCCAGCGCGCGTTCGGCTATCGCGAAACCCACGGCATTGGTGATGCCTTGACCTAACGGTcc                                                                                                                             <  1:809091/149‑1 (MQ=255)
                           gagaTGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCg                                                                                                                            >  2:274525/1‑150 (MQ=255)
                             gaTGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCggt                                                                                                                          >  1:40317/1‑150 (MQ=255)
                                  cTTCCATTAAGCAGCCGTCGCCCAGGAAGGCGTAGGTGTAGTGGTCGACGATTTCGTGGCCGGGGCGGTTGAATTGGCCCGCCAGCGCGCGTTCGGCTATCGCGAAACCCACGGCATTGGGGATGCCTTTACCTAACGGTCCggtcgtcg                                                                                                                     >  2:461402/1‑148 (MQ=255)
                                  cTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCggtggtgg                                                                                                                     <  1:652650/150‑1 (MQ=255)
                                   ttCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGCTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCggtggtggt                                                                                                                    >  2:605558/1‑150 (MQ=255)
                                     ccATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCggtggtggt                                                                                                                    <  1:555674/148‑1 (MQ=255)
                                      cattaagCAGCCGTCGCCCAGGAAGGCGTAGGTGTAGTGGTCGACGATTTCGTGGCCGGGGCGGTTGAATTGGCCCGCCAGCGCGCGTTCGGCTATCGCGAAACCCACGGCATTGGTGATGCCTTGACCTAACGGTCCGGTCGTCGTTTc                                                                                                                 <  1:567787/144‑1 (MQ=255)
                                      catcatGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTc                                                                                                                 >  2:192550/1‑150 (MQ=255)
                                       atcatGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTc                                                                                                                 >  1:772343/1‑149 (MQ=255)
                                        ttaagCAGCCGTCGCCCAGGAAGGCGTAGGTGTAGTGGTCGACGATTTCGTGGCCGGGGCGGTTGAATTGGCCCGCCAGCGCGCGTTCGGCTATCGCGAAACCCACGGCATTGGGGATGCCCTGACCTAACGGTCCGGTCGTCGGTTcga                                                                                                               >  2:849048/5‑148 (MQ=255)
                                        tcatGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCca                                                                                                               >  2:563642/1‑150 (MQ=255)
                                          aagCAGCCGTCGCCCAGGAAGGCGTAGGTGTAGTGGTCGACGATTTCGTGGCCGGGGCGGTTGAATTGGCCCGCCAGCGCGCGTTCGGCTATCGCGAAACCCACGGCATTGGTGATGCCTTGACCTAACGGTCCGGTCGTCGTTTCgac                                                                                                              <  1:110087/147‑1 (MQ=255)
                                          aagCAGCCGTCGCCCAGGAAGGCGTAGGTGTAGTGGTCGACGATTTCGTGGCCGGGGCGGTTGAATTGGCCCGCCAGCGCGCGTTCGGCTATCGCGAAACCCACGGCATTGGTGATGCCTTGACCTAACGGTCCGGTCGTCGTTTCgac                                                                                                              <  1:320780/147‑1 (MQ=255)
                                            gCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCacac                                                                                                            >  1:456434/1‑149 (MQ=255)
                                             cAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACAcc                                                                                                           <  1:949365/149‑1 (MQ=255)
                                             cAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCa                                                                                                          <  1:895040/150‑1 (MQ=255)
                                              aGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAg                                                                                                         <  2:161321/150‑1 (MQ=255)
                                               gCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGc                                                                                                        <  1:1010835/150‑1 (MQ=255)
                                               gCCGTCGCCCAGGAAGGCGTAGGTGTAGTGGTCGACGATTTCGTGGCCGGGGCGGTTGAATTGGCCCGCCAGCGCGCGTTCGGCTATCGCGAAACCCACGGCATTGGTGATGCCTTGACCTAACGGTCCGGTCGTCGTTTCGACGcccgg                                                                                                        <  1:173341/150‑4 (MQ=255)
                                                ccGTCGCCCAGGAAGGCGTAGGTGTAGTGGTCGACGATTTCGTGGCCGGGGCGGTTGAATTGGCCCGCCAGCGCGCGTTCGGCTATCGCGAAACCCACGGCATTGGTGATGCCTTGACCTAACGGTCCGGTCGTCGTTTCGACGcccggc                                                                                                       <  1:743526/150‑5 (MQ=255)
                                                 cGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCgg                                                                                                      >  2:245065/1‑150 (MQ=255)
                                                 cGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCgg                                                                                                      >  2:813124/1‑150 (MQ=255)
                                                  gTCGCCCAGGAAGGCGTAGGTGTAGTGGTCGACGATTTCGTGGCCGGGGCGGTTGAATTGGCCCGCCAGCGCGCGTTCGGCTATCGCGAAACCCACGGCATTGGTGATGCCTTGACCTAACGGTCCGGTCGTCGTTTCGACGcccggcg                                                                                                      >  2:426442/1‑144 (MQ=255)
                                                  gTCGCCCAGGAAGGCGTAGGTGTAGTGGTCGACGATTTCGTGGCCGGGGCGGTTGAATTGGCCCGCCAGCGCGCGTTCGGCTATCGCGAAACCCACGGCATTGGTGATGCCTTGACCTAACGGTCCGGTCGTCGTTTCGACGcccggcg                                                                                                      >  2:62954/1‑144 (MQ=255)
                                                    cGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGgtg                                                                                                    >  1:841985/1‑149 (MQ=255)
                                                     gCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTa                                                                                                  <  1:985396/150‑1 (MQ=255)
                                                     gCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGCTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGgtgt                                                                                                   >  2:509046/1‑149 (MQ=255)
                                                        caggAAGGGGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGCACGGTTAAGCTGCGCCGCCAGCCTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAAcc                                                                                               <  2:956136/147‑1 (MQ=255)
                                                        cATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAAcc                                                                                               >  1:68217/1‑150 (MQ=255)
                                                        cATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAAcc                                                                                               <  2:771762/150‑1 (MQ=255)
                                                         aTGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAAcc                                                                                               >  1:728066/1‑149 (MQ=255)
                                                           gAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCAc                                                                                            >  2:285542/1‑150 (MQ=255)
                                                            aaGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACt                                                                                           <  2:984546/150‑1 (MQ=255)
                                                            aaGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACt                                                                                           >  1:311109/1‑150 (MQ=255)
                                                             aggcgtaggcgtAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACtt                                                                                          >  2:905938/1‑150 (MQ=255)
                                                             aGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACtt                                                                                          <  2:667180/150‑1 (MQ=255)
                                                                cGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCg                                                                                       >  1:1029073/1‑150 (MQ=255)
                                                                   aGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGt                                                                                    >  2:304143/1‑150 (MQ=255)
                                                                    ggTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGt                                                                                    <  2:521350/149‑1 (MQ=255)
                                                                    ggTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTg                                                                                   <  1:930010/150‑1 (MQ=255)
                                                                       gtAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGAcc                                                                                >  1:330466/1‑150 (MQ=255)
                                                                       gtAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGAcc                                                                                >  2:100462/1‑150 (MQ=255)
                                                                         aGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCg                                                                              >  1:451580/1‑150 (MQ=255)
                                                                           tGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGa                                                                            <  2:38298/150‑1 (MQ=255)
                                                                             gTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGTCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGt                                                                          >  2:397356/1‑150 (MQ=255)
                                                                                  cAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTag                                                                     <  2:109353/150‑1 (MQ=255)
                                                                                  cAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTag                                                                     <  2:1029073/150‑1 (MQ=255)
                                                                                     tGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGt                                                                  <  2:651087/150‑1 (MQ=255)
                                                                                      gTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTg                                                                 >  2:299231/1‑150 (MQ=255)
                                                                                       tCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTg                                                                 <  2:915841/149‑1 (MQ=255)
                                                                                       tCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTg                                                                 <  2:973320/149‑1 (MQ=255)
                                                                                          tGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGc                                                             <  2:828585/150‑1 (MQ=255)
                                                                                            gCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCt                                                            >  1:249878/1‑149 (MQ=255)
                                                                                            gCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTg                                                           <  1:192550/150‑1 (MQ=255)
                                                                                            gCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTg                                                           <  1:865816/150‑1 (MQ=255)
                                                                                                    ggTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGt                                                   >  2:925395/1‑150 (MQ=255)
                                                                                                    ggTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGt                                                   <  2:450637/150‑1 (MQ=255)
                                                                                                     gTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGtt                                                  >  1:720347/1‑150 (MQ=255)
                                                                                                        aaaCTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGttttt                                               >  2:550234/1‑150 (MQ=255)
                                                                                                        aaaCTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGtttt                                                <  1:1000151/149‑1 (MQ=255)
                                                                                                         aaCTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGttttt                                               <  2:750295/149‑1 (MQ=255)
                                                                                                         aaCTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTc                                              <  2:579278/150‑1 (MQ=255)
                                                                                                         aaCTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTc                                              <  2:332364/150‑1 (MQ=255)
                                                                                                              cgccgcCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGttc                                         >  1:11528/1‑150 (MQ=255)
                                                                                                              cgccgcCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGt                                           <  2:592756/148‑1 (MQ=255)
                                                                                                               gccgccAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGttct                                        <  1:530753/150‑1 (MQ=255)
                                                                                                               gccgccAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGttct                                        >  1:513653/1‑150 (MQ=255)
                                                                                                                 cgccAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGttcttc                                      >  2:712857/1‑150 (MQ=255)
                                                                                                                   ccAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCa                                    <  2:736922/150‑1 (MQ=255)
                                                                                                                    cAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCAt                                   <  1:297713/150‑1 (MQ=255)
                                                                                                                    cAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCAt                                   >  1:351390/1‑150 (MQ=255)
                                                                                                                    cAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCAt                                   >  2:315659/1‑150 (MQ=255)
                                                                                                                      gCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCg                                 <  1:151830/150‑1 (MQ=255)
                                                                                                                         ttttttCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCGGCa                              <  1:274525/150‑1 (MQ=255)
                                                                                                                          tttttCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCgg                                >  2:53236/1‑147 (MQ=255)
                                                                                                                           ttttCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCGGCAg                             <  1:376725/149‑1 (MQ=255)
                                                                                                                                 gCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCGGCAGATCGTaa                      >  1:59650/1‑150 (MQ=255)
                                                                                                                                        ccATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCGGCAGATCGTAACCGGTGa               <  1:489440/150‑1 (MQ=255)
                                                                                                                                         cATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCGGCAGATCGTAACCGGTGAg              >  2:147414/1‑150 (MQ=255)
                                                                                                                                          aTACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCGGCAGATCGTAACCGGTGAgg             <  2:351390/150‑1 (MQ=255)
                                                                                                                                                aCTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCGGCAGATCGTAACCGGTGAGGTgcag        >  2:22846/1‑149 (MQ=255)
                                                                                                                                                    cGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCGGCAGATCGTAACCGGTGAGGTGCAGCAGGc   <  2:884914/150‑1 (MQ=255)
                                                                                                                                                    cGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCGGCAGATCGTAACCGGTGAGGTGCAGCAGGc   >  2:173515/1‑150 (MQ=255)
                                                                                                                                                     gTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCGGCAGATCGTAACCGGTGAGGTGCAGCAGGCt  >  1:133862/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCGGCAGATCGTAACCGGTGAGGTGCAGCAGGCT  >  CP009273/3074480‑3074778

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