Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,574,0510GT40.2% 166.2 / 46.1 102V384V (GTC→GTApqqLputative periplasmic M16 family zinc metalloendopeptidase
Reads supporting (aligned to +/- strand):  ref base G (8/53);  new base T (41/0);  total (49/53)
Fisher's exact test for biased strand distribution p-value = 7.96e-21
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 106 total aligned reads displayed.

CAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGTGCAATTTTAACGCTGCGC  >  NC_000913/1573901‑1574201
                                                                                                                                                      |                                                                                                                                                        
caGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTg                                                                                                                                                          <  2:214153/151‑1 (MQ=255)
      tttCGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGGAGATTTTAcat                                                                                                                                                      >  1:1942796/1‑149 (MQ=255)
       ttCGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGAcatca                                                                                                                                                    <  1:1461337/150‑1 (MQ=255)
       ttCGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGAcatca                                                                                                                                                    <  1:1425944/150‑1 (MQ=255)
       ttCGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGGGGATTTTAcatc                                                                                                                                                     >  1:257760/1‑149 (MQ=255)
        tCGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGAcatcat                                                                                                                                                   <  2:805822/150‑1 (MQ=255)
         cGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACatcatc                                                                                                                                                  <  2:1100512/150‑1 (MQ=255)
          gAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGa                                                                                                                                                <  2:1735775/151‑1 (MQ=255)
            aaGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGt                                                                                                                                              <  1:251662/151‑1 (MQ=255)
                tGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGGAGATTTGACATCATCGAGttttt                                                                                                                                          >  1:1167253/1‑151 (MQ=255)
                     tGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGGTTTTACATCATTGAGtttttt                                                                                                                                         >  1:1859703/1‑146 (MQ=255)
                     tGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGGAGATTTTACATCATCGGGttctt                                                                                                                                          >  1:1781295/1‑146 (MQ=255)
                       tctcTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGGATTTACATCATCGAggtttt                                                                                                                                          >  2:1612297/1‑139 (MQ=255)
                        ctctTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTTAccacattgtgttatt                                                                                                                                          >  1:1659658/1‑129 (MQ=255)
                        ctctTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTTACATCATCGAGttttttatc                                                                                                                                      >  2:580864/1‑143 (MQ=255)
                        ctctTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGGAGATTTTTCATCATCGAGtttttttgca                                                                                                                                     >  2:1847120/1‑143 (MQ=255)
                            tCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGaaa                                                                                                                                 <  2:31588/148‑1 (MQ=255)
                            tCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGaaa                                                                                                                                 <  1:1026725/148‑1 (MQ=255)
                             ccGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAAc                                                                                                                              <  2:1558312/150‑1 (MQ=255)
                               ggCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTTATATCATCGGGTTTTTTAGCAGaaaaaac                                                                                                                             >  2:1004538/1‑147 (MQ=255)
                               ggCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTTACTACATCGCGTTTTTTAGCAGaaaaaa                                                                                                                              >  1:1346595/1‑147 (MQ=255)
                               ggCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGGAGATTTGACATTATCCAGGTTTTTAGCAGaaaaaac                                                                                                                             >  1:1280132/1‑147 (MQ=255)
                                   aCAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTg                                                                                                                       <  1:507246/151‑1 (MQ=255)
                                    cAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGGAGATTTTACATCATCGAGTTCTTCAGCCGAAAAAACAAGCt                                                                                                                        >  1:1118329/1‑149 (MQ=255)
                                       gAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCa                                                                                                                   <  1:1041420/151‑1 (MQ=255)
                                        aaCGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCGCGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCa                                                                                                                   <  1:1545586/150‑1 (MQ=255)
                                        aaCGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTTACATCAACGAGGTCTTTATCAAAAAAAACaagatg                                                                                                                       >  1:1982672/1‑141 (MQ=255)
                                          cGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATg                                                                                                               <  2:1946728/151‑1 (MQ=255)
                                           ggCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATgg                                                                                                              <  1:1403396/151‑1 (MQ=255)
                                            gCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGt                                                                                                             <  1:80252/151‑1 (MQ=255)
                                             cGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGGAGATTTTACATCATCGAGTTCTTTAGCAGAAAAAACAAGCTGATCAAAgggt                                                                                                            >  1:309161/1‑149 (MQ=255)
                                                  tattTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCa                                                                                                       <  1:1800065/150‑1 (MQ=255)
                                                  tattTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCa                                                                                                       <  1:1980175/150‑1 (MQ=255)
                                                  tattTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGAATTTACATTATCGAGTTTTTTTGGAGAAAAAACAAGCTGTTCAAAGGGTgtc                                                                                                         >  2:674596/1‑147 (MQ=255)
                                                   attTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTTGAATCTTAGAGGTCTTTTAAAGAAAAAACATGCCGATCAAAGGGTGccca                                                                                                      >  1:1017157/1‑148 (MQ=255)
                                                   attTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCaa                                                                                                      <  1:2092627/150‑1 (MQ=255)
                                                   attTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGAATTTACATCATCCAGGTCTTTAGCAGAAAAAACATGCTGATCAAAAGTTGGCaa                                                                                                      >  2:1727338/1‑150 (MQ=255)
                                                     ttAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAAct                                                                                                    <  2:2059882/150‑1 (MQ=255)
                                                     ttAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAAct                                                                                                    <  1:685900/150‑1 (MQ=255)
                                                      tAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTTACATCATCGGGTTTTTTTGCAGAAAAAACAAGCCGATCAATAGTTGCCCAcc                                                                                                    >  2:1299079/1‑148 (MQ=255)
                                                         tGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGGAGATTTTTCATCTTCCCGGTTTTCTGTAGAAAAAACACTCTGATCAATAGTTggctaaccttcc                                                                                                >  1:1368355/1‑140 (MQ=255)
                                                           aGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCAt                                                                                             <  1:1260693/151‑1 (MQ=255)
                                                             cTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGGAGATTTTTCATCATCCAGTTCTTTAGCCGAAAAAACAAGCTGGTCAAAGGGTGCCAACTCTTCCCtt                                                                                            >  1:605444/1‑150 (MQ=255)
                                                               ggATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCTGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCATTa                                                                                           <  2:1723033/149‑1 (MQ=255)
                                                               ggATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGGAGATTTTAAATCCTTGGGGTTTTCAGCAAAAAAAACAAGCTGATCAAAGGTTGCCAACTCTGCCATTAAt                                                                                         >  2:1347942/1‑151 (MQ=255)
                                                               ggATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGGAGAATTTTAATCAACGAGGTCTTCAGCAGAAAAAACAAGCTGATCAATAGTTGCCAAATCTTCCTTTa                                                                                           >  1:1112900/1‑149 (MQ=255)
                                                                 aTGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCATTAATGc                                                                                       <  1:565705/151‑1 (MQ=255)
                                                                   gCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCATTAATGCa                                                                                      <  1:537995/150‑1 (MQ=255)
                                                                       ggCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCATTAATGCATTCgc                                                                                 <  1:1223566/151‑1 (MQ=255)
                                                                        gCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATTCTCGAGTTCTTCTGCAGGAAAAACAAGCTGATTAATGGTTGCTAACTCTTCCCTTAATGCAtt                                                                                    >  2:1638146/1‑147 (MQ=255)
                                                                        gCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGGAGATTTTTCATCATCCAGGTCTTCCTCCGAAAAACCAAGCTTATTAATGGGTGCCCAatctttcctt                                                                                            >  1:1651073/1‑133 (MQ=255)
                                                                        gCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGGAGAATTTACAACATCGAGTTCTTTAGCAGAAAAAACAAGCTGATCAATGGTTGCCCAACCCTCCATTAATGCATTCgc                                                                                 >  2:1436964/1‑150 (MQ=255)
                                                                         cGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCATTAATGCATTCgcag                                                                               <  1:1815766/151‑1 (MQ=255)
                                                                          gACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTTTCATCAACGGGTTCTTTAGCAGAAAAAACATGCTGTTCAATTGGTGCCCACTCTTCCATTTAAGCATTCgccg                                                                               >  2:1977294/1‑148 (MQ=255)
                                                                          gACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGGAGATTTTTCATCAACCAGTTCTTTACCAGAAAAAACAAGTTGATCAATTGGTTCCCAATCTCCCATTAATTCaattac                                                                                 >  1:973083/1‑143 (MQ=255)
                                                                              ggTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCATTAATGCATTCGCAGCAt                                                                            <  1:1175813/149‑1 (MQ=255)
                                                                                 cAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTTTCACTAACGAGTTCTTCTGCCGAAAAAACATGGCTATTAATTGTTGTCCACACTGCCCATAAATCATTTGCAGCATCCTGc                                                                       >  2:1532860/1‑151 (MQ=255)
                                                                                 cAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTTACATCATCGAGGTTTTTAGCAGAAAAAACATGCTGATCAATGGGTGCCAACTCTTCCATTTATTCATTTGCAGCATCCTGc                                                                       >  2:269515/1‑151 (MQ=255)
                                                                                  aGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCa                                                                      <  1:1540816/151‑1 (MQ=255)
                                                                                   gCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTTACATCATCGCGTTTTTTAGCAGAAAAAACATGCTGATCAATGGTTGCCAACTCTTCCATTAATTCATTCGCAGCATCCTGCa                                                                      >  2:1675831/1‑150 (MQ=255)
                                                                                    cATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCata                                                                    <  1:558536/151‑1 (MQ=255)
                                                                                    cATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCata                                                                    <  2:603565/151‑1 (MQ=255)
                                                                                    cATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGGAGATTTTTCAACAATGAGGTTTTTCGGAAAAAAAACAAACTGATTAATTGGTGCCAAACCTGCCCTTAAATCAATCccagcaaacct                                                                        >  2:600988/1‑142 (MQ=255)
                                                                                      tACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGGAGATTTTTCATCAACCAGGTTTTTAACCGAAAAAACAAGCTGATCAATGGTTTCCAACTCTCCCCTTTATGCATTCTCAGCATccctcctatt                                                                  >  1:980802/1‑143 (MQ=255)
                                                                                        cGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCAtattat                                                                <  2:1912592/151‑1 (MQ=255)
                                                                                           aaaTCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCg                                                              <  2:1799904/150‑1 (MQ=255)
                                                                                             aTCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGtct                                                           <  1:2051994/151‑1 (MQ=255)
                                                                                                   ttCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTTACATCATCGAGGTTTTTAGCAGAAAAAACATGCTGATCAAAGGGTGCCAACTCTGCCATTAATGCATTTGCAGCATCCTGCATATTATCGtctc                                                          >  2:210706/1‑146 (MQ=255)
                                                                                                       gCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGGAGATTTTTTATCAACGAGTTCTTTAGCAGAAAAAACATGCTGGTCCATGGTTTGCCACTCTTCCCTTTATGGATTTTCAGCCTCCTGCATATTATTTTTTCGTGCActt                                                  >  1:714152/1‑150 (MQ=255)
                                                                                                           gCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTacac                                              <  2:2043686/150‑1 (MQ=255)
                                                                                                           gCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTacac                                              <  2:1466219/150‑1 (MQ=255)
                                                                                                           gCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGAAAGATTTTACATCATCGGGATTTTCAGCCGCAAAACCACGCTGATTAACGGCTGCCCACCCTTCCAATAATGAATTACCAGCATCCTGCATATTATCATGTCGTGCATTTtca                                               >  2:241307/1‑149 (MQ=255)
                                                                                                             tGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGGAGATTTTACATCAACGAGGTCTTCAGCAGAAAAAACAAGCTGATCAAAGGGTTGCAACTCCTCCATTAAAGCATTTGCCGCATCCTGCCTATTATCGTCTCGGGCATTTacac                                              >  2:1474563/1‑148 (MQ=255)
                                                                                                              gATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACgg                                            <  1:1788226/149‑1 (MQ=255)
                                                                                                              gATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGa                                           <  2:1604695/150‑1 (MQ=255)
                                                                                                              gATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGa                                           <  1:1169302/150‑1 (MQ=255)
                                                                                                              gATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGa                                           <  2:1428008/150‑1 (MQ=255)
                                                                                                                 cAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGaaaa                                        <  1:2032322/150‑1 (MQ=255)
                                                                                                                  aaCCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAAc                                      <  2:692469/151‑1 (MQ=255)
                                                                                                                   aCCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACa                                     <  2:294266/151‑1 (MQ=255)
                                                                                                                         tttttCAGCCAGGTGAGGCGGGGAGAATTTACATCATCGAGTTTTTTTGCAGAAAAAACATGCTGGTCAAAGGGTGCCCACTCCGGCCTTAATGCAATTGCAGCATCCCGCATATTTTCGTCTCGTGCATTTATAACGaaaaaaa                                     >  1:675759/1‑143 (MQ=255)
                                                                                                                           tttCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTg                              <  1:2078442/150‑1 (MQ=255)
                                                                                                                            ttCAGCCAGGTGAGGCGGGTAGAATTTTCATCCTCGGGGTCTTTTTCAGAAAAAACAAGCCGGTCAATGGTTGGCAACTCTGCCCTTAATTTATTCTGCGCATTCTTCCTAATATTTTCTCGTGCATTTTTACGGAAAAAAAACGACTTAt                            >  1:1783292/1‑151 (MQ=255)
                                                                                                                             tCAGCCAGGTGAGGCGGGTAGATTTTAAATCATCGAGTTCTTCAGCAGAAAAAACATGCTGATCAAAGGGTGCCCACCCTTCCCTTAATGCATTCCCCGCATCCTGCATATTATTGTCTCGTGCATTTACACCGAAAAAAAACGACTGAt                            >  2:1912907/1‑150 (MQ=255)
                                                                                                                             tCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGAt                            <  2:218892/150‑1 (MQ=255)
                                                                                                                              cAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATa                           <  1:209489/150‑1 (MQ=255)
                                                                                                                              cAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATa                           <  2:625148/150‑1 (MQ=255)
                                                                                                                                gCCAGGTGAGGCGGGTAGATTTTACAACATCGAGTTCTTTAGCAGAAAAAACATGCTGATCCAAGGTTGCCCACTCTTCCATTAATTCATTTGCAGCATCCCACCTATTATCGTCCCCTGCATTTTCACGGAAAAAAAACGACTGATaa                          >  2:1445654/1‑149 (MQ=255)
                                                                                                                                 ccAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCg                       <  2:230144/151‑1 (MQ=255)
                                                                                                                                     gTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGTg                   <  2:1398418/151‑1 (MQ=255)
                                                                                                                                        aGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGTGCa                 >  1:1384253/1‑150 (MQ=255)
                                                                                                                                         ggCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGTGCAAt               >  1:1358702/1‑151 (MQ=255)
                                                                                                                                            gggTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGTGCAAtttt            <  1:390232/151‑1 (MQ=255)
                                                                                                                                                 gATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGTGCAATTTTAACGc       >  1:1458741/1‑151 (MQ=255)
                                                                                                                                                   tttGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCTACTGATAATCGGGTGCAATTTTAACGCt      <  1:1647784/150‑1 (MQ=255)
                                                                                                                                                    ttGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGTGCAATTTTAACGCTgc    >  2:853316/1‑151 (MQ=255)
                                                                                                                                                     tGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGTGCAATTTTAACGCTgc    <  1:2085311/150‑1 (MQ=255)
                                                                                                                                                      gACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTTCATTTACACGGAAAAACAGCGACTGATAATCGGGTGCAATTTTAACGCTgcg   <  1:585893/150‑1 (MQ=255)
                                                                                                                                                      gACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGTGCAATTTTAACGCTgcgc  <  1:1847027/151‑1 (MQ=255)
                                                                                                                                                      gACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGTGCAATTTTAACGCTgcg   <  1:749257/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                        
CAGACGTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGTGCAATTTTAACGCTGCGC  >  NC_000913/1573901‑1574201

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