Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,373,985 C→A 84.5% intergenic (+37/‑117) yhcB → / → degQ DUF1043 family inner membrane‑anchored protein/serine endoprotease, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,373,9850CA84.5% 782.4 / 169.2 336intergenic (+37/‑117)yhcB/degQDUF1043 family inner membrane‑anchored protein/serine endoprotease, periplasmic
Reads supporting (aligned to +/- strand):  ref base C (18/34);  new base A (116/168);  total (134/202)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.94e-02

Only 100 of 350 total aligned reads displayed.

AATCGTCTGGCAGAGTCTGAAGCCAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTCTCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGAAAAAACAAACCCAGCTGTTGAGTGCATT  >  CP009273/3373838‑3374133
                                                                                                                                                   |                                                                                                                                                    
aaTCGTCTGGCAGAGTCTGAAGCCAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCtcttct                                                                                                                                                     >  1:4093315/1‑149 (MQ=255)
 aTCGTCTGGCAGAGTCTGAAGCCAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTtatc                                                                                                                                                    >  2:3046065/1‑149 (MQ=255)
        ggCAGAGTCTGAAGCCAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATctccctcc                                                                                                                                             <  1:2486300/149‑1 (MQ=255)
         gCAGAGTCTGAAGCCAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCtccctccc                                                                                                                                            >  1:3216321/1‑149 (MQ=255)
             agTCTGAAGCCAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTCTCTCCCTCCCCGAc                                                                                                                                        <  1:119566/149‑1 (MQ=255)
               tCTGAAGCCAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTa                                                                                                                                      <  1:3975244/149‑1 (MQ=255)
                 tGAAGCCAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATc                                                                                                                                    >  2:4331694/1‑149 (MQ=255)
                  gAAGCCAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCa                                                                                                                                   <  1:779956/149‑1 (MQ=255)
                  gAAGCCAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCa                                                                                                                                   >  1:4683028/1‑149 (MQ=255)
                   aaGCCAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCAt                                                                                                                                  <  2:3645246/149‑1 (MQ=255)
                    aGCCAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCAtt                                                                                                                                 <  2:671616/149‑1 (MQ=255)
                    aGCCAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCAtt                                                                                                                                 >  2:4587977/1‑149 (MQ=255)
                     gCCAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCAttt                                                                                                                                >  2:2802631/1‑149 (MQ=255)
                      ccAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTa                                                                                                                               <  2:189464/149‑1 (MQ=255)
                       cAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTaa                                                                                                                              <  2:1947539/149‑1 (MQ=255)
                                tCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGtc                                                                                                                     <  1:4638317/149‑1 (MQ=255)
                                   ggCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCa                                                                                                                  >  2:474374/1‑149 (MQ=255)
                                   ggCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCa                                                                                                                  >  2:208096/1‑149 (MQ=255)
                                     cACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCAtt                                                                                                                <  1:4399362/149‑1 (MQ=255)
                                      aCCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTg                                                                                                               >  2:1330454/1‑149 (MQ=255)
                                      aCCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTg                                                                                                               >  2:2701894/1‑149 (MQ=255)
                                      aCCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTg                                                                                                               <  2:3165073/149‑1 (MQ=255)
                                       ccGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTCTCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGt                                                                                                              <  2:504632/149‑1 (MQ=255)
                                       ccGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGt                                                                                                              <  2:445643/149‑1 (MQ=255)
                                        cGGTGCAGATGCCTCGCGTCTATTCTGAAGGCTCATCCGGGCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCTTTGtt                                                                                                             <  2:1510475/149‑1 (MQ=255)
                                         ggTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTCTCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTa                                                                                                            <  1:3571387/149‑1 (MQ=255)
                                          gTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAg                                                                                                           <  1:3292914/149‑1 (MQ=255)
                                          gTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAg                                                                                                           <  1:4447815/149‑1 (MQ=255)
                                          gTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAg                                                                                                           >  1:2214290/1‑149 (MQ=255)
                                           tGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGc                                                                                                          <  1:3014753/149‑1 (MQ=255)
                                            gCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGcc                                                                                                         >  1:3162685/1‑149 (MQ=255)
                                            gCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACGAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTCTCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGcc                                                                                                         <  2:1139088/149‑1 (MQ=255)
                                               gATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTc                                                                                                      <  2:2038458/149‑1 (MQ=255)
                                                aTGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCt                                                                                                     <  1:3049881/149‑1 (MQ=255)
                                                     tCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGaaaa                                                                                                >  2:264494/1‑149 (MQ=255)
                                                        cgACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTc                                                                                             <  2:2120480/149‑1 (MQ=255)
                                                          aCTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTCTCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCaa                                                                                           >  2:3329315/1‑149 (MQ=255)
                                                          aCTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTCTCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCaa                                                                                           >  1:2313167/1‑149 (MQ=255)
                                                           cTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAAt                                                                                          <  2:1674343/149‑1 (MQ=255)
                                                            tATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATa                                                                                         <  2:282950/149‑1 (MQ=255)
                                                            tATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATa                                                                                         <  2:815842/149‑1 (MQ=255)
                                                             aTTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATaa                                                                                        <  1:3790633/149‑1 (MQ=255)
                                                               tCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACa                                                                                      <  1:1185410/149‑1 (MQ=255)
                                                                cTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACAt                                                                                     <  2:3968478/149‑1 (MQ=255)
                                                                   aaGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCACTAACATCaa                                                                                  >  2:4094994/1‑149 (MQ=255)
                                                                   aaGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCaa                                                                                  >  2:251613/1‑149 (MQ=255)
                                                                    aGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCaaa                                                                                 <  2:2882215/149‑1 (MQ=255)
                                                                    aGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCaaa                                                                                 >  2:348625/1‑149 (MQ=255)
                                                                      gcgcATCCGGCCTGCTGCGTACTGGCGCGACGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACt                                                                               >  1:3642764/1‑149 (MQ=255)
                                                                        gcATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTCTCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGt                                                                             <  2:2333384/149‑1 (MQ=255)
                                                                          aTCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGttt                                                                           >  1:159175/1‑149 (MQ=255)
                                                                           tCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGtttt                                                                          <  2:1274626/149‑1 (MQ=255)
                                                                                ccTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTCTCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAAtc                                                                     <  1:993834/149‑1 (MQ=255)
                                                                                 ctgctgCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAAtct                                                                    <  1:258564/149‑1 (MQ=255)
                                                                                 ctgctgCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAAtct                                                                    <  2:592477/149‑1 (MQ=255)
                                                                                  tgctgcGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAAtctc                                                                   <  1:4339527/149‑1 (MQ=255)
                                                                                  tgctgcGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAAtctc                                                                   <  1:4339427/149‑1 (MQ=255)
                                                                                  tgctgcGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAAtctc                                                                   <  1:4110154/149‑1 (MQ=255)
                                                                                   gctgcGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATctct                                                                  <  1:4017527/149‑1 (MQ=255)
                                                                                   gctgcGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATctct                                                                  <  2:4397776/149‑1 (MQ=255)
                                                                                    ctgcGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCtt                                                                 <  1:297377/149‑1 (MQ=255)
                                                                                             ggCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTCTCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTATGAATCTCTTTTCTTATCa                                                        >  1:1485077/1‑149 (MQ=255)
                                                                                              gcgcGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCAt                                                       <  2:214016/149‑1 (MQ=255)
                                                                                               cgcgAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCAtt                                                      <  1:191838/149‑1 (MQ=255)
                                                                                               cgcgAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCAtt                                                      <  1:2095068/149‑1 (MQ=255)
                                                                                                  gAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTCTCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAg                                                   <  2:4011543/149‑1 (MQ=255)
                                                                                                  gAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAg                                                   <  1:564267/149‑1 (MQ=255)
                                                                                                  gAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAg                                                   >  2:4211763/1‑149 (MQ=255)
                                                                                                  gAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAg                                                   >  1:604123/1‑149 (MQ=255)
                                                                                                    aGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGt                                                 >  1:1246612/1‑149 (MQ=255)
                                                                                                    aGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGt                                                 >  2:3129888/1‑149 (MQ=255)
                                                                                                     gcgcGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTCTCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTa                                                <  2:680022/149‑1 (MQ=255)
                                                                                                       gcgACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACg                                              >  1:2376657/1‑149 (MQ=255)
                                                                                                        cgACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACga                                             <  1:3340308/149‑1 (MQ=255)
                                                                                                         gACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACgag                                            <  2:895012/149‑1 (MQ=255)
                                                                                                         gACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACgag                                            <  2:2395614/149‑1 (MQ=255)
                                                                                                          aCTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTCTCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACgaga                                           >  1:4247912/1‑149 (MQ=255)
                                                                                                           cTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGagag                                          >  1:4150236/1‑149 (MQ=255)
                                                                                                              atttatttTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAg                                       >  2:1593315/1‑149 (MQ=255)
                                                                                                              atttatttTTCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAg                                       <  1:1522348/149‑1 (MQ=255)
                                                                                                                  atttTTCGGGCGCAGCCATTGCGCCCTCCTCTTCTCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGaat                                   <  1:3460303/149‑1 (MQ=255)
                                                                                                                    ttttCGGGCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGaataa                                 <  2:73550/149‑1 (MQ=255)
                                                                                                                         gggCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGaaa                            <  1:1305443/149‑1 (MQ=255)
                                                                                                                          ggCGCAGCCATTGCGCCCTCCTCTTCTCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGaaaa                           >  1:2909154/1‑149 (MQ=255)
                                                                                                                          ggCGCAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGaaaa                           <  1:4257652/149‑1 (MQ=255)
                                                                                                                           gcgcAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGaaaaa                          <  2:3902602/149‑1 (MQ=255)
                                                                                                                            cgcAGCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGaaaaaa                         <  1:2670188/149‑1 (MQ=255)
                                                                                                                                gCCATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGAAAaaacaaa                     <  2:3692933/149‑1 (MQ=255)
                                                                                                                                 ccATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGAAAaaacaaac                    >  2:4681158/1‑149 (MQ=255)
                                                                                                                                 ccATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGAAAaaacaaac                    <  2:381195/149‑1 (MQ=255)
                                                                                                                                 ccATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGAAAaaacaaac                    <  1:2735153/149‑1 (MQ=255)
                                                                                                                                  cATTGCGCCCTCCTCTTCTCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGAAAAAACAAAcc                   >  1:1686712/1‑149 (MQ=255)
                                                                                                                                  cATTGCGCCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGAAAAAACTAAcc                   >  2:4131473/1‑149 (MQ=255)
                                                                                                                                       cgcCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGAAAAAACAAACCCAGCt              <  1:2812646/149‑1 (MQ=255)
                                                                                                                                        gcCCTCCTCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGAAAAAACAAACCCAGCTg             >  2:1541087/1‑149 (MQ=255)
                                                                                                                                          cctcctCTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGAAAAAACAAACCCAGCTGtt           >  1:2443276/1‑149 (MQ=255)
                                                                                                                                           ctcctcTTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGAAAAAACAAACCCAGCTGTTg          >  2:2204952/1‑149 (MQ=255)
                                                                                                                                            tcctcTTATCTCCCTCCCCGACTATCATTTATTCTGGTTTCTCATTGTTAGCAGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGAAAAAACAAACCCAGCTGTTGa         <  1:4038659/149‑1 (MQ=255)
                                                                                                                                              ctctTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGAAAAAACAAACCCAGCTGTTGAGt       <  1:2253608/149‑1 (MQ=255)
                                                                                                                                              ctctTATCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGAAAAAACAAACCCAGCTGTTGAGt       <  1:4632760/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
AATCGTCTGGCAGAGTCTGAAGCCAGCAACGATCAGGCACCGGTGCAGATGCCTCGCGACTATTCTGAAGGCGCATCCGGCCTGCTGCGTACTGGCGCGAAGCGCGACTAATTTATTTTTCGGGCGCAGCCATTGCGCCCTCCTCTTCTCTCCCTCCCCGACTATCATTTAATCTGGTGTCTCATTGTTAGCCGTCTGAAAATTCAATAACATCAAACTGTTTTGAATCTCTTTTCTTATCATTCAGGTACGAGAGCAGGAATAATGAAAAAACAAACCCAGCTGTTGAGTGCATT  >  CP009273/3373838‑3374133

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: