Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,074,591 G→T 100% A135E (GCG→GAG)  tktA ← transketolase 1, thiamine triphosphate‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,074,5910GT98.8% 538.8 / ‑2.8 168A135E (GCG→GAG) tktAtransketolase 1, thiamine triphosphate‑binding
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base T (89/77);  total (89/79)
Fisher's exact test for biased strand distribution p-value = 2.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 171 total aligned reads displayed.

ATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTT  >  CP009273/3074446‑3074733
                                                                                                                                                 |                                                                                                                                              
  cGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGctccgc                                                                                                                                            <  2:181811/148‑1 (MQ=255)
       aaTGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCTCCGCCAGCg                                                                                                                                       <  2:1641021/148‑1 (MQ=255)
       aaTGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCg                                                                                                                                       <  1:1990242/148‑1 (MQ=255)
        aTGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCTCCGCCAGCGt                                                                                                                                      >  2:192453/1‑148 (MQ=255)
        aTGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGGAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCTCCGCCAGCGt                                                                                                                                      >  1:3206446/1‑148 (MQ=255)
          gCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCTCCGCCAGCGttt                                                                                                                                    <  2:1549336/148‑1 (MQ=255)
           cAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCTCCGCCAGCGtttt                                                                                                                                   >  2:901559/1‑148 (MQ=255)
             aTCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCTCCGCCAGCGtttttt                                                                                                                                 >  2:1250460/1‑148 (MQ=255)
             aTCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCTCCGCCAGCGtttttt                                                                                                                                 <  2:2494552/148‑1 (MQ=255)
             aTCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCTCCGCCAGCGtttttt                                                                                                                                 <  2:1593056/148‑1 (MQ=255)
              tCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCTCCGCCAGCGTTTTTTc                                                                                                                                >  1:391499/1‑148 (MQ=255)
              tCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCTCCGCCAGCGTTTTTTc                                                                                                                                >  2:2357755/1‑148 (MQ=255)
               cAGTTTACCCAGCTTCAGAGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCAGTCTCCCATGAATGCGTAGGTGTAGTGGTCGACAATGTCGTGTCCCGGACGGTTAATCTGCTCCGCCAGCGTTTTTTCt                                                                                                                               >  2:1209397/1‑148 (MQ=255)
                    tACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCTCCGCCAGCGTTTTTTCTGCAAt                                                                                                                          <  1:92416/148‑1 (MQ=255)
                         aGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCTCCGCCAGCGTTTTTTCTGCAATCGCCa                                                                                                                     >  2:1523276/1‑148 (MQ=255)
                           cTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCTCCGCCAGCGTTTTTTCTGCAATCGCCATa                                                                                                                   <  2:147704/148‑1 (MQ=255)
                            ttCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCTCCGCCAGCGTTTTTTCTGCAATCGCCATAc                                                                                                                  <  2:1631206/148‑1 (MQ=255)
                             tCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCTCCGCCAGCGTTTTTTCTGCAATCGCCATAcc                                                                                                                 >  1:181272/1‑148 (MQ=255)
                                 cGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCTCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACt                                                                                                             >  1:1637346/1‑148 (MQ=255)
                                 cGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCTCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACt                                                                                                             >  2:2077659/1‑148 (MQ=255)
                                    aCCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCTCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCg                                                                                                          >  2:1719954/1‑148 (MQ=255)
                                    aCCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCTCCGCCAGCGTTGTTTCTGCAATCGCCATACCGACTGCg                                                                                                          >  2:1017827/1‑148 (MQ=255)
                                     cccGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCTCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGt                                                                                                         >  1:76088/1‑148 (MQ=255)
                                     cccGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCTCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGt                                                                                                         >  1:2336415/1‑148 (MQ=255)
                                         ccAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCTCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGc                                                                                                     >  1:1464394/1‑148 (MQ=255)
                                           agagagCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCTCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCaa                                                                                                   >  1:79578/1‑148 (MQ=255)
                                             agagCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCTCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATa                                                                                                 >  1:2584410/1‑148 (MQ=255)
                                             agagCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCTCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATa                                                                                                 >  2:2499006/1‑148 (MQ=255)
                                             agagCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCTCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATa                                                                                                 >  1:1597887/1‑148 (MQ=255)
                                             agagCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCTCCGCCAGCGTTTTTTATGCAATCGCAAAAACGACTTCGATGGCaaaa                                                                                                 >  1:2073792/1‑146 (MQ=255)
                                              gagCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCTCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATAc                                                                                                >  1:719590/1‑148 (MQ=255)
                                               agCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCTCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATAcc                                                                                               <  1:2065094/148‑1 (MQ=255)
                                                 cAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCTCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCt                                                                                             >  2:2131341/1‑148 (MQ=255)
                                                 cAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCTCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCt                                                                                             >  1:1342436/1‑148 (MQ=255)
                                                    aCTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCTCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGAc                                                                                          >  1:2834209/1‑148 (MQ=255)
                                                     cTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCTCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGAcc                                                                                         <  1:2984703/148‑1 (MQ=255)
                                                     cTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCTCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGAcc                                                                                         <  1:1219054/148‑1 (MQ=255)
                                                       tCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCTCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCa                                                                                       >  2:3160917/1‑148 (MQ=255)
                                                            ggAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCTCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGa                                                                                  <  2:647897/148‑1 (MQ=255)
                                                            ggAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCTCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGa                                                                                  >  2:2251163/1‑148 (MQ=255)
                                                            ggAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCTCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGa                                                                                  >  1:108189/1‑148 (MQ=255)
                                                            ggAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCTCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGGTGGCAATACCCTGACCCAGCGGa                                                                                  >  1:1965885/1‑148 (MQ=255)
                                                             gagaTGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCTCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGAc                                                                                 <  2:2394674/148‑1 (MQ=255)
                                                               gaTGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCTCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCg                                                                               <  2:684375/148‑1 (MQ=255)
                                                                  gCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCTCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCggtg                                                                            <  1:331287/148‑1 (MQ=255)
                                                                   ccTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCTCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCggtgg                                                                           <  1:1363201/148‑1 (MQ=255)
                                                                   ccTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCTCCGCCAGCGTTGTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCggtgg                                                                           <  2:2770757/148‑1 (MQ=255)
                                                                    cTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCTCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCggtggt                                                                          <  1:3515310/148‑1 (MQ=255)
                                                                      tCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCTCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCggtggtgg                                                                        <  2:2434576/148‑1 (MQ=255)
                                                                       ccATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCTCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCggtggtggt                                                                       <  2:425319/148‑1 (MQ=255)
                                                                           catGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCTCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTcc                                                                   >  1:1982548/1‑148 (MQ=255)
                                                                           catGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCTCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTcc                                                                   >  2:3392284/1‑148 (MQ=255)
                                                                                 gCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCTCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCa                                                             >  2:1445816/1‑148 (MQ=255)
                                                                                        cccATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCTCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGgtgt                                                      >  2:215572/1‑148 (MQ=255)
                                                                                          cATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCTCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTaa                                                    <  1:2357755/148‑1 (MQ=255)
                                                                                           aTGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCTCCCCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAAc                                                   >  1:2647409/1‑148 (MQ=255)
                                                                                            tGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCTCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAAcc                                                  >  2:2048749/1‑148 (MQ=255)
                                                                                             gAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCTCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGAACCAGCGGACCGGTGGTGTTTTCCACACCAGCGGTGTAAccc                                                 >  1:1892916/1‑148 (MQ=255)
                                                                                              aaGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCTCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCa                                                <  1:264546/148‑1 (MQ=255)
                                                                                              aaGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCTCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCa                                                <  1:2036594/148‑1 (MQ=255)
                                                                                               aGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCTCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCAc                                               >  2:2994221/1‑148 (MQ=255)
                                                                                                 gCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCTCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACtt                                             <  1:2443299/148‑1 (MQ=255)
                                                                                                     aGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCTCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCgg                                         <  2:66292/148‑1 (MQ=255)
                                                                                                     aGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCTCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCgg                                         <  2:1843903/148‑1 (MQ=255)
                                                                                                     aGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCTCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCgg                                         <  1:759773/148‑1 (MQ=255)
                                                                                                      ggTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCTCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCggg                                        >  2:47049/1‑148 (MQ=255)
                                                                                                       gtgtAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCTCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGt                                       <  1:1075054/148‑1 (MQ=255)
                                                                                                           aGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCTCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGAcc                                   <  2:3206446/148‑1 (MQ=255)
                                                                                                           aGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCTCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGAcc                                   >  1:1745207/1‑148 (MQ=255)
                                                                                                           aGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCTCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGAcc                                   >  1:3483499/1‑148 (MQ=255)
                                                                                                           aGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCTCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGAcc                                   <  1:623867/148‑1 (MQ=255)
                                                                                                               gTCGACAATGTCGTGGCCCGGACGGTTAAACTGCTCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGa                               <  2:391499/148‑1 (MQ=255)
                                                                                                               gTCGACAATGTCGTGGCCCGGACGGTTAAACTGCTCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGa                               <  2:2407424/148‑1 (MQ=255)
                                                                                                                 cGACAATGTCGTGGCCCGGACGGTTAAACTGCTCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGt                             <  1:1021282/148‑1 (MQ=255)
                                                                                                                 cGACAATGTCGTGGCCCGGACGGTTAAACTGCTCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGt                             <  1:3164310/148‑1 (MQ=255)
                                                                                                                 cGACAATGTCGTGGCCCGGACGGTTAAACTGCTCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGt                             <  2:1357721/148‑1 (MQ=255)
                                                                                                                 cGACAATGTCGTGGCCCGGACGGTTAAACTGCTCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGt                             <  1:1021030/148‑1 (MQ=255)
                                                                                                                  gACAATGTCGTGGCCCGGACGGTTAAACTGCTCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGtt                            <  1:3243207/148‑1 (MQ=255)
                                                                                                                  gACAATGTCGTGGCCCGGACGGTTAAACTGCTCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCAACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGtt                            <  1:2214052/148‑1 (MQ=255)
                                                                                                                   aCAATGTCGTGGCCCGGACGGTTAAACTGCTCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGttt                           >  2:2846842/1‑148 (MQ=255)
                                                                                                                    cAATGTCGTGGCCCGGACGGTTAAACTGCTCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGtttt                          >  1:502676/1‑148 (MQ=255)
                                                                                                                      aTGTCGTGGCCCGGACGGTTAAACTGCTCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTag                        <  2:2374042/148‑1 (MQ=255)
                                                                                                                        gTCGTGGCCCGGACGGTTAAACTGCTCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTagag                      <  1:3489951/148‑1 (MQ=255)
                                                                                                                            tGGCCCGGACGGTTAAACTGCTCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCa                  <  2:2600890/148‑1 (MQ=255)
                                                                                                                            tGGCCCGGACGGTTAAACTGCTCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCa                  <  2:2031385/148‑1 (MQ=255)
                                                                                                                             ggCCCGGACGGTTAAACTGCTCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAg                 <  2:194561/148‑1 (MQ=255)
                                                                                                                               cccGGACGGTTAAACTGCTCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCt               >  2:1635255/1‑148 (MQ=255)
                                                                                                                                 ccgACGGCGAAACCGCTCCGCCAGCGTTTTGTCTGCAATCGCCATACCGCCTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGa             <  2:3498866/146‑1 (MQ=255)
                                                                                                                                 cGGACGGTTAAACTGCTCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGa             >  1:1100470/1‑148 (MQ=255)
                                                                                                                                     cGGTTAAACTGCTCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGa         <  2:872352/148‑1 (MQ=255)
                                                                                                                                     cGGTTAAACTGCTCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGa         <  1:2077659/148‑1 (MQ=255)
                                                                                                                                      ggTTAAACTGCTCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGaa        >  2:2015274/1‑148 (MQ=255)
                                                                                                                                      ggTTAAACTGCTCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGaa        <  2:2340817/148‑1 (MQ=255)
                                                                                                                                       gTTAAACTGCTCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAg       <  2:2589073/148‑1 (MQ=255)
                                                                                                                                       gTTAAACGGCTCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAg       <  2:1891759/148‑1 (MQ=255)
                                                                                                                                        ttAAACTGCTCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGt      >  1:2581919/1‑148 (MQ=255)
                                                                                                                                        ttAAACTGCTCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGt      >  2:1785449/1‑148 (MQ=255)
                                                                                                                                          aaaCTGCTCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGttt    >  1:2116752/1‑148 (MQ=255)
                                                                                                                                            aCTGCTCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGttttt  >  1:193446/1‑148 (MQ=255)
                                                                                                                                            aCTGCTCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGttttt  >  2:3420096/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
ATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTT  >  CP009273/3074446‑3074733

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: