Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,573,360 C→A 5.3% A122A (GCC→GCA tktB → transketolase 2, thiamine triphosphate‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,573,3600CA 5.3% 859.3 / 16.1 223A122A (GCC→GCAtktBtransketolase 2, thiamine triphosphate‑binding
Reads supporting (aligned to +/- strand):  ref base C (124/87);  new base A (8/4);  total (132/91)
Fisher's exact test for biased strand distribution p-value = 7.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 231 total aligned reads displayed.

GCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGGCTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTGGCAGGCACG  >  CP009273/2573213‑2573507
                                                                                                                                                   |                                                                                                                                                   
 cTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGGCTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCg                                                                                                                                                    >  2:2649853/1‑148 (MQ=255)
 cTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGGCTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCg                                                                                                                                                    <  2:2003424/148‑1 (MQ=255)
   aCATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGGCTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTc                                                                                                                                                  >  1:2298294/1‑148 (MQ=255)
       cTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGGCTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGc                                                                                                                                              >  1:805175/1‑148 (MQ=255)
         gACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGGCTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTg                                                                                                                                            <  2:407375/148‑1 (MQ=255)
          aCCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGGCTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTgg                                                                                                                                           <  1:24059/148‑1 (MQ=255)
          aCCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGGCTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTgg                                                                                                                                           >  2:2990463/1‑148 (MQ=255)
          aCCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGGCTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTgg                                                                                                                                           <  1:2002550/148‑1 (MQ=255)
             ggTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGGCTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGa                                                                                                                                        >  2:2836409/1‑148 (MQ=255)
              gTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGGCTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGAt                                                                                                                                       >  1:2569013/1‑148 (MQ=255)
              gTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGGCTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGAt                                                                                                                                       <  1:1187162/148‑1 (MQ=255)
              gTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGGCTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGAt                                                                                                                                       <  2:1665719/148‑1 (MQ=255)
               ttACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGGCTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATa                                                                                                                                      >  2:3419216/1‑148 (MQ=255)
                 aCGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGGCTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGc                                                                                                                                    >  1:1149894/1‑148 (MQ=255)
                     ccTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGGCTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCagag                                                                                                                                <  1:2851655/148‑1 (MQ=255)
                     ccTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGGCTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCagag                                                                                                                                <  1:1886239/148‑1 (MQ=255)
                       tGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGGCTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCg                                                                                                                              >  2:1551077/1‑148 (MQ=255)
                         ccGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGGCTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTa                                                                                                                            <  2:2278724/148‑1 (MQ=255)
                         ccGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGGCTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTa                                                                                                                            <  1:944950/148‑1 (MQ=255)
                         ccGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGGCTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTa                                                                                                                            <  1:944975/148‑1 (MQ=255)
                         ccGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGGCTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTa                                                                                                                            >  2:3324908/1‑148 (MQ=255)
                         ccGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGGCTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTa                                                                                                                            >  2:1307564/1‑148 (MQ=255)
                           gCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGGCTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTaca                                                                                                                          <  1:2460260/148‑1 (MQ=255)
                              ggAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGGCTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTg                                                                                                                       >  2:2136335/1‑148 (MQ=255)
                               gaagaactgaagaactTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGGCTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTgg                                                                                                                      <  2:1100843/148‑1 (MQ=255)
                                aagaactgaagaactTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGGCTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTggc                                                                                                                     >  2:485093/1‑148 (MQ=255)
                                 agaactgaagaactTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGGCTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTaca                                                                                                                          <  2:246240/142‑1 (MQ=255)
                                  gaactgaagaactTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGGCTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTggcgg                                                                                                                   >  1:2563088/1‑148 (MQ=255)
                                   aactgaagaactTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGGCTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTggcggc                                                                                                                  >  1:2682546/1‑148 (MQ=255)
                                    actgaagaactTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGGCTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGgcggcg                                                                                                                 >  1:3432596/1‑148 (MQ=255)
                                        aagaactTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGGCTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGt                                                                                                             >  2:841080/1‑148 (MQ=255)
                                        aagaactTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGGCTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGt                                                                                                             >  1:130528/1‑148 (MQ=255)
                                        aagaactTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGGCTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGt                                                                                                             >  1:164078/1‑148 (MQ=255)
                                         agaactTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGGCTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGtt                                                                                                            >  2:1632011/1‑148 (MQ=255)
                                           aactTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGGCTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTa                                                                                                          >  2:1711135/1‑148 (MQ=255)
                                           aactTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGGCTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTa                                                                                                          >  2:1929384/1‑148 (MQ=255)
                                            actTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGGCTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTaa                                                                                                         <  1:1687564/148‑1 (MQ=255)
                                              ttCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGGCTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAAcc                                                                                                       <  2:118807/148‑1 (MQ=255)
                                               tCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGGCTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAAcca                                                                                                      >  1:173967/1‑148 (MQ=255)
                                                 cGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGGCTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAAccagc                                                                                                    <  2:1129085/148‑1 (MQ=255)
                                                 cGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGGCTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAAccagc                                                                                                    <  1:1469865/148‑1 (MQ=255)
                                                  gTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGGCTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAAccagcc                                                                                                   <  1:1025011/148‑1 (MQ=255)
                                                       ttGCATTCGAAAACCCCAGGCCACCCGGAGATTGGCTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGAcc                                                                                              >  2:2460981/1‑148 (MQ=255)
                                                        tGCATTCGAAAACCCCAGGCCACCCGGAGATTGGCTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCa                                                                                             >  1:1265010/1‑148 (MQ=255)
                                                          cATTCGAAAACCCCAGGCCACCCGGAGATTGGCTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATg                                                                                           >  2:1453736/1‑148 (MQ=255)
                                                            ttCGAAAACCCCAGGCCACCCGGAGATTGGCTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATgag                                                                                         >  2:756806/1‑148 (MQ=255)
                                                            ttCGAAAACCCCAGGCCACCCGGAGATTGGCTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATgag                                                                                         <  2:2507263/148‑1 (MQ=255)
                                                             tCGAAAACCCCAGGCCACCCGGAGATTGGCTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATgaga                                                                                        <  1:1775192/148‑1 (MQ=255)
                                                              cGAAAACCCCAGGCCACCCGGAGATTGGCTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGAt                                                                                       <  1:3389330/148‑1 (MQ=255)
                                                               gAAAACCCCAGGCCACCCGGAGATTGGCTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGAtc                                                                                      >  2:2300367/1‑148 (MQ=255)
                                                                aaaaCCCCAGGCCACCCGGAGATTGGCTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGAtcg                                                                                     <  2:3253160/148‑1 (MQ=255)
                                                                aaaaCCCCAGGCCACCCGGAGATTGGCTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGAtcg                                                                                     >  1:2616950/1‑148 (MQ=255)
                                                                 aaaCCCCAGGCCACCCGGAGATTGGCTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGAtcgt                                                                                    >  2:704693/1‑148 (MQ=255)
                                                                   aCCCCAGGCCACCCGGAGATTGGCTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGAtcgtcg                                                                                  >  2:3074207/1‑148 (MQ=255)
                                                                   aCCCCAGGCCACCCGGAGATTGGCTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGAtcgtcg                                                                                  >  1:217528/1‑148 (MQ=255)
                                                                    ccccAGGCCACCCGGAGATTGGCTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGa                                                                                 >  2:661770/1‑148 (MQ=255)
                                                                      ccAGGCCACCCGGAGATTGGCTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATc                                                                               >  1:3530335/1‑148 (MQ=255)
                                                                       cAGGCCACCCGGAGATTGGCTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCa                                                                              <  2:2754973/148‑1 (MQ=255)
                                                                        aGGCCACCCGGAGATTGGCTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCAc                                                                             >  1:917565/1‑148 (MQ=255)
                                                                        aGGCCACCCGGAGATTGGCTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCAc                                                                             >  1:3254808/1‑148 (MQ=255)
                                                                        aGGCCACCCGGAGATTGGCTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCAc                                                                             >  1:1487991/1‑148 (MQ=255)
                                                                        aGGCCACCCGGAGATTGGCTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCAGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCAc                                                                             <  2:3515456/148‑1 (MQ=255)
                                                                         ggCCACCCGGAGATTGGCTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACt                                                                            >  1:1160292/1‑148 (MQ=255)
                                                                          gCCACCCGGAGATTGGCTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACtt                                                                           >  2:3376041/1‑148 (MQ=255)
                                                                             aCCCGGAGATTGGCTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCAc                                                                        <  2:1817156/148‑1 (MQ=255)
                                                                                 ggAGATTGGCTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTAt                                                                    <  1:2208173/148‑1 (MQ=255)
                                                                                    gaTTGGCTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTAtgtg                                                                 >  2:2799843/1‑148 (MQ=255)
                                                                                     aTTGGCTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATgtgt                                                                <  2:1284339/148‑1 (MQ=255)
                                                                                       tGGCTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGttt                                                              >  2:357637/1‑148 (MQ=255)
                                                                                       tGGCTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGttt                                                              >  2:3485779/1‑148 (MQ=255)
                                                                                       tGGCTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCAGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGttt                                                              <  1:1597876/148‑1 (MQ=255)
                                                                                        ggCTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGTTTa                                                             <  2:3227814/148‑1 (MQ=255)
                                                                                         gCTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGTTTAt                                                            <  1:1926735/148‑1 (MQ=255)
                                                                                         gCTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGTTTAt                                                            >  1:2251767/1‑148 (MQ=255)
                                                                                          cTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGTTTATg                                                           <  2:889772/148‑1 (MQ=255)
                                                                                               cGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGa                                                      >  1:753575/1‑148 (MQ=255)
                                                                                                 ccAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACg                                                    <  1:3104449/148‑1 (MQ=255)
                                                                                                  cAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACgg                                                   <  1:689364/148‑1 (MQ=255)
                                                                                                        ttGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGTCGACGGCTGCCt                                             <  2:2876735/148‑1 (MQ=255)
                                                                                                        ttGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCt                                             >  1:2516039/1‑148 (MQ=255)
                                                                                                           aaaCCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGAt                                          <  2:1923898/148‑1 (MQ=255)
                                                                                                           aaaCCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGAt                                          <  1:2528/148‑1 (MQ=255)
                                                                                                                accaccGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAg                                     <  2:2811574/148‑1 (MQ=255)
                                                                                                                accaccGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAg                                     >  1:567637/1‑148 (MQ=255)
                                                                                                                 ccaccGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAgg                                    >  2:708513/1‑148 (MQ=255)
                                                                                                                     cGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTAtt                                >  1:1477989/1‑148 (MQ=255)
                                                                                                                          cGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCa                           <  1:3145088/148‑1 (MQ=255)
                                                                                                                               ggACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAg                      <  1:1386654/148‑1 (MQ=255)
                                                                                                                               ggACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAg                      <  2:2266738/148‑1 (MQ=255)
                                                                                                                                   aaGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTg                  >  2:2862251/1‑148 (MQ=255)
                                                                                                                                   aaGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTg                  >  2:358346/1‑148 (MQ=255)
                                                                                                                                     ggTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGtt                >  1:2864732/1‑148 (MQ=255)
                                                                                                                                       tttGGCGAACGCAGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCg              >  1:2747884/1‑148 (MQ=255)
                                                                                                                                       tttGGCGAACGCAGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCg              >  1:1319201/1‑148 (MQ=255)
                                                                                                                                           gCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTgg          <  2:3250363/148‑1 (MQ=255)
                                                                                                                                                cGCAGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTggcaggc     >  1:1246066/1‑148 (MQ=255)
                                                                                                                                                   cGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTGGCAGGCACg  <  1:209058/148‑1 (MQ=255)
                                                                                                                                                   cGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTGGCAGGCACg  <  2:1487991/148‑1 (MQ=255)
                                                                                                                                                   cGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTGGCAGGCACg  <  2:3254808/148‑1 (MQ=255)
                                                                                                                                                   cGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTGGCAGGCACg  <  2:917606/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGGCTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTGGCAGGCACG  >  CP009273/2573213‑2573507

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: