Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,114,520 C→A 6.1% intergenic (‑168/+30) 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,5200CA 6.1% 320.7 / 9.3 114intergenic (‑168/+30)pppA/sslEbifunctional prepilin leader peptidase/ methylase/putative secreted and surface‑associated lipoprotein mucinase
Reads supporting (aligned to +/- strand):  ref base C (55/51);  new base A (3/4);  total (58/56)
Fisher's exact test for biased strand distribution p-value = 7.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.54e-01

Only 100 of 120 total aligned reads displayed.

TAACTTGCTGGATGCGGCGTAAAACGCCTTATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGCTCACACCGCCCTCGAA  >  NC_000913/3114372‑3114669
                                                                                                                                                    |                                                                                                                                                     
 aaCTTGCTGGCTGCGGCGTAAAACGCCTTATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATgcgg                                                                                                                                                     <  1:715555/150‑1 (MQ=255)
 aaCTTGCTGGATGCGGCGTAAAACGCCTTATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATgcgg                                                                                                                                                     >  1:1064052/1‑150 (MQ=255)
     tGCTGGATGCGGCGTAAAACGCCTTATCCGTCCTACGGGTGTCTGCGGGCGCAATTAATAGCGGCTTTCCCCGTGTGTTGGGGAGATTAGCAGGAGCATCGGGCAAATGTGTTCACCCGGTGCGCCTTATTTCATGCCGGATGCGGcgcg                                                                                                                                                 >  1:324461/1‑150 (MQ=255)
      gCTGGATGCGGCGTAAAACGCCTTATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGa                                                                                                                                                >  1:1529487/1‑150 (MQ=255)
       cTGGATGCGGCGTAAAACGCCTTATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAg                                                                                                                                               <  1:749103/150‑1 (MQ=255)
       cTGGATGCGGCGTAAAACGCCTTATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTCTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCAAGCCGGATGCGGCGCGAg                                                                                                                                               <  2:86317/150‑1 (MQ=255)
          gATGCGGCGTAAACCGCCTTATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTTTTTCATGCCGGATGCGGCGCGAgcgc                                                                                                                                            >  2:1463717/1‑150 (MQ=255)
          gATGCGGCGTAAAACGCCTTATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAgcgc                                                                                                                                            >  2:886341/1‑150 (MQ=255)
             gcggcgTAAAACGCCTTATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCtt                                                                                                                                         <  1:507948/150‑1 (MQ=255)
              cggcgTAAAACGCCTTATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTa                                                                                                                                        >  2:108183/1‑150 (MQ=255)
              cggcgTAAAACGCCTTATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTa                                                                                                                                        >  2:76210/1‑150 (MQ=255)
              cggcgTAAAACGCCTTATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCACGAGCGCCTTa                                                                                                                                        <  2:756329/150‑1 (MQ=255)
               ggcgTAAAACGCCTTATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTAt                                                                                                                                       >  2:1577023/1‑150 (MQ=255)
                gcgTAAAACGCCTTATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATc                                                                                                                                      >  2:47116/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)
                       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)
                             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:134696/150‑1 (MQ=255)
                                                      cgcAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCgg                                                                                                <  2:1708778/150‑1 (MQ=255)
                                                      cgcAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCgg                                                                                                >  2:936612/1‑150 (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                                                                                 <  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:1091485/1‑150 (MQ=255)
                                                                               cTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGCGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTcc                                                                       <  1:47116/150‑1 (MQ=255)
                                                                                 gTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTcccg                                                                     >  1:794666/1‑149 (MQ=255)
                                                                                   aGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTcc                                                                       >  2:1183065/1‑146 (MQ=255)
                                                                                   aGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTcc                                                                       <  1:1183065/146‑1 (MQ=255)
                                                                                   aGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGt                                                                   <  1:178732/150‑1 (MQ=255)
                                                                                    ggCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGtt                                                                  >  1:27810/1‑150 (MQ=255)
                                                                                        cGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCgg                                                              >  2:128923/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)
                                                                                               cAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGc                                                       <  1:1657916/150‑1 (MQ=255)
                                                                                                 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:980571/1‑150 (MQ=255)
                                                                                                          ggCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGa                                            >  1:1219965/1‑150 (MQ=255)
                                                                                                          ggCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGa                                            <  1:1452260/150‑1 (MQ=255)
                                                                                                            caaagccgtTACCCCGGCGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGAcg                                          <  2:1569173/143‑1 (MQ=255)
                                                                                                              aaTGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGAcgcg                                        >  1:782973/1‑150 (MQ=255)
                                                                                                                tgtgTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAg                                      >  2:1243477/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)
                                                                                                                       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)
                                                                                                                                 ccTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACt                     <  1:1352406/150‑1 (MQ=255)
                                                                                                                                 ccTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACt                     >  2:1172581/1‑150 (MQ=255)
                                                                                                                                 ccTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACt                     >  2:51699/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)
                                                                                                                                          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)
                                                                                                                                                  ggcggggGGAGCGCGTTATCCGGCTTAGGGGCTACCCGGGCCGACATCTTATCCTCGGAACCCTGATAAATGGTTTCGGGT‑CCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGCTCACACCGCCCTCGa   <  2:162170/149‑1 (MQ=255)
                                                                                                                                                    cggcgCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGCTCACACCGCCCTCGaa  >  1:364044/1‑150 (MQ=255)
                                                                                                                                                    cggcgCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGCTCACACAGCCCTCGac  >  1:186239/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TAACTTGCTGGATGCGGCGTAAAACGCCTTATCCGTCCTACGGGTGTCTGCCAGCGCAAATAATCGCGGCTTTCCCCGTCTGTAGGCCCGATAAGCAGGCGCATCGGGCAAATGTGTTAACCCGGTGCGCCTTATTTCATGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGCTTACTCGGCAGACATCTTATGCTCGGTAACCTGATTAATGGTTTCCGGTCCCTGTTCCGGTTTCGGCAGATCGAGTGACGCGAGCGTGTTGTAAGCCGACTGGCTCACACCGCCCTCGAA  >  NC_000913/3114372‑3114669

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