Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,074,538 T→G 32.3% M153L (ATG→CTG)  tktA ← transketolase 1, thiamine triphosphate‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,074,5380TG32.3% 80.7 / 67.5 65M153L (ATG→CTG) tktAtransketolase 1, thiamine triphosphate‑binding
Reads supporting (aligned to +/- strand):  ref base T (26/18);  new base G (10/11);  total (36/29)
Fisher's exact test for biased strand distribution p-value = 4.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

CGGTGTCGTCGGTGAACCAGCCTTCAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACC  >  CP009273/3074391‑3074685
                                                                                                                                                   |                                                                                                                                                   
cGGTGTCGTCGGTGAACCAGCCTTCAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATg                                                                                                                                                    <  1:77725/149‑1 (MQ=255)
       gtcgGTGAACCAGCCTTCAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCg                                                                                                                                              >  1:132846/1‑148 (MQ=255)
       gtcgGTGAACCAGCCTTCAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTa                                                                                                                                            <  2:104852/150‑1 (MQ=255)
            tGAACCAGCCTTCAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGgtgt                                                                                                                                       <  1:185967/150‑1 (MQ=255)
                   gCCTTCAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTc                                                                                                                                >  1:152882/1‑150 (MQ=255)
                    ccTTCAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTc                                                                                                                                >  2:49841/1‑149 (MQ=255)
                        cAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGAc                                                                                                                             <  2:81343/148‑1 (MQ=255)
                        cAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACa                                                                                                                            <  1:170380/149‑1 (MQ=255)
                         aaCGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTa                                                                                                                                      >  2:115508/1‑138 (MQ=255)
                          aCGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAg                                                                                                                                     <  1:115508/138‑1 (MQ=255)
                          aCGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATg                                                                                                                         <  1:83282/150‑1 (MQ=255)
                               aCCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGt                                                                                                                     <  2:32942/149‑1 (MQ=255)
                                   tCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGccc                                                                                                                <  2:98425/150‑1 (MQ=255)
                                   tCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGccc                                                                                                                >  1:69388/1‑150 (MQ=255)
                                    cGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCg                                                                                                               >  2:60793/1‑150 (MQ=255)
                                       tAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGAc                                                                                                            >  2:151025/1‑150 (MQ=255)
                                        aGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACg                                                                                                           >  2:73896/1‑150 (MQ=255)
                                                   tGTCGTCGTAGACGGCAATCAGTTTGCCCAGTTTCATCGAGCCTGCCAATGAGCAAGCTTCATGCGAGATGCCTTCCATTAAGCAGCCGTCGCCCAGGAAGGCGTAGGTGTAGTGGTCGACGATTTCGTGGCCGGGGCGGTTGAATTggc                                                                                                >  1:75445/1‑148 (MQ=255)
                                                    gtcgtcGTAGACGGCAATCAGTTTGCCCAGTTTCATCGAGCCTGCCAATGAGCAAGCTTCATGCGAGATGCCTTCCATTAAGCAGCCGTCGCCCAGGAAGGCGTAGGTGTAGTGGTCGACGATTTCGTGGCCGGGGCGGTTGAATTggc                                                                                                <  1:194233/149‑3 (MQ=255)
                                                    gTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTgcgc                                                                                               >  2:157359/1‑150 (MQ=255)
                                                     tcgtCGTAGACGGCAATCAGTTTGCCCAGTTTCATCGAGCCTGCCAATGAGCAAGCTTCATGCGAGATGCCTTCCATTAAGCAGCCGTCGCCCAGGAAGGCGTAGGTGTAGTGGTCGACGATTTCGTGGCCGGGGCGGTTGAATTggccc                                                                                              <  2:165289/147‑5 (MQ=255)
                                                       gtCGTAGACGGCAATCAGTTTGCCCAGTTTCATCGAGCCTGCCAATGAGCAAGCTTCATGCGAGATGCCTTCCATTAAGCAGCCGTCGCCCAGGAAGGCGTAGGTGTAGTGGTCGACGATTTCGTGGCCGGGGCGGTTGAATTggcccg                                                                                             <  1:80840/148‑6 (MQ=255)
                                                          gTAGACGGCAATCAGTTTGCCCAGTTTCATCGAGCCTGCCAATGAGCAAGCTTCATGCGAGATGCCTTCCATTAAGCAGCCGTCGCCCAGGAAGGCGTAGGTGTAGTGGTCGACGATTTCGTGGCCGGGGCGGTTGAATTGGCCCGCCAg                                                                                         >  2:176703/1‑150 (MQ=255)
                                                           tagacggCAATCAGTTTGCCCAGTTTCATCGAGCCTGCCAATGAGCAAGCTTCATGCGAGATGCCTTCCATTAAGCAGCCGTCGCCCAGGAAGGCGTAGGTGTAGTGGTCGACGATTTCGTGGCCGGGGCGGTTGAATTGGCCCGCCAGc                                                                                        <  1:139838/144‑1 (MQ=255)
                                                              acggCAATCAGTTTGCCCAGTTTCATCGAGCCTGCCAATGAGCAAGCTTCATGCGAGATGCCTTCCATTAAGCAGCCGTCGCCCAGGAAGGCGTAGGTGTAGTGGTCGACGATTTCGTGGCCGGGGCGGTTGAATTGGCCCGCCAGCgcg                                                                                     <  2:94201/147‑3 (MQ=255)
                                                               aTGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGttt                                                                                    >  2:99941/1‑150 (MQ=255)
                                                                ggCAATCAGTTTGCCCAGTTTCATCGAGCCTGCCAATGAGCAAGCTTCATGCGAGATGCCTTCCATTAAGCAGCCGTCGCCCAGGAAGGCGTAGGTGTAGTGGTCGACGATTTCGTGGCCGGGGCGGTTGAATTGGCCCGCCcgcgcgcg                                                                                   >  1:109079/2‑146 (MQ=255)
                                                                 gCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGttttt                                                                                  >  2:155328/1‑150 (MQ=255)
                                                                   aaTCAGTTTGCCCAGTTTCATCGAGCCTGCCAATGAGCAAGCTTCATGCGAGATGCCTTCCATTAAGCAGCCGTCGCCCAGGAAGGCGTAGGTGTAGTGGTCGACGATTTCGTGGCCGGGGCGGTTGAATTGGCCCGCCAGCgcgcgttc                                                                                <  2:39535/150‑8 (MQ=255)
                                                                   aaTCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTc                                                                                >  2:20401/1‑150 (MQ=255)
                                                                       aGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGc                                                                             <  2:190504/149‑1 (MQ=255)
                                                                         tttACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAAt                                                                          <  1:70251/150‑1 (MQ=255)
                                                                         tttACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAAt                                                                          >  2:11163/1‑150 (MQ=255)
                                                                           tACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATc                                                                         <  2:11936/149‑1 (MQ=255)
                                                                              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)
                                                                                                                                                   |                                                                                                                                                   
CGGTGTCGTCGGTGAACCAGCCTTCAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACC  >  CP009273/3074391‑3074685

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