Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 8,642 T→C 83.3% S135S (TCT→TCC talB → transaldolase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092738,6420TC83.3% 779.2 / 138.5 311S135S (TCT→TCCtalBtransaldolase B
Reads supporting (aligned to +/- strand):  ref base T (26/26);  new base C (145/114);  total (171/140)
Fisher's exact test for biased strand distribution p-value = 4.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.22e-01

Only 100 of 311 total aligned reads displayed.

AAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCAATTGCGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGAC  >  CP009273/8496‑8789
                                                                                                                                                  |                                                                                                                                                   
  aCTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCAATTGCGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTg                                                                                                                                                >  1:294261/1‑150 (MQ=255)
  aCTGATCCCGGGCCGTATCTCCACCGAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTg                                                                                                                                                <  2:941941/150‑1 (MQ=255)
   cTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCAATTGCGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTgg                                                                                                                                               >  1:195061/1‑150 (MQ=255)
     gatcccGGGCCGTATCTCCACCGAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTgg                                                                                                                                               >  1:507165/5‑148 (MQ=255)
     gatcccGGGCCGTATCTCCACCGAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTgg                                                                                                                                               >  2:456752/5‑148 (MQ=255)
       tcccGGGCCGTATCTCCACCGAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAgg                                                                                                                                           <  2:125522/148‑1 (MQ=255)
         ccGGGCCGTATCTCCACCGAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAggg                                                                                                                                          >  1:739142/1‑149 (MQ=255)
              ccGTATCTCCACCGAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCg                                                                                                                                    >  2:743267/1‑150 (MQ=255)
                  aTCTCCACCGAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCgcgcc                                                                                                                                >  2:686114/1‑149 (MQ=255)
                     tccaccgAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGcgccgct                                                                                                                             >  1:430962/7‑149 (MQ=255)
                         ccgAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGaa                                                                                                                          >  1:148727/3‑149 (MQ=255)
                            aaGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCAATTGCGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACa                                                                                                                        <  1:339137/148‑1 (MQ=255)
                             aGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACa                                                                                                                        <  2:661483/147‑1 (MQ=255)
                                  acgcacgTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGaa                                                                                                                 >  2:316990/6‑149 (MQ=255)
                                     cacgTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGaaa                                                                                                                >  1:565545/3‑147 (MQ=255)
                                        gTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAaagaag                                                                                                          <  1:517803/150‑1 (MQ=255)
                                        gTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAaagaag                                                                                                          <  2:212815/150‑1 (MQ=255)
                                             tCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGc                                                                                                        <  1:760840/147‑1 (MQ=255)
                                             tCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCAt                                                                                                      >  1:809444/1‑149 (MQ=255)
                                               cTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCAtta                                                                                                    <  2:790213/149‑3 (MQ=255)
                                               cTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCAtt                                                                                                     <  1:490675/148‑2 (MQ=255)
                                               cTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCAt                                                                                                      <  2:826803/147‑1 (MQ=255)
                                               cTATGACACCGAAGCGTCAATTGCGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCaa                                                                                                   >  1:179810/1‑150 (MQ=255)
                                                tATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAAc                                                                                                  >  1:107375/1‑150 (MQ=255)
                                                 aTGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACt                                                                                                 >  2:881370/1‑150 (MQ=255)
                                                  tGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACt                                                                                                 <  1:921723/149‑1 (MQ=255)
                                                  tGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTg                                                                                                >  1:1017891/1‑150 (MQ=255)
                                                  tGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTg                                                                                                >  1:493621/1‑150 (MQ=255)
                                                   gACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTg                                                                                                >  2:700187/1‑149 (MQ=255)
                                                    acacCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGt                                                                                               <  2:723384/149‑1 (MQ=255)
                                                    acacCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGt                                                                                               >  1:2010/1‑149 (MQ=255)
                                                      acCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAAc                                                                                            >  2:242342/1‑150 (MQ=255)
                                                       ccGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAAcc                                                                                           <  1:845437/150‑1 (MQ=255)
                                                        cGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCt                                                                                          <  1:975681/150‑1 (MQ=255)
                                                         gAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCt                                                                                          >  1:30672/1‑149 (MQ=255)
                                                         gAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCt                                                                                          >  2:633745/1‑149 (MQ=255)
                                                         gAAGCGTCAATTGCGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCt                                                                                          >  2:1005181/1‑149 (MQ=255)
                                                             cGTCAATTGCGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGAccc                                                                                     <  2:197901/150‑1 (MQ=255)
                                                              gTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGct                                                                                    <  2:772962/150‑1 (MQ=255)
                                                                ctatcgCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGctgc                                                                                  >  1:93706/6‑150 (MQ=255)
                                                                ctatcgCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGctgc                                                                                  >  1:915497/6‑150 (MQ=255)
                                                                  aTTGCGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCctgctg                                                                                >  2:521775/1‑150 (MQ=255)
                                                                    cgCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGCAACCTGACGCTGCTGtt                                                                              >  2:832131/2‑150 (MQ=9)
                                                                      cGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTtct                                                                            >  1:481802/1‑150 (MQ=255)
                                                                      cGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTtct                                                                            >  2:25864/1‑150 (MQ=255)
                                                                      cGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTtct                                                                            >  1:718990/1‑150 (MQ=255)
                                                                       gAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTtctc                                                                           >  2:294093/1‑150 (MQ=255)
                                                                       gAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTtctc                                                                           <  1:116643/150‑1 (MQ=255)
                                                                       gAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCACCAACTGTAACCTGACCCTGCTGTtctc                                                                           >  2:1068004/1‑150 (MQ=255)
                                                                          aGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCtt                                                                        <  2:1046948/150‑1 (MQ=255)
                                                                                 cgcatcaTCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAg                                                                 >  2:272678/7‑150 (MQ=255)
                                                                                 cgcatcaTCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAg                                                                 <  2:227187/144‑1 (MQ=255)
                                                                                  gcatcaTCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAgg                                                                >  2:1056015/6‑150 (MQ=255)
                                                                                   catcaTCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGc                                                               >  1:259018/5‑150 (MQ=255)
                                                                                    atcaTCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCt                                                              <  1:884351/147‑1 (MQ=255)
                                                                                     tcaTCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTc                                                             <  2:75502/148‑1 (MQ=255)
                                                                                     tGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTc                                                             >  1:880374/1‑150 (MQ=255)
                                                                                       aTCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCg                                                            <  2:609633/149‑1 (MQ=255)
                                                                                           agctgtACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCtt                                                       >  1:1052052/6‑150 (MQ=255)
                                                                                             ctgtACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTtgc                                                     <  2:141636/147‑2 (MQ=255)
                                                                                             ctgtACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTtgc                                                     <  2:885197/147‑2 (MQ=255)
                                                                                               ctACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGc                                                   >  2:324217/1‑150 (MQ=255)
                                                                                                 aCAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTgcgcc                                                  >  1:416574/1‑148 (MQ=255)
                                                                                                  cAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCgg                                                 >  1:947691/1‑149 (MQ=255)
                                                                                                  cAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGa                                                <  2:190701/150‑1 (MQ=255)
                                                                                                  cAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTgcgccg                                                 >  1:861519/1‑147 (MQ=255)
                                                                                                          cTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCg                                        >  2:465559/1‑150 (MQ=255)
                                                                                                           gggcatcaGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCgt                                       <  2:84662/143‑1 (MQ=255)
                                                                                                             gcatcaGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCgtct                                     <  1:142013/145‑3 (MQ=255)
                                                                                                              catcaGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCtt                                    >  1:242453/5‑150 (MQ=255)
                                                                                                               atcaGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTc                                   >  1:675866/4‑150 (MQ=255)
                                                                                                               atcaGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTc                                   >  1:439249/4‑150 (MQ=255)
                                                                                                               atcaGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTc                                   >  1:782061/4‑150 (MQ=255)
                                                                                                               atcaGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTc                                   >  2:37420/4‑150 (MQ=255)
                                                                                                                tcaGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTcc                                  >  2:585217/3‑150 (MQ=255)
                                                                                                                  aGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCgg                                           >  1:991194/1‑139 (MQ=255)
                                                                                                                   gCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGa                               >  1:692873/1‑150 (MQ=255)
                                                                                                                    cAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGAt                              <  2:565545/150‑1 (MQ=255)
                                                                                                                     aaCGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGAtc                             >  2:868038/1‑150 (MQ=255)
                                                                                                                       cGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGAtttc                           <  2:652025/150‑1 (MQ=255)
                                                                                                                        gATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCg                          >  1:64106/1‑150 (MQ=255)
                                                                                                                        gATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCg                          >  1:58795/1‑150 (MQ=255)
                                                                                                                         aTCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGc                         >  2:926084/1‑150 (MQ=255)
                                                                                                                          tCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGcc                        <  1:212756/150‑1 (MQ=255)
                                                                                                                           cgcaTTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCg                       >  2:541546/4‑150 (MQ=255)
                                                                                                                           cgcaTTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGAAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCg                       >  2:102105/4‑150 (MQ=255)
                                                                                                                            gcaTTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGt                      >  1:461798/3‑150 (MQ=255)
                                                                                                                              aTTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGttc                    >  2:408478/1‑149 (MQ=255)
                                                                                                                                tCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTtcgt                  >  2:365035/1‑150 (MQ=255)
                                                                                                                                tCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTtcgt                  >  1:227772/1‑150 (MQ=255)
                                                                                                                                   gATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTgg               >  1:219918/1‑150 (MQ=255)
                                                                                                                                    aTCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTgg               <  1:133776/149‑1 (MQ=255)
                                                                                                                                    aTCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGCTCGTTGGc              >  2:966612/1‑150 (MQ=255)
                                                                                                                                     tCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGcc             >  2:502590/1‑150 (MQ=255)
                                                                                                                                      cAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCg            >  2:949163/1‑150 (MQ=255)
                                                                                                                                        aaCTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTa          >  1:753484/1‑150 (MQ=255)
                                                                                                                                            ggCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCt      >  1:376397/1‑150 (MQ=255)
                                                                                                                                            ggCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCt      >  2:903134/1‑150 (MQ=255)
                                                                                                                                              cTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTg    <  1:681837/150‑1 (MQ=255)
                                                                                                                                                 cTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGAc  <  2:42212/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
AAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCAATTGCGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGAC  >  CP009273/8496‑8789

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