Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,114,534 T→C 5.2% intergenic (‑182/+16) pppA ← / ← sslE bifunctional prepilin leader peptidase/ methylase/putative secreted and surface‑associated lipoprotein mucinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,114,5340TC 5.2% 728.7 / 20.0 212intergenic (‑182/+16)pppA/sslEbifunctional prepilin leader peptidase/ methylase/putative secreted and surface‑associated lipoprotein mucinase
Reads supporting (aligned to +/- strand):  ref base T (100/100);  new base C (6/5);  total (107/105)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.75e-01

Only 100 of 213 total aligned reads displayed.

GCGTAAAACGCCTTATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGCTCACACCGCCCTCGAAGCTCATCTCGCTC  >  NC_000913/3114388‑3114682
                                                                                                                                                  |                                                                                                                                                    
gcgTAAAACGCCTTATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATc                                                                                                                                                   <  2:121830/150‑1 (MQ=255)
gcgTAAAACGCCTTATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATGTCATGCCGGATGCGGCGCGAGCGCCTTATc                                                                                                                                                   >  1:1633478/1‑150 (MQ=255)
 cgTAAAACGCCTTATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATcc                                                                                                                                                  >  1:545765/1‑150 (MQ=255)
   tAAAACGCCTTATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCgg                                                                                                                                                >  1:2818168/1‑150 (MQ=255)
   tAAAACGCCTTATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCgg                                                                                                                                                <  2:1724712/150‑1 (MQ=255)
   tAAAACGCCTTATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCgg                                                                                                                                                <  2:1823044/150‑1 (MQ=255)
    aaaaCGCCTTATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGc                                                                                                                                               <  2:462623/150‑1 (MQ=255)
    aaaaCGCCTTATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGc                                                                                                                                               <  1:4114/150‑1 (MQ=255)
      aaCGCCTTATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCt                                                                                                                                             >  1:327639/1‑150 (MQ=255)
      aaCGCCTTATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATAAGGCGCTCGCGCCGCATCCGGcat                                                                                                                                             <  2:1606212/150‑3 (MQ=255)
      aaCGCCTTATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATAAGGCGCTCGCGCCGCATCCGGcat                                                                                                                                             <  2:1255880/150‑3 (MQ=255)
       aCGCCTTATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTa                                                                                                                                            <  1:2688944/150‑1 (MQ=255)
        cGCCTTATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTAc                                                                                                                                           <  2:1500603/150‑1 (MQ=255)
          ccTTATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACgg                                                                                                                                         <  2:815649/150‑1 (MQ=255)
          ccTTATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACgg                                                                                                                                         <  1:2494839/150‑1 (MQ=255)
           cTTATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACggg                                                                                                                                        <  1:1575832/150‑1 (MQ=255)
           cTTATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACggg                                                                                                                                        <  1:1455312/150‑1 (MQ=255)
               tCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCtaa                                                                                                                                    >  2:2718061/1‑148 (MQ=255)
               tCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTa                                                                                                                                    <  2:508392/150‑1 (MQ=255)
               tCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTa                                                                                                                                    <  1:2205412/150‑1 (MQ=255)
                 cGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACt                                                                                                                                  >  2:1225696/1‑150 (MQ=255)
                  gTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTc                                                                                                                                 >  2:379472/1‑150 (MQ=255)
                    ccTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCgg                                                                                                                               <  2:359388/150‑1 (MQ=255)
                    ccTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCgg                                                                                                                               <  1:1170508/150‑1 (MQ=255)
                      tACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATCAGCAGGCGCATCGGGCAAATGTGGTAACCCGGTGCGCCTTAGGTCATGCCGGAGGCGGCGCGAGCGCCGGCGCCGGCCTACGGGCTGACTCGGCa                                                                                                                             >  1:1984485/1‑150 (MQ=255)
                      tACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCa                                                                                                                             <  1:1271037/150‑1 (MQ=255)
                            tgtCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATc                                                                                                                       >  2:2083559/1‑150 (MQ=255)
                               cTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTa                                                                                                                    >  1:2172342/1‑150 (MQ=255)
                                tGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTAt                                                                                                                   <  2:468183/150‑1 (MQ=255)
                                  ccAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGc                                                                                                                 >  1:2050061/1‑150 (MQ=255)
                                  ccAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGc                                                                                                                 >  1:1901447/1‑150 (MQ=255)
                                  ccAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGc                                                                                                                 >  2:624453/1‑150 (MQ=255)
                                    aGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTc                                                                                                               <  2:331924/150‑1 (MQ=255)
                                     gcgcAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCg                                                                                                              >  1:2824861/1‑150 (MQ=255)
                                     gcgcAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCg                                                                                                              >  2:197405/1‑150 (MQ=255)
                                          aataatCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAAc                                                                                                         <  2:2777970/150‑1 (MQ=255)
                                            tcatCGCGGCTTTCCCCGTCTGTAGGCCCGAGAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCt                                                                                                       <  2:2721864/148‑1 (MQ=255)
                                              atCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGa                                                                                                     <  2:1152413/150‑1 (MQ=255)
                                                  cgGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTaa                                                                                                 <  2:1519789/150‑1 (MQ=255)
                                                  cgGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTaa                                                                                                 <  2:1380133/150‑1 (MQ=255)
                                                  cgGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTaa                                                                                                 <  2:832320/150‑1 (MQ=255)
                                                    gCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATg                                                                                               <  2:929182/150‑1 (MQ=255)
                                                       ttCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGtt                                                                                            >  2:2749785/1‑150 (MQ=255)
                                                           ccGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCg                                                                                        >  1:1215204/1‑150 (MQ=255)
                                                             gTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGt                                                                                      >  2:2521245/1‑150 (MQ=255)
                                                              tCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTc                                                                                     <  2:2106566/150‑1 (MQ=255)
                                                                  tAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTg                                                                                 >  2:1657349/1‑150 (MQ=255)
                                                                     gCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTc                                                                              <  2:2818168/150‑1 (MQ=255)
                                                                     gCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTc                                                                              >  1:489342/1‑150 (MQ=255)
                                                                       ccGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCg                                                                            >  2:241252/1‑150 (MQ=255)
                                                                        cGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCgg                                                                           >  1:1080052/1‑150 (MQ=255)
                                                                         gATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGt                                                                          >  1:626510/1‑150 (MQ=255)
                                                                         gATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGt                                                                          >  2:1445041/1‑150 (MQ=255)
                                                                          aTAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGtt                                                                         >  2:1929042/1‑150 (MQ=255)
                                                                             aGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCg                                                                      >  1:1561557/1‑150 (MQ=255)
                                                                              gCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCgg                                                                     >  2:645013/1‑150 (MQ=255)
                                                                              gCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCgg                                                                     >  2:2917570/1‑150 (MQ=255)
                                                                              gCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCgg                                                                     >  2:2059903/1‑150 (MQ=255)
                                                                              gCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCgg                                                                     >  2:1059620/1‑150 (MQ=255)
                                                                                aGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTCTCGGCa                                                                   <  1:727680/150‑1 (MQ=255)
                                                                                   cgcATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGAt                                                                >  2:1649949/1‑150 (MQ=255)
                                                                                      aTCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGa                                                             <  1:833825/150‑1 (MQ=255)
                                                                                       tCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAg                                                            >  2:2275026/1‑150 (MQ=255)
                                                                                          ggCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGa                                                         <  2:743348/150‑1 (MQ=255)
                                                                                          ggCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGa                                                         <  2:1411454/150‑1 (MQ=255)
                                                                                          ggCAAATGTGTTAACCCGGTGCGCCTTATTTCAGGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGa                                                         <  1:358553/150‑1 (MQ=255)
                                                                                            cAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGAcg                                                       >  2:697064/1‑150 (MQ=255)
                                                                                             aaaTGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGAcgc                                                      <  2:2939751/150‑1 (MQ=255)
                                                                                             aaaTGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGAcgc                                                      <  2:2154815/150‑1 (MQ=255)
                                                                                                      aaCCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGtgtt                                             <  1:323858/150‑1 (MQ=255)
                                                                                                       aCCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGtgttg                                            >  1:506162/1‑150 (MQ=255)
                                                                                                         ccGGTGCGCCTTATGTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCGTACTCGGCAGACATCTTATGCTCGGTAACCTGAGTCATGGTTTCCGGTCCCTGTTCCGGGTTCGGCAGATCGAGTGACGCGAGCGTGTTGTa                                          >  2:552479/1‑150 (MQ=39)
                                                                                                          cGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTaa                                         >  2:1242357/1‑150 (MQ=255)
                                                                                                           ggTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAg                                        >  1:1130757/1‑150 (MQ=255)
                                                                                                           ggTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAg                                        >  1:1130662/1‑150 (MQ=255)
                                                                                                                gcCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGAc                                   >  1:2895485/1‑150 (MQ=255)
                                                                                                                 ccTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACt                                  <  1:468841/150‑1 (MQ=255)
                                                                                                                  cTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTg                                 >  1:358052/1‑150 (MQ=255)
                                                                                                                   ttATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTgg                                >  1:1648414/1‑150 (MQ=255)
                                                                                                                     aTTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGCt                              <  2:2111236/150‑1 (MQ=255)
                                                                                                                       ttCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGCTca                            <  2:2127494/150‑1 (MQ=255)
                                                                                                                       ttCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGCTca                            <  1:197405/150‑1 (MQ=255)
                                                                                                                            gtcGGATAAGGCGCTCGCGCCGCATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTGCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGCTCACACCg                       >  2:2178757/3‑150 (MQ=255)
                                                                                                                            gCCGGATAAGGCGCTCGCGCCGCATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGCTCACACCg                       >  1:896862/1‑150 (MQ=255)
                                                                                                                             ccGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGCTCACACCGc                      >  2:443066/1‑150 (MQ=255)
                                                                                                                              cGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGCTCACACCGcc                     >  1:1859315/1‑150 (MQ=255)
                                                                                                                                gataaggcgCTCGCGCCGCATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGCTCACACCGCCCt                   >  1:2400186/6‑150 (MQ=255)
                                                                                                                                 aTGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGCTCACACCGCCCTc                  <  1:2749785/150‑1 (MQ=255)
                                                                                                                                    cggcgCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGCTCACACCGCCCTCGaa               <  2:2824861/150‑1 (MQ=255)
                                                                                                                                     ggcgCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGCTCACACCGCCCTCGAAg              <  1:1957859/150‑1 (MQ=255)
                                                                                                                                     ggcgCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGCTCACACCGCCCTCGAAg              <  1:2417042/150‑1 (MQ=255)
                                                                                                                                     ggcgCGAGCGCCTTATCCGGCCTACGGGCTTACTCGACAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGCTCACACCGCCCTCGAAg              <  2:719428/150‑1 (MQ=255)
                                                                                                                                      gcgcGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGCTCACACCGCCCTCGAAGc             >  1:1060954/1‑150 (MQ=255)
                                                                                                                                        gcgTGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGCTCACACCGCCCTCGAAGCTc           <  1:2926079/150‑1 (MQ=255)
                                                                                                                                        gcgAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGCTCACACCGCCCTCGAAGCTc           <  1:1466891/150‑1 (MQ=255)
                                                                                                                                            gcgcCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGCTCACACCGCCCTCGAAGCTCAtct       <  1:1382429/150‑1 (MQ=255)
                                                                                                                                            gcgcCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGCTCACACCGCCCTCGAAGCTCAtct       <  1:1480678/150‑1 (MQ=255)
                                                                                                                                                cTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGCTCACACCGCCCTCGAAGCTCATctcgct   >  1:2871545/1‑150 (MQ=255)
                                                                                                                                                cTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGCTCACACCGCCCTCGAAGCTCATctcgct   >  2:796689/1‑150 (MQ=255)
                                                                                                                                                cTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGCTCACACCGCCCTCGAAGCTCATctcgct   >  1:1691376/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GCGTAAAACGCCTTATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGCTCACACCGCCCTCGAAGCTCATCTCGCTC  >  NC_000913/3114388‑3114682

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: ‑