Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,630,956 A→G 6.1% D75G (GAT→GGT)  slt → soluble lytic murein transglycosylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,630,9560AG 6.1% 247.4 / 3.3 81D75G (GAT→GGT) sltsoluble lytic murein transglycosylase
Reads supporting (aligned to +/- strand):  ref base A (41/35);  new base G (2/3);  total (43/38)
Fisher's exact test for biased strand distribution p-value = 6.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGACTCACTGGATGAGCAGCGTAGTCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAG  >  NC_000913/4630812‑4631097
                                                                                                                                                |                                                                                                                                             
ccGACTCACTGGATGAGCAGCGTAGTCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATc                                                                                                                                          >  1:973890/1‑150 (MQ=255)
  gACTCACTGGATGAGCAGCGTAGTCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCtg                                                                                                                                        <  1:226334/150‑1 (MQ=255)
     tCACTGTATGAGCATCGTACTCGTTACGCTGAAATCAAGCATGACTGGGATAATCGACAAATGGATGTGGTCTAACAAATGATGCCTGTACTGAAGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCATCGGTGATCtggcg                                                                                                                                     <  1:700193/150‑4 (MQ=255)
      cACTGGGTGAGCACCGTAGTCCTTACGCGCAAATCAAGCAGGCCTGTGATAATCGAAAAATGGATGTGGTAGAAAAAATGATGCCTGGACTGAAGGATTATCCTCTTTATCCCAACCTCGAATACCTCCAGATAACAGATCATCtcatga                                                                                                                                    >  2:1252532/1‑150 (MQ=255)
      cACTGGATGAGCAGCGTAGTCGTTACGCGCAAATGAAGCAGGCCTGGGATAATCGACAAATGGATGTGGGCGAACAAATGATGCCTGGACTGAAGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCATCGGTGATCtggcgg                                                                                                                                    <  1:701729/150‑5 (MQ=255)
      cACTGGATGAGCAGCGTAGTCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCATCGGTGATCtggcgg                                                                                                                                    <  1:1163524/150‑5 (MQ=255)
      cACTGGATGAGCAGCGTAGTCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCtgatga                                                                                                                                    >  1:840877/1‑150 (MQ=255)
        cTGGATGAGCAGCGTAGTCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAAt                                                                                                                                  >  1:746143/1‑150 (MQ=255)
         tGGATGTGCAGCGTAGTCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATc                                                                                                                                 <  2:1318914/150‑1 (MQ=255)
          ggATGAGCAGCGTAGTCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCa                                                                                                                                >  2:1014558/1‑150 (MQ=255)
               agcagcGTAGTCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTTTATTGAAGTATTATCCGGTTTATCCCTTCCTTGAATACCGCCTTATAATCGATGATCTGATGAATCTActtt                                                                                                                           <  2:1562226/150‑4 (MQ=255)
               agcagcGTAGTCGTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACcggc                                                                                                                          <  2:1547430/150‑1 (MQ=255)
                gcagcGTAGTCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACcggc                                                                                                                          >  1:24215/1‑150 (MQ=255)
                   gcGTAGTCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGt                                                                                                                       >  2:1202985/1‑150 (MQ=255)
                    cGTAGTCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTg                                                                                                                      <  2:1610868/150‑1 (MQ=255)
                      tAGTCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGAc                                                                                                                    >  2:227427/1‑150 (MQ=255)
                      tAGTCGTAACGCGCAAATCAAGCAGGCCTGTGATAATCGACAAATGGATGTGGTCCAACAAATGGTGCCTGGACTGAAGGATAATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGAAGATCTGATGAATCAACCGGCGGTGAc                                                                                                                    <  1:1565927/150‑1 (MQ=255)
                        gTCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACgg                                                                                                                  >  1:1096219/1‑150 (MQ=255)
                              aCGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTNACGGTCACTa                                                                                                            <  1:1014558/150‑1 (MQ=255)
                              aCGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTa                                                                                                            <  2:336752/150‑1 (MQ=255)
                               cgcgCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTaa                                                                                                           >  1:166028/1‑150 (MQ=255)
                                      aTCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGtt                                                                                                    >  1:831218/1‑150 (MQ=255)
                                        cAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTcg                                                                                                  <  2:380870/150‑1 (MQ=255)
                                         aaGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTcgc                                                                                                 <  2:973890/150‑1 (MQ=255)
                                           ggaGGCGTGGGATAATCGACAAAGGGATGTGGTCGAAGAAATGTTGCCTGGACTGAAGGATTATCCGCTTTATCGCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAAGCGGCGGTGACGGTCACTAACTTTGTTCgcgc                                                                                               <  2:496040/148‑1 (MQ=255)
                                              ggCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTaa                                                                                            >  2:814684/1‑150 (MQ=255)
                                                    ggATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAGGATTATCCACTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCAc                                                                                      >  1:1552996/1‑150 (MQ=255)
                                                         aTCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTtc                                                                                 >  1:773663/1‑150 (MQ=255)
                                                                    gATGTGGTCGAACAAATGATGCCTGGACTGAAGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGc                                                                      >  2:805399/1‑150 (MQ=255)
                                                                        tgGTCGAACAAATGATGCCTGGACTGAAGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGc                                                                  >  1:1445803/1‑150 (MQ=255)
                                                                          gTCGAACAAATGATGCCTGGACTGAAGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTg                                                                >  1:1210262/1‑150 (MQ=255)
                                                                             gAACAAATGATGCCTGGACTGAAGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCaa                                                             >  2:179250/1‑150 (MQ=255)
                                                                                cAAATGATGCCTGGACTGAAGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACAGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTAGCACGCTGAAAtct                                                          >  2:616966/1‑150 (MQ=255)
                                                                                 aaaTGATGCCTGGACTGAAGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAAtctc                                                         <  1:1438010/150‑1 (MQ=255)
                                                                                       tgCCTGGACTGAAGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCg                                                   >  2:1578627/1‑150 (MQ=255)
                                                                                        gCCTGGACTGAAGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGt                                                  <  1:218133/150‑1 (MQ=255)
                                                                                         ccTGGACTGAAGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTc                                                 <  2:24215/150‑1 (MQ=255)
                                                                                         ccTGGACTGAAGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTc                                                 <  2:1516550/150‑1 (MQ=255)
                                                                                         ccTGGACTGAAGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTc                                                 <  2:603402/150‑1 (MQ=255)
                                                                                               cTGAAGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGaa                                           >  1:685817/1‑150 (MQ=255)
                                                                                               cTGAAGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGaa                                           >  1:679990/1‑150 (MQ=255)
                                                                                                tGAAGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAAc                                          <  2:1096219/150‑1 (MQ=255)
                                                                                                 gAAGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACt                                         >  1:648354/1‑150 (MQ=255)
                                                                                                  aaGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCCTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTg                                        >  1:1299059/1‑150 (MQ=255)
                                                                                                   aGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTgg                                       <  1:875277/150‑1 (MQ=255)
                                                                                                         aTCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCgcggc                                 >  1:1601329/1‑150 (MQ=255)
                                                                                                         aTCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCgcggc                                 >  2:659311/1‑150 (MQ=255)
                                                                                                         aTCCGCTTTATCACTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATAAACCGGCGGTGACGGTCACTAACTTTGTTAGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCCTTTCGTCAATGAACTGGCgcggc                                 >  1:463181/1‑150 (MQ=255)
                                                                                                            cGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTg                              <  1:1109650/150‑1 (MQ=255)
                                                                                                             gCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGa                             <  1:1252532/150‑1 (MQ=255)
                                                                                                              cTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGaa                            >  2:730275/1‑150 (MQ=255)
                                                                                                               tttATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAg                           <  2:831218/150‑1 (MQ=255)
                                                                                                                ttATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGa                          <  1:1473178/150‑1 (MQ=255)
                                                                                                                 tATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTAACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGAc                         >  1:888938/1‑150 (MQ=255)
                                                                                                                  aTCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGt                                                                                                                       <  1:1135887/55‑1 (MQ=255)
                                                                                                                  aTCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGt                                                                                                                       >  2:1135887/1‑55 (MQ=255)
                                                                                                                  aTCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACt                        >  2:1048769/1‑150 (MQ=255)
                                                                                                                    cccTACCTGGAATACGGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTgg                      <  2:1069214/150‑1 (MQ=255)
                                                                                                                    cccTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTgg                      <  1:814684/150‑1 (MQ=255)
                                                                                                                     ccTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGc                     >  1:1089630/1‑150 (MQ=255)
                                                                                                                     ccTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGc                     >  1:162864/1‑150 (MQ=255)
                                                                                                                      cTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCg                    >  1:703043/1‑150 (MQ=255)
                                                                                                                      cTACCTGGAATACCGCCAGATCACCGACGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCg                    >  2:671016/1‑150 (MQ=255)
                                                                                                                        aCCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTg                  <  2:840877/150‑1 (MQ=255)
                                                                                                                         ccTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTgg                 <  2:162864/150‑1 (MQ=255)
                                                                                                                          cTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGc                >  1:343658/1‑150 (MQ=255)
                                                                                                                          cTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGc                <  2:773663/150‑1 (MQ=255)
                                                                                                                              aaTACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGt            <  1:1202985/150‑1 (MQ=255)
                                                                                                                              aaTACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGt            >  1:1521937/1‑150 (MQ=255)
                                                                                                                              aaTACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTAACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGt            <  2:888938/150‑1 (MQ=255)
                                                                                                                               aTACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGtt           <  1:1266449/150‑1 (MQ=255)
                                                                                                                                  tcatcAGATCATCGGTGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGc        >  1:17967/5‑150 (MQ=255)
                                                                                                                                  ccGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGc        <  2:1129881/150‑1 (MQ=255)
                                                                                                                                   catcAGATCATCGGTGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGcc       >  2:965797/4‑150 (MQ=255)
                                                                                                                                   cGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGcc       >  2:802849/1‑150 (MQ=255)
                                                                                                                                   cGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGcc       <  1:659311/150‑1 (MQ=255)
                                                                                                                                    gCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCt      >  1:15303/1‑150 (MQ=255)
                                                                                                                                    gCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCt      >  1:300298/1‑150 (MQ=255)
                                                                                                                                    gCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCt      <  2:1601329/150‑1 (MQ=255)
                                                                                                                                       aGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCCCGCTTCCTCCCGCTCGCCCGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTa   >  1:1514753/1‑150 (MQ=255)
                                                                                                                                        gATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAg  <  1:805399/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                             
CCGACTCACTGGATGAGCAGCGTAGTCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAG  >  NC_000913/4630812‑4631097

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