Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 477,435 T→C T395A (ACT→GCT)  ylaB ← putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913477,4350TC100.0% 310.9 / NA 85T395A (ACT→GCT) ylaBputative membrane‑anchored cyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (39/46);  total (39/46)

CGACATCCAGATCCTGTAAATAACTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGTGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTCCGAGGTGAGCACGGGGGATTCAAGATTGATCGAAATATGCTGCTGTGGATGCTGACGC  >  NC_000913/477300‑477568
                                                                                                                                       |                                                                                                                                     
cGACATCCAGATCCTGTAAATAACTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAg                                                                                                                                    <  1:383737/139‑1 (MQ=255)
     tCCAGATCCTGTAAATAACTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCaa                                                                                                                               >  1:65748/1‑139 (MQ=255)
          aTCCTGTAAATAACTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGcgcg                                                                                                                          <  1:1005333/139‑1 (MQ=255)
            ccTGTAAATAACTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAt                                                                                                                        <  2:849617/139‑1 (MQ=255)
            ccTGTAAATAACTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAt                                                                                                                        <  2:65748/139‑1 (MQ=255)
             cTGTAAATAACTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAtc                                                                                                                       <  2:229810/139‑1 (MQ=255)
               gTAAATAACTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAtctg                                                                                                                     >  2:777919/1‑139 (MQ=255)
                tAAATAACTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAtctgt                                                                                                                    >  2:1739712/1‑139 (MQ=255)
                tAAATAACTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAtctgt                                                                                                                    >  2:921272/1‑139 (MQ=255)
                tAAATAACTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAt                                                                                                                        >  2:275596/1‑135 (MQ=255)
                tAAATAACTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAt                                                                                                                        <  1:275596/135‑1 (MQ=255)
                 aaaTAACTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAtctgtc                                                                                                                   <  1:1171757/139‑1 (MQ=255)
                       cTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGtt                                                                                                                                  >  2:306391/1‑118 (MQ=255)
                       cTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGtt                                                                                                                                  <  1:306391/118‑1 (MQ=255)
                       cTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGAt                                                                                                             <  1:777919/139‑1 (MQ=255)
                          aaaCTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAt                                                                                                                        <  1:118964/125‑1 (MQ=255)
                          aaaCTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAt                                                                                                                        >  2:118964/1‑125 (MQ=255)
                           aaCTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAAc                                                                                                         <  2:1028596/139‑1 (MQ=255)
                                 aaTACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAgcg                                                                                                                            >  2:1204147/1‑114 (MQ=255)
                                 aaTACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAgcg                                                                                                                            <  1:1204147/114‑1 (MQ=255)
                                 aaTACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGcgcg                                                                                                                          >  1:344413/1‑116 (MQ=255)
                                 aaTACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGcgcg                                                                                                                          <  2:344413/116‑1 (MQ=255)
                                 aaTACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAt                                                                                                                        <  1:1213097/118‑1 (MQ=255)
                                 aaTACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAt                                                                                                                        >  2:1213097/1‑118 (MQ=255)
                                     ccccGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGtga                                                                                               <  1:921272/139‑1 (MQ=255)
                                        cGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAt                                                                                                                        <  1:1829800/111‑1 (MQ=255)
                                        cGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAt                                                                                                                        >  2:1829800/1‑111 (MQ=255)
                                        cGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGtga                                                                                               <  2:1061069/136‑1 (MQ=255)
                                        cGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGtga                                                                                               >  1:1061069/1‑136 (MQ=255)
                                           aCCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATc                                                                                         >  1:642165/1‑139 (MQ=255)
                                           aCCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATc                                                                                         >  2:410665/1‑139 (MQ=255)
                                                    atcaAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAt                                                                                                                        >  1:1344637/1‑99 (MQ=255)
                                                    atcaAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAt                                                                                                                        >  2:296906/1‑99 (MQ=255)
                                                    atcaAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAt                                                                                                                        <  2:1344637/99‑1 (MQ=255)
                                                    atcaAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAt                                                                                                                        <  1:296906/99‑1 (MQ=255)
                                                    atcaAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGAt                                                                                                    >  2:569632/1‑119 (MQ=255)
                                                    atcaAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGAt                                                                                                    <  1:569632/119‑1 (MQ=255)
                                                     tcaAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAt                                                                                                                        <  1:525012/98‑1 (MQ=255)
                                                     tcaAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAt                                                                                                                        >  2:525012/1‑98 (MQ=255)
                                                       aaGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTTAAGCGCGAt                                                                                                                        >  2:1582710/1‑96 (MQ=255)
                                                       aaGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAt                                                                                                                        <  1:1582710/96‑1 (MQ=255)
                                                              aTTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGtt                                                                                                                                  >  1:513372/1‑79 (MQ=255)
                                                              aTTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGtt                                                                                                                                  <  2:513372/79‑1 (MQ=255)
                                                               tttCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGAt                                                                                                             <  1:1235042/99‑1 (MQ=255)
                                                               tttCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGAt                                                                                                             >  2:1235042/1‑99 (MQ=255)
                                                                   aTGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAtctgtc                                                                                                                   >  2:639683/1‑89 (MQ=255)
                                                                   aTGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAtctgtc                                                                                                                   <  1:639683/89‑1 (MQ=255)
                                                                    tGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGtgatt                                                                                             >  1:1672793/1‑110 (MQ=255)
                                                                    tGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGtgatt                                                                                             <  2:1672793/110‑1 (MQ=255)
                                                                       cccGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAt                                                                                                                        <  2:701639/80‑1 (MQ=255)
                                                                       cccGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAt                                                                                                                        >  1:701639/1‑80 (MQ=255)
                                                                       cccGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGAtttttt                                                             >  2:518889/1‑139 (MQ=255)
                                                                          gCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTCCg                                                          >  2:1419047/1‑139 (MQ=255)
                                                                           ccTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGtt                                                                                                                                  <  1:1578623/66‑1 (MQ=255)
                                                                           ccTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGtt                                                                                                                                  >  2:1578623/1‑66 (MQ=255)
                                                                           ccTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGAt                                                                                                             >  1:1183045/1‑87 (MQ=255)
                                                                           ccTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGAt                                                                                                             <  2:1183045/87‑1 (MQ=255)
                                                                              cccGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTCCGAGGt                                                      <  1:518889/139‑1 (MQ=255)
                                                                              cccGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTCCGAGGt                                                      <  2:338437/139‑1 (MQ=255)
                                                                                 ggTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAt                                                                                                                        >  1:298667/1‑70 (MQ=255)
                                                                                 ggTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAt                                                                                                                        <  2:298667/70‑1 (MQ=255)
                                                                                   tAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTCCGAGGTGAGCa                                                 <  2:1462242/139‑1 (MQ=255)
                                                                                     gCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTCCGAGGTGAGCACg                                               >  2:1501512/1‑139 (MQ=255)
                                                                                     gCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTCCGAGGTGAGCACg                                               <  1:410665/139‑1 (MQ=255)
                                                                                      cGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAg                                                                                                                              <  2:1596879/59‑1 (MQ=255)
                                                                                      cGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAg                                                                                                                              >  1:1596879/1‑59 (MQ=255)
                                                                                      cGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAt                                                                                                                        >  2:47948/1‑65 (MQ=255)
                                                                                      cGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAt                                                                                                                        <  1:47948/65‑1 (MQ=255)
                                                                                       gagaAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTCCGAGGTGAGCACggg                                             <  1:1739712/139‑1 (MQ=255)
                                                                                        agaAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCaa                                                                                                                               >  2:1715832/1‑56 (MQ=255)
                                                                                        agaAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCaa                                                                                                                               <  1:1715832/56‑1 (MQ=255)
                                                                                        agaAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAt                                                                                                                        <  2:45916/63‑1 (MQ=255)
                                                                                        agaAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAt                                                                                                                        >  1:45916/1‑63 (MQ=255)
                                                                                           aaTTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAt                                                                                                                        >  1:652731/1‑60 (MQ=255)
                                                                                           aaTTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAt                                                                                                                        <  2:652731/60‑1 (MQ=255)
                                                                                                    ggCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTCCGAGGTGAGCACg                                               <  2:822239/124‑1 (MQ=255)
                                                                                                    ggCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTCCGAGGTGAGCACg                                               >  1:822239/1‑124 (MQ=255)
                                                                                                     gcgcTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTcc                                                           >  1:1487026/1‑111 (MQ=255)
                                                                                                     gcgcTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTcc                                                           <  2:1487026/111‑1 (MQ=255)
                                                                                                     gcgcTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTCCGAGGTGAGCACGGGGGATTCAAGATTGa                               <  1:1419047/139‑1 (MQ=255)
                                                                                                     gcgcTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTCCGAGGTGAGCACGGGGGATTCAAGATTGa                               <  2:833632/139‑1 (MQ=255)
                                                                                                                  gATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTCCGAGGTGAGCACGGGGGATTCAAGATTGATCGAAATAtgctg                  <  1:174656/139‑1 (MQ=255)
                                                                                                                      ggCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTCCGAGGt                                                      >  2:120747/1‑99 (MQ=255)
                                                                                                                      ggCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTCCGAGGt                                                      <  1:120747/99‑1 (MQ=255)
                                                                                                                                  ttCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTCCGAGGTGAGCACGGGGGATTCAAGATTGATCGAAATATGCTGCTGTGGATGCTGACGc  >  2:122473/1‑139 (MQ=255)
                                                                                                                                       |                                                                                                                                     
CGACATCCAGATCCTGTAAATAACTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGTGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTCCGAGGTGAGCACGGGGGATTCAAGATTGATCGAAATATGCTGCTGTGGATGCTGACGC  >  NC_000913/477300‑477568

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑