Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,114,526 G→T 7.0% intergenic (‑174/+24) 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,5260GT 7.0% 356.7 / 8.1 117intergenic (‑174/+24)pppA/sslEbifunctional prepilin leader peptidase/ methylase/putative secreted and surface‑associated lipoprotein mucinase
Reads supporting (aligned to +/- strand):  ref base G (55/52);  new base T (3/6);  total (59/58)
Fisher's exact test for biased strand distribution p-value = 4.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 123 total aligned reads displayed.

TGCTGGATGCGGCGTAAAACGCCTTATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGCTCACACCGCCCTCGAAGCTCAT  >  NC_000913/3114377‑3114675
                                                                                                                                                     |                                                                                                                                                     
tGCTGGATGCGGCGTAAAACGCCTTATCCGTCCTACGGGTGTCTGCGGGCGCAATTAATAGCGGCTTTCCCCGTGTGTTGGGGAGATTAGCAGGAGCATCGGGCAAATGTGTTCACCCGGTGCGCCTTATTTCATGCCGGATGCGGcgcg                                                                                                                                                       >  1:324461/1‑150 (MQ=255)
 gCTGGATGCGGCGTAAAACGCCTTATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCTGGAGCATCGGGCAAATGTGTTACCCCGGTGCGCCTTATTTCATGCCGGATGCGGcgcgc                                                                                                                                                      >  2:869480/1‑149 (MQ=255)
 gCTGGATGCGGCGTAAAACGCCTTATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGa                                                                                                                                                      >  1:1529487/1‑150 (MQ=255)
  cTGGATGCGGCGTAAAACGCCTTATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAg                                                                                                                                                     <  1:749103/150‑1 (MQ=255)
  cTGGATGCGGCGTAAAACGCCTTATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTCTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCAAGCCGGATGCGGCGCGAg                                                                                                                                                     <  2:86317/150‑1 (MQ=255)
     gATGCGGCGTAAAACGCCTTATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAgcgc                                                                                                                                                  >  2:886341/1‑150 (MQ=255)
        gcggcgTAAAACGCCTTATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCtt                                                                                                                                               <  1:507948/150‑1 (MQ=255)
         cggcgTAAAACGCCTTATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTa                                                                                                                                              >  2:76210/1‑150 (MQ=255)
         cggcgTAAAACGCCTTATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTa                                                                                                                                              >  2:108183/1‑150 (MQ=255)
          ggcgTAAAACGCCTTATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTAt                                                                                                                                             >  2:1577023/1‑150 (MQ=255)
           gcgTAAAACGCCTTATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATc                                                                                                                                            >  2:47116/1‑150 (MQ=255)
             gTAAAACGCCTTATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCg                                                                                                                                          >  2:178732/1‑150 (MQ=255)
              tAAAACGCCTTATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCgg                                                                                                                                         <  1:1199512/150‑1 (MQ=255)
                aaaCGCCTTATCCGTCCTACGGGTGTCTGCCAGAGCAAATAATCGCTGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCATATGTGTTAACCCGGTGCGCCTTATTTCATACCGGATGCGGCGAGAGCGCATTATACGGcc                                                                                                                                       >  1:1113264/1‑150 (MQ=255)
                 aaCGCCTTATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATAAGGCGCTCGCGTCGCATCCGGcat                                                                                                                                      <  2:1035846/150‑3 (MQ=255)
                  ccGCCTTCTCCGTCCTCCGGGTGTCTGCCCGCGCAAATCATCGCGGCTCTCCCCGTCTGTCGGCCCGATCAGCCGGCGCAGCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATAAGGCGCTCGCGCCGCATCCGGcatg                                                                                                                                     <  1:324600/149‑4 (MQ=255)
                    gCCTTATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACg                                                                                                                                   <  2:1399258/150‑1 (MQ=255)
                      cTTATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACggg                                                                                                                                 <  2:1529487/150‑1 (MQ=255)
                      cTTATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACggg                                                                                                                                 >  1:1377435/1‑150 (MQ=255)
                       ttATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGc                                                                                                                                <  2:1223175/150‑1 (MQ=255)
                        tATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCt                                                                                                                               <  2:1275064/150‑1 (MQ=255)
                        tATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCt                                                                                                                               <  1:1009072/150‑1 (MQ=255)
                            cGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACt                                                                                                                           >  2:629383/1‑150 (MQ=255)
                            cGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACt                                                                                                                           >  1:552948/1‑150 (MQ=255)
                                cTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGc                                                                                                                       <  1:163379/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:134696/150‑1 (MQ=255)
                                                 cgcAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCgg                                                                                                      >  2:936612/1‑150 (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                                                                                       <  1:327832/150‑1 (MQ=255)
                                                                cTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATgg                                                                                       >  2:471532/1‑150 (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)
                                                                          cTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGCGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTcc                                                                             <  1:47116/150‑1 (MQ=255)
                                                                            gTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTcccg                                                                           >  1:794666/1‑149 (MQ=255)
                                                                              aGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTcc                                                                             <  1:1183065/146‑1 (MQ=255)
                                                                              aGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTcc                                                                             >  2:1183065/1‑146 (MQ=255)
                                                                              aGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGt                                                                         <  1:178732/150‑1 (MQ=255)
                                                                                gCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGACTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGTCAGACATCTTATGCTCGGTACGCGGATTAATGGTTTCCGGTCCCTGtta                                                                       >  2:314777/1‑149 (MQ=39)
                                                                                   cGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCgg                                                                    >  2:128923/1‑150 (MQ=255)
                                                                                    gATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCTGGTCCGGtt                                                                  >  1:162170/1‑150 (MQ=255)
                                                                                     aTAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGtt                                                                  <  1:629383/150‑1 (MQ=255)
                                                                                        aGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCg                                                               >  1:599206/1‑150 (MQ=255)
                                                                                        aGCAGGCGCATCGGGCAAATGGGGTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACGCGGCAGACATCTTATGCTCGGTAACGTGAGTAATGGTTTCCGGTCCCTGTTCCGGTTTCg                                                               <  2:911881/150‑1 (MQ=39)
                                                                                            ggCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAg                                                           >  2:673430/1‑150 (MQ=255)
                                                                                             gcgcATCGGGCAAATGTGTTAACCGGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGa                                                          <  2:1380385/150‑1 (MQ=255)
                                                                                                  tCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAg                                                     >  1:1694667/1‑150 (MQ=255)
                                                                                                   cGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGACTTATCCGGCCTTCGGGCTTACTCGGCAGACCTCTTATGCTCGGGAACCTGATTAGTGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAgg                                                    >  2:82234/1‑149 (MQ=255)
                                                                                                    gggCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACCGCTTCTGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTg                                                   >  1:1530931/1‑150 (MQ=255)
                                                                                                    gggCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGGTTCGGCAGATCGAGTg                                                   >  2:279946/1‑150 (MQ=255)
                                                                                                     ggCAAATGTGTTTACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGa                                                  <  1:87704/150‑1 (MQ=255)
                                                                                                     ggCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGa                                                  <  1:1452260/150‑1 (MQ=255)
                                                                                                     ggCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGa                                                  >  1:980571/1‑150 (MQ=255)
                                                                                                     ggCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGa                                                  >  1:1219965/1‑150 (MQ=255)
                                                                                                       caaagccgtTACCCCGGCGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGAcg                                                <  2:1569173/143‑1 (MQ=255)
                                                                                                         aaTGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGAcgcg                                              >  1:326375/1‑150 (MQ=255)
                                                                                                         aaTGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGAcgcg                                              >  1:782973/1‑150 (MQ=255)
                                                                                                                taccgcGGGGCGGCTTATTTCATGCCGGATGCGGCGGTAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCgtgt                                       <  1:139008/145‑1 (MQ=255)
                                                                                                                tAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTTCTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGATTGACGCGAGCgtgt                                       >  2:1058548/1‑150 (MQ=255)
                                                                                                                  aCCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGtgttg                                     <  1:1027518/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                           <  1:1352406/150‑1 (MQ=255)
                                                                                                                            ccTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACt                           >  2:1172581/1‑150 (MQ=255)
                                                                                                                               tATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGc                        <  2:1404817/150‑1 (MQ=255)
                                                                                                                               tATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGc                        <  2:1694667/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)
                                                                                                                                     aTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGCTCacac                  <  1:673430/150‑1 (MQ=255)
                                                                                                                                      ggCCGGATAAGGCGCTCGCGCCGCATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGTAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGCTCacac                  >  2:144200/2‑149 (MQ=255)
                                                                                                                                      ggCCGGATAAGGCGCTCGCGCCGCATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGTAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGCTCacac                  <  1:144200/148‑1 (MQ=255)
                                                                                                                                        ccGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGCTCACACCGc               >  1:1704997/1‑150 (MQ=255)
                                                                                                                                        ccGGATAAGGCGCTCGCGCCGCATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGCTCACACCGc               >  1:1658536/1‑150 (MQ=255)
                                                                                                                                          ggATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGCTCACACCGccc             <  1:1652849/150‑1 (MQ=255)
                                                                                                                                               cggcgCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGCTCACACCGCCCTCGaa        >  1:364044/1‑150 (MQ=255)
                                                                                                                                               cggcgCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGCTCACACAGCCCTCGac        >  1:186239/1‑149 (MQ=255)
                                                                                                                                                ggcgCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGCTCACACCGCCCTCGAAg       >  1:939974/1‑150 (MQ=255)
                                                                                                                                                ggcgCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGCTCACACCGCCCTCGAAg       <  2:1064052/150‑1 (MQ=255)
                                                                                                                                                  cgcgAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGCTCACACCGCCCTCGAAGCt     >  1:109475/1‑150 (MQ=255)
                                                                                                                                                     gAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGCTCACACCGCCCTCGAAGCTCAt  <  1:1474570/150‑1 (MQ=255)
                                                                                                                                                     gAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGCTCACACCGCCCTCGAAGCTCAt  >  1:378191/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGCTGGATGCGGCGTAAAACGCCTTATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGCTCACACCGCCCTCGAAGCTCAT  >  NC_000913/3114377‑3114675

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: