Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,074,617 C→G 17.2% M126I (ATG→ATC) ‡ tktA ← transketolase 1, thiamine triphosphate‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,074,6170CG17.2% 191.5 / 28.2 69M126I (ATG→ATC) ‡tktAtransketolase 1, thiamine triphosphate‑binding
Reads supporting (aligned to +/- strand):  ref base C (27/30);  new base G (7/5);  total (34/35)
Fisher's exact test for biased strand distribution p-value = 5.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCGGCAGATCGTAACCGGT  >  CP009273/3074469‑3074763
                                                                                                                                                    |                                                                                                                                                  
ccAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCa                                                                                                                                                   >  1:174197/1‑150 (MQ=255)
 cAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCa                                                                                                                                                   >  2:117660/1‑149 (MQ=255)
         gCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAACCTGCGCCGCCAGCGTTTTTTCTGCCATCGCCATACCGACTg                                                                                                                                          <  2:158364/150‑1 (MQ=255)
             aCCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGtt                                                                                                                                      >  2:37575/1‑150 (MQ=255)
                       gagCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATAccc                                                                                                                            <  1:42808/150‑1 (MQ=255)
                        agCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCt                                                                                                                           >  1:145666/1‑150 (MQ=255)
                         gCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTg                                                                                                                          <  1:73896/150‑1 (MQ=255)
                          cAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGa                                                                                                                         >  2:155057/1‑150 (MQ=255)
                           aagcTTCATGCGAGATGCCTTCCATTAAGCAGCCGTCGCCCAGGAAGGCGTAGGTGTAGTGGTCGACGATTTCGTGGCCGGGGCGGTTGAATTGGCCCGCCAGCGCGCGTTCGGCTATCGCGAAACCCACGGCATTGGTGATGCCTTGAc                                                                                                                        >  1:122710/4‑150 (MQ=255)
                            agcTTCATGCGAGATGCCTTCCATTAAGCAGCCGTCGCCCAGGAAGGCGTAGGTGTAGTGGTCGACGATTTCGTGGCCGGGGCGGTTGAATTGGCCCGCCAGCGCGCGTTCGGCTATCGCGAAACCCACGGCATTGGTGATGCCTTGAcc                                                                                                                       >  1:8257/3‑150 (MQ=255)
                             aCTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGAcc                                                                                                                       >  1:140509/1‑149 (MQ=255)
                              cTTCATGCGAGATGCCTTCCATTAAGCAGCCGTCGCCCAGGAAGGCGTAGGTGTAGTGGTCGACGATTTCGTGGCCGGGGCGGTTGAATTGGCCCGCCAGCGCGCGTTCGGCTATCGCGAAACCCACGGCATTGGTGATGCCTTGAccta                                                                                                                     <  2:182026/150‑3 (MQ=255)
                              cTTCATGCGAGATGCCTTCCATTAAGCAGCCGTCGCCCAGGAAGGCGTAGGTGTAGTGGTCGACGATTTCGTGGCCGGGGCGGTTGAATTGGCCCGCCAGCGCGCGTTCGGCTATCGCGAAACCCACGGCATTGGTGATGCCTTGAccta                                                                                                                     <  2:75445/150‑3 (MQ=255)
                                     cgagaTGCCTTCCATTAAGCAGCCGTCGCCCAGGAAGGCGTAGGTGTAGTGGTCGACGATTTCGTGGCCGGGGCGGTTGAATTGGCCCGCCAGCGCGCGTTCGGCTATCGCGAAACCCACGGCATTGGTGATGCCTTGACCTAACGGTcc                                                                                                              <  1:4254/149‑1 (MQ=255)
                                      gagaTGCCTTCCATTAAGCAGCCGTCGCCCAGGAAGGCGTAGGTGTAGTGGTCGACGATTTCGTGGCCGGGGCGGTTGAATTGGCCCGCCAGCGCGCGTTCGGCTATCGCGAAACCCACGGCATTGGTGATGCCTTGACCTAACGGTCCg                                                                                                             >  2:17983/1‑150 (MQ=255)
                                        gaTGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCggt                                                                                                           >  1:175010/1‑150 (MQ=255)
                                                 cattaagCAGCCGTCGCCCAGGAAGGCGTAGGTGTAGTGGTCGACGATTTCGTGGCCGGGGCGGTTGAATTGGCCCGCCAGCGCGCGTTCGGCTATCGCGAAACCCACGGCATTGGTGATGCCTTGACCTAACGGTCCGGTCGTCGTTTc                                                                                                  <  2:92247/144‑1 (MQ=255)
                                                 catcatGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGACCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTc                                                                                                  >  2:169259/1‑150 (MQ=255)
                                                     aagCAGCCGTCGCCCAGGAAGGCGTAGGTGTAGTGGTCGACGATTTCGTGGCCGGGGCGGTTGAATTGGCCCGCCAGCGCGCGTTCGGCTATCGCGAAACCCACGGCATTGGTGATGCCTTGACCTAACGGTCCGGTCGTCGTTTCgacg                                                                                              >  1:122517/3‑149 (MQ=255)
                                                     aagCAGCCGTCGCCCAGGAAGGCGTAGGTGTAGTGGTCGACGATTTCGTGGCCGGGGCGGTTGAATTGGCCCGCCAGCGCGCGTTCGGCTATCGCGAAACCCACGGCATTGGTGATGCCTTGACCTAACGGTCCGGTCGTCGTTTCgac                                                                                               <  1:149637/147‑1 (MQ=255)
                                                       gCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCacac                                                                                             >  2:141851/1‑149 (MQ=255)
                                                        cAGCCGTCGCCCAGGAAGGCGTAGGTGTAGTGGTCGACGATTTCGTGGCCGGGGCGGTTGAATTGGCCCGCCCGCGCGCGTTCGGCTATCGCGAAACCCACGGCATTGGTGATGCCTTGACCTAACGGTCCGGTCGTCGTTTCGACGccc                                                                                           >  2:71120/1‑149 (MQ=255)
                                                        cAGCCGTCGCCCAGGAAGGCGTAGGTGTAGTGGTCGACGATTTCGTGGCCGGGGCGGTTGAATTGGCCCGCCAGCGCGCGTTCGGCTATCGCGAAACCCACGGCATTGGTGATGCCTTGACCTAACGGTCCGGTCGTCGTTTCGACGccc                                                                                           >  1:5291/1‑149 (MQ=255)
                                                             gTCGCCCAGGAAGGCGTAGGTGTAGTGGTCGACGATTTCGTGGCCGGGGCGGTTGAATTGGCCCGCCCGCGCGCGTTCGGCTATCGCGAAACCCACCGCATTGGTGATGCCTTGACCTAACGGTCCGGTCGTCGTTTCGACGcccggcg                                                                                       >  1:69933/1‑144 (MQ=255)
                                                              tCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGgtg                                                                                     <  2:175010/150‑1 (MQ=255)
                                                              tCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGgtg                                                                                     <  2:50836/150‑1 (MQ=255)
                                                                gCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGgtgt                                                                                    >  2:125804/1‑149 (MQ=255)
                                                                 cccATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTaa                                                                                  >  2:164962/1‑150 (MQ=255)
                                                                  ccATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAAc                                                                                 <  1:618/150‑1 (MQ=255)
                                                                   cATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAAcc                                                                                >  1:103261/1‑150 (MQ=255)
                                                                    aTGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAAcc                                                                                >  1:87668/1‑149 (MQ=255)
                                                                      gAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCAc                                                                             >  1:75752/1‑150 (MQ=255)
                                                                      gAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCAc                                                                             >  2:170237/1‑150 (MQ=255)
                                                                              aGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGt                                                                     >  2:66394/1‑150 (MQ=255)
                                                                               ggTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGt                                                                     <  1:99941/149‑1 (MQ=255)
                                                                                   tAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGAcc                                                                 <  1:25703/149‑1 (MQ=255)
                                                                                        gTCGACAATGTCGTGGCCCGGACGGTTAAGCTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATGCCCTGACCCAGCGGACCGGTGGTGGTGTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGt                                                           <  2:41066/150‑1 (MQ=255)
                                                                                           gACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGtttt                                                        <  1:155057/150‑1 (MQ=255)
                                                                                              aaTGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTaga                                                     <  1:116709/150‑1 (MQ=255)
                                                                                                tGTCGTGGCCCGGGCGGTTAAACTGCGCCGCCATTTTTTTTTCTACAGTCGCCATACCGACTGCGTTGGCAATCCCCTGACCCAGCGCGCCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGt                                                   <  1:157359/150‑1 (MQ=255)
                                                                                                tGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGt                                                   <  2:103261/150‑1 (MQ=255)
                                                                                                       gCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTg                                            <  2:145666/150‑1 (MQ=255)
                                                                                                          cGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACg                                         <  1:11163/150‑1 (MQ=255)
                                                                                                          cGGAAGGGTAAACGGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACg                                         <  1:117660/150‑1 (MQ=255)
                                                                                                              cGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAg                                     >  1:105958/1‑150 (MQ=255)
                                                                                                               ggTTAAGCTGCGCCGCCAGCGTTTTTTCTGCAGTCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCGCTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGt                                    <  2:83421/150‑1 (MQ=255)
                                                                                                               ggTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGt                                    >  2:10619/1‑150 (MQ=255)
                                                                                                                    aaCTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTc                               <  1:151025/150‑1 (MQ=255)
                                                                                                                         ggccgcaAGCGTGTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGttc                          >  2:38399/2‑150 (MQ=255)
                                                                                                                         cgccgcCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGtt                           <  1:132315/149‑1 (MQ=255)
                                                                                                                            cgccAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGttcttc                       >  2:133284/1‑150 (MQ=255)
                                                                                                                             gccAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTcc                      <  2:87668/150‑1 (MQ=255)
                                                                                                                             gccAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTcc                      <  2:176723/150‑1 (MQ=255)
                                                                                                                              ccAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCa                     <  2:140509/150‑1 (MQ=255)
                                                                                                                              ccAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCa                     >  2:168703/1‑150 (MQ=255)
                                                                                                                               cAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCAt                    >  1:145218/1‑150 (MQ=255)
                                                                                                                                aGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATc                   >  2:10998/1‑150 (MQ=255)
                                                                                                                                 gCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCg                  <  1:170237/150‑1 (MQ=255)
                                                                                                                                   gTTTTTTCTGCAGTCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCgg                 >  1:118285/1‑149 (MQ=255)
                                                                                                                                   gTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCgg                 >  1:21110/1‑149 (MQ=255)
                                                                                                                                    ttttttCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCgg                 >  2:46530/1‑148 (MQ=255)
                                                                                                                                      ttttCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCGGCAg              <  1:141851/149‑1 (MQ=255)
                                                                                                                                         tCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCGGCAGATc           <  1:49841/149‑1 (MQ=255)
                                                                                                                                           tGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCGGCAGATCGt         >  2:100517/1‑149 (MQ=255)
                                                                                                                                             cAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCGGCAGATCGTa        <  2:135054/148‑1 (MQ=255)
                                                                                                                                             cAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCGGCAGATCGTa        <  1:155328/148‑1 (MQ=255)
                                                                                                                                              aaTCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCGGCAGATCGTAAcc     <  1:37575/150‑1 (MQ=255)
                                                                                                                                               aTCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCGGCAGATCGTAACCg    <  1:125804/150‑1 (MQ=255)
                                                                                                                                                  gCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCGGCAGATCGTAACCGGt  <  2:152842/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCGGCAGATCGTAACCGGT  >  CP009273/3074469‑3074763

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: