Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,114,533 T→G 5.9% intergenic (‑181/+17) 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,5330TG 5.9% 338.9 / 9.8 119intergenic (‑181/+17)pppA/sslEbifunctional prepilin leader peptidase/ methylase/putative secreted and surface‑associated lipoprotein mucinase
Reads supporting (aligned to +/- strand):  ref base T (58/54);  new base G (3/4);  total (61/58)
Fisher's exact test for biased strand distribution p-value = 7.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.19e-01

Only 100 of 121 total aligned reads displayed.

GCGGCGTAAAACGCCTTATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGCTCACACCGCCCTCGAAGCTCATCTCGCTC  >  NC_000913/3114385‑3114682
                                                                                                                                                    |                                                                                                                                                     
gcggcgTAAAACGCCTTATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCtt                                                                                                                                                      <  1:507948/150‑1 (MQ=255)
 cggcgTAAAACGCCTTATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTa                                                                                                                                                     >  2:76210/1‑150 (MQ=255)
 cggcgTAAAACGCCTTATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTa                                                                                                                                                     >  2:108183/1‑150 (MQ=255)
 cggcgTAAAACGCCTTATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCACGAGCGCCTTa                                                                                                                                                     <  2:756329/150‑1 (MQ=255)
  ggcgTAAAACGCCTTATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTAt                                                                                                                                                    >  2:1577023/1‑150 (MQ=255)
    cgTAAAACGCCTTATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGAAGGACCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATcc                                                                                                                                                  >  1:67978/1‑150 (MQ=255)
     gTAAAACGCCTTATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCg                                                                                                                                                 >  2:178732/1‑150 (MQ=255)
      tAAAACGCCTTATCCTTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATAAGGCGCTCGCGCCGCATCCgg                                                                                                                                                <  1:340931/150‑1 (MQ=255)
      tAAAACGCCTTATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCgg                                                                                                                                                <  1:1199512/150‑1 (MQ=255)
        aaaCGCCTTATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGcc                                                                                                                                              <  2:883383/150‑1 (MQ=255)
        aaaCGCCTTATCCGTCCTACGGGTGTCTGCCAGAGCAAATAATCGCTGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCATATGTGTTAACCCGGTGCGCCTTATTTCATACCGGATGCGGCGAGAGCGCATTATACGGcc                                                                                                                                              >  1:1113264/1‑150 (MQ=255)
         aaCGCCTTATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATAAGGCGCTCGCGTCGCATCCGGcat                                                                                                                                             <  2:1035846/150‑3 (MQ=255)
            gCCTTATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACg                                                                                                                                          <  2:1399258/150‑1 (MQ=255)
              cTTATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACggg                                                                                                                                        >  1:1377435/1‑150 (MQ=255)
               ttATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGc                                                                                                                                       <  2:1223175/150‑1 (MQ=255)
                tATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCt                                                                                                                                      <  1:1009072/150‑1 (MQ=255)
                tATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCt                                                                                                                                      <  2:1275064/150‑1 (MQ=255)
                    cGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACt                                                                                                                                  >  1:552948/1‑150 (MQ=255)
                    cGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACt                                                                                                                                  >  2:629383/1‑150 (MQ=255)
                          aCGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAg                                                                                                                            <  1:1447900/150‑1 (MQ=255)
                           cGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCTGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGa                                                                                                                           <  2:707686/150‑1 (MQ=255)
                             ggTGTCTGCCAGCGCAAATAATCTCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTCTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACa                                                                                                                         <  2:40343/150‑1 (MQ=255)
                             ggTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACa                                                                                                                         >  2:1452260/1‑150 (MQ=255)
                                gtCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCt                                                                                                                      <  2:1385781/150‑1 (MQ=255)
                                      cAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACTGGCTTACTCGGCAGACATCTTATGCt                                                                                                                >  2:1365337/1‑150 (MQ=255)
                                       agggcAAATAATCGCTGCTTTCCCCTTCTGTAGTCCCGATAAGCAGGCGCATCTGTCAAATGTTTTAACCCGGTGCTCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGATc                                                                                                               <  2:1113264/147‑1 (MQ=38)
                                       aGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTc                                                                                                               >  2:1027518/1‑150 (MQ=255)
                                         cgcAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGCGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCgg                                                                                                             >  2:1651670/1‑150 (MQ=255)
                                         cgcAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCgg                                                                                                             <  2:1708778/150‑1 (MQ=255)
                                         cgcAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCgg                                                                                                             >  2:936612/1‑150 (MQ=255)
                                         cgcAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCgg                                                                                                             <  2:134696/150‑1 (MQ=255)
                                          gcAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGt                                                                                                            <  2:539889/150‑1 (MQ=255)
                                               taatCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCt                                                                                                       <  1:1320496/150‑1 (MQ=255)
                                                  tCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGAt                                                                                                    <  1:1599052/150‑1 (MQ=255)
                                                    gcgGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTa                                                                                                  >  2:1352406/1‑150 (MQ=255)
                                                     cgGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACTGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTaa                                                                                                 <  1:76210/150‑1 (MQ=255)
                                                     cgGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTaa                                                                                                 >  1:1150774/1‑150 (MQ=255)
                                                     cgGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTaa                                                                                                 <  1:108183/150‑1 (MQ=255)
                                                        cTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATgg                                                                                              >  2:471532/1‑150 (MQ=255)
                                                        cTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATgg                                                                                              <  1:327832/150‑1 (MQ=255)
                                                        cTTTCCCCGTCTGTAAGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATgg                                                                                              >  2:89948/1‑150 (MQ=255)
                                                             cccGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTcc                                                                                         >  1:1354432/1‑150 (MQ=255)
                                                             cccGGCTGTCGGCCCGCTAAGCAGGCGCATAGGTCAAATGTGTTAACCCGGGGCGACTTTTTTCATGCGGGTTGGGGGGATAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTcc                                                                                         <  2:324461/150‑1 (MQ=255)
                                                               cGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCgg                                                                                       >  2:696598/1‑150 (MQ=255)
                                                                      aGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTcc                                                                                    <  1:1183065/146‑1 (MQ=255)
                                                                      aGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTcc                                                                                    >  2:1183065/1‑146 (MQ=255)
                                                                      aGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGt                                                                                <  1:178732/150‑1 (MQ=255)
                                                                       ggCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGtt                                                                               >  1:27810/1‑150 (MQ=255)
                                                                        gCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGACTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGTCAGACATCTTATGCTCGGTACGCGGATTAATGGTTTCCGGTCCCTGtta                                                                              >  2:314777/1‑149 (MQ=39)
                                                                           cGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCgg                                                                           >  2:128923/1‑150 (MQ=255)
                                                                            gATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCTGGTCCGGtt                                                                         >  1:162170/1‑150 (MQ=255)
                                                                                aGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCg                                                                      <  1:829112/150‑1 (MQ=255)
                                                                                aGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCg                                                                      >  1:599206/1‑150 (MQ=255)
                                                                                aGCAGGCGCATCGGGCAAATGGGGTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACGCGGCAGACATCTTATGCTCGGTAACGTGAGTAATGGTTTCCGGTCCCTGTTCCGGTTTCg                                                                      <  2:911881/150‑1 (MQ=39)
                                                                                  cAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGc                                                                    <  1:1657916/150‑1 (MQ=255)
                                                                                     gcgcATCGGGCAAATGTGTTAACCGGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGa                                                                 <  2:1380385/150‑1 (MQ=255)
                                                                                          tCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAg                                                            >  1:1694667/1‑150 (MQ=255)
                                                                                            gggCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACCGCTTCTGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTg                                                          >  1:1530931/1‑150 (MQ=255)
                                                                                             ggCAAATGTGTTTACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGa                                                         <  1:87704/150‑1 (MQ=255)
                                                                                             ggCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGa                                                         >  1:1219965/1‑150 (MQ=255)
                                                                                             ggCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGa                                                         >  1:980571/1‑150 (MQ=255)
                                                                                               caaagccgtTACCCCGGCGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGAcg                                                       <  2:1569173/143‑1 (MQ=255)
                                                                                                 aaTGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGAcgcg                                                     >  1:782973/1‑150 (MQ=255)
                                                                                                 aaTGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGAcgcg                                                     >  1:326375/1‑150 (MQ=255)
                                                                                                     tgtTAACCCGGTGCGCCTTATTTCATGCCGGATTCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCg                                                 <  2:235592/150‑1 (MQ=255)
                                                                                                     tgtTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCg                                                 <  1:1577023/150‑1 (MQ=255)
                                                                                                        taccgcGGGGCGGCTTATTTCATGCCGGATGCGGCGGTAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCgtgt                                              <  1:139008/145‑1 (MQ=255)
                                                                                                        tAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTTCTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGATTGACGCGAGCgtgt                                              >  2:1058548/1‑150 (MQ=255)
                                                                                                          aCCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGtgttg                                            <  1:1027518/150‑1 (MQ=255)
                                                                                                           cccGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGtgttgt                                           <  2:38663/150‑1 (MQ=255)
                                                                                                            ccGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTa                                          <  2:552948/150‑1 (MQ=255)
                                                                                                                ggcgcCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGcc                                      <  1:936612/149‑1 (MQ=255)
                                                                                                                  cgcCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGAGTGTTGTAAGCCgc                                    >  2:1266296/1‑149 (MQ=255)
                                                                                                                    gcTTATTTCATTCCGGATGCGGCGCGAGCGCCTTATCCGGCCGACGGGCTTTCTGGGCAGACATGTTATGGTGGGTTACCTGATTATTGGTTTGCGGTGCCTGTTGCGGGTTCGGGAGATCGAGTGACGCGCGCGTGTTGTAAGCCGGCt                                  >  2:373952/2‑150 (MQ=255)
                                                                                                                    ccTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGCCt                                  >  2:1465935/1‑150 (MQ=255)
                                                                                                                    ccTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACt                                  >  2:51699/1‑150 (MQ=255)
                                                                                                                    ccTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACt                                  >  2:1172581/1‑150 (MQ=255)
                                                                                                                       tATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGc                               <  2:1694667/150‑1 (MQ=255)
                                                                                                                       tATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGc                               <  2:1404817/150‑1 (MQ=255)
                                                                                                                        cccgtaggCCGGATAAGGCGCTCGCGCCGCATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGCt                              >  2:1248404/8‑150 (MQ=38)
                                                                                                                         tttCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGCTc                             <  1:471532/150‑1 (MQ=255)
                                                                                                                           tCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCGGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGCTcac                           >  2:1667367/1‑150 (MQ=255)
                                                                                                                           tCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGCTcac                           >  1:499865/1‑150 (MQ=255)
                                                                                                                             aTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGCTCacac                         <  1:673430/150‑1 (MQ=255)
                                                                                                                              ggCCGGATAAGGCGCTCGCGCCGCATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGTAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGCTCacac                         >  2:144200/2‑149 (MQ=255)
                                                                                                                              ggCCGGATAAGGCGCTCGCGCCGCATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGTAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGCTCacac                         <  1:144200/148‑1 (MQ=255)
                                                                                                                                ccGGATAAGGCGCTCGCGCCGCATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGCTCACACCGc                      >  1:1658536/1‑150 (MQ=255)
                                                                                                                                  ggATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGCTCACACCGccc                    <  1:1652849/150‑1 (MQ=255)
                                                                                                                                     ggcggggGGAGCGCGTTATCCGGCTTAGGGGCTACCCGGGCCGACATCTTATCCTCGGAACCCTGATAAATGGTTTCGGGTCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGCTCACACCGCCCTCGa                <  2:162170/149‑1 (MQ=255)
                                                                                                                                       cggcgCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGCTCACACCGCCCTCGaa               >  1:364044/1‑150 (MQ=255)
                                                                                                                                       cggcgCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGCTCACACAGCCCTCGac               >  1:186239/1‑149 (MQ=255)
                                                                                                                                        ggcgCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGCTCACACCGCCCTCGAAg              <  2:1064052/150‑1 (MQ=255)
                                                                                                                                        ggcgCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGCTCACACCGCCCTCGAAg              >  1:939974/1‑150 (MQ=255)
                                                                                                                                          cgcgAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGCTCACACCGCCCTCGAAGCt            >  1:109475/1‑150 (MQ=255)
                                                                                                                                            ctaGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGCTCACACCGCCCTCGAAGCTCa          <  2:67978/148‑1 (MQ=255)
                                                                                                                                             gAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGCTCACACCGCCCTCGAAGCTCAt         >  1:378191/1‑150 (MQ=255)
                                                                                                                                             gAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGCTCACACCGCCCTCGAAGCTCAt         <  1:1474570/150‑1 (MQ=255)
                                                                                                                                                 gcCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGCTCACACCGCCCTCGAAGCTCATctcg     >  1:266461/1‑150 (MQ=255)
                                                                                                                                                   cTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGCTCACACCGCCCTCGAAGCTCATctcgct   <  2:1219965/150‑1 (MQ=255)
                                                                                                                                                    ttATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGCTCACACCGCCCTCGAAGCTCATctcgctc  <  1:279946/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCGGCGTAAAACGCCTTATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGCTCACACCGCCCTCGAAGCTCATCTCGCTC  >  NC_000913/3114385‑3114682

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