Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,197,162 A→G V124A (GTT→GCT)  ymfD ← e14 prophage, putative SAM‑dependent methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,197,1620AG100.0% 323.9 / NA 89V124A (GTT→GCT) ymfDe14 prophage, putative SAM‑dependent methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (44/45);  total (44/45)

GTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGTCGACATCCTCGAAAGCAACTGTATTTGCATTTTTATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAAT  >  NC_000913/1197034‑1197289
                                                                                                                                |                                                                                                                               
gTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAGCGGCagagag                                                                                                                       <  2:1483866/139‑1 (MQ=255)
        ccTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAGCGGCagagag                                                                                                                       <  1:1157249/131‑1 (MQ=255)
        ccTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAGCGGCagagag                                                                                                                       >  2:1157249/1‑131 (MQ=255)
                    tttttGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAGCGGCAGAGAGAACATTAGAGCAAAGGATGa                                                                                                   <  2:1590910/139‑1 (MQ=255)
                               tGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAGCGGCAGAGAGAACATTAGAGCa                                                                                                           <  2:962043/120‑1 (MQ=255)
                               tGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAGCGGCAGAGAGAACATTAGAGCa                                                                                                           >  1:962043/1‑120 (MQ=255)
                                       tataTTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAGCGGCaga                                                                                                                          <  1:888789/97‑1 (MQ=255)
                                       tataTTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAGCGGCaga                                                                                                                          >  2:888789/1‑97 (MQ=255)
                                              aTTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAGCGGCagagag                                                                                                                       <  2:1519330/93‑1 (MQ=255)
                                              aTTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAGCGGCagagag                                                                                                                       >  1:1519330/1‑93 (MQ=255)
                                                ttACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAGCGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGTCGAc                                                                       <  1:229518/139‑1 (MQ=255)
                                                 tACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAGCGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGTCGACa                                                                      >  1:1059689/1‑139 (MQ=255)
                                                   cAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAGCGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGTCGACATc                                                                    >  2:1578915/1‑139 (MQ=255)
                                                         gagTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAGCGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGTCGACATCCTCGaa                                                              <  2:1059689/139‑1 (MQ=255)
                                                              tcACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAGCGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGTCGAc                                                                       >  2:469410/1‑125 (MQ=255)
                                                              tcACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAGCGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGTCGAc                                                                       <  1:469410/125‑1 (MQ=255)
                                                              tcACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAGCGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGTCGACATCCTCGAAAGCaa                                                         >  1:1132570/1‑139 (MQ=255)
                                                              tcACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAGCGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGTCGACATCCTCGAAAGCaa                                                         >  2:588289/1‑139 (MQ=255)
                                                              tcACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAGCGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGTCGACATCCTCGAAAGCaa                                                         <  1:588289/139‑1 (MQ=255)
                                                                aCCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAGCGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGTCGACATCCTCGAAAGCAACt                                                       >  1:432398/1‑139 (MQ=255)
                                                                  ctcctGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAGCGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGTCGACATCCTCGAAAGCaa                                                         >  2:1552196/1‑135 (MQ=255)
                                                                  ctcctGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAGCGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGTCGACATCCTCGAAAGCaa                                                         <  1:1552196/135‑1 (MQ=255)
                                                                    cctGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAGCGGCa                                                                                                                            >  2:58987/1‑66 (MQ=255)
                                                                    cctGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAGCGGCa                                                                                                                            <  1:58987/66‑1 (MQ=255)
                                                                    cctGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAGCGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGTCGACATCCTCGAAAGc                                                           >  1:718501/1‑131 (MQ=255)
                                                                    cctGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAGCGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGTCGACATCCTCGAAAGc                                                           <  2:718501/131‑1 (MQ=255)
                                                                     ctGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAGCGGCa                                                                                                                            >  1:269940/1‑65 (MQ=255)
                                                                     ctGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAGCGGCa                                                                                                                            <  2:269940/65‑1 (MQ=255)
                                                                      tGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAGCGGCagagag                                                                                                                       >  1:1374375/1‑69 (MQ=255)
                                                                      tGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAGCGGCagagag                                                                                                                       <  2:1374375/69‑1 (MQ=255)
                                                                      tGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAGCGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGTCGACATCCTCGAAAGCAACTGTAtt                                                  >  2:550390/1‑138 (MQ=255)
                                                                      tGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAGCGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGTCGACATCCTCGAAAGCAACTGTAtt                                                  <  1:550390/138‑1 (MQ=255)
                                                                          tttAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAGCGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGTCGACATCCTCGAAAGCAACTGTATTTGCAt                                             <  2:1285101/139‑1 (MQ=255)
                                                                            tAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAGCGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGTCGACATCCTCGAAAGCAACTGTATTTGCAttt                                           >  2:1576287/1‑139 (MQ=255)
                                                                                  tctcTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAGCGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGTCGACATCCTCGAAAGCAACTGTATTTGCATTTTTataa                                     >  1:1493776/1‑139 (MQ=255)
                                                                                  tctcTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAGCGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGTCGACATCCTCGAAAGCAACTGTATTTGCATTTTTataa                                     <  1:1576287/139‑1 (MQ=255)
                                                                                  tctcTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAGCGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGTCGACATCCTCGAAAGCAACTGTATTTGCATTTTTataa                                     <  2:432398/139‑1 (MQ=255)
                                                                                       tGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAGCGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGTCGACATCCTCGAAAGCAACTGTATTTGCATTTTTATAATATCg                                <  1:846495/139‑1 (MQ=255)
                                                                                         aTGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAGCGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGTCGACATCCTCGAAAGCAACTGTATTTGCATTTTTATAATATCGTg                              >  2:1575335/1‑139 (MQ=255)
                                                                                             tAAGAACTATTTTGTCGATTGTATCCCGACAAGGAGCGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCGCAATTATTTTGTCGACa                                                                      <  1:90821/95‑1 (MQ=255)
                                                                                             tAAGAACTATTTTGTCGATTGTATCCCGACAAGGAGCGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGTCGACa                                                                      >  2:90821/1‑95 (MQ=255)
                                                                                              aaGAACTATTTTGTCGATTGTATCCCGACAAGGAGCGGCAGAGAGAACATTAGAGCAAAGGATGaa                                                                                                  >  2:580366/1‑66 (MQ=255)
                                                                                              aaGAACTATTTTGTCGATTGTATCCCGACAAGGAGCGGCAGAGAGAACATTAGAGCAAAGGATGaa                                                                                                  <  1:580366/66‑1 (MQ=255)
                                                                                              aaGAACTATTTTGTCGATTGTATCCCGACAAGGAGCGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGTCGACa                                                                      <  2:333682/94‑1 (MQ=255)
                                                                                              aaGAACTATTTTGTCGATTGTATCCCGACAAGGAGCGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGTCGACa                                                                      <  2:1369430/94‑1 (MQ=255)
                                                                                              aaGAACTATTTTGTCGATTGTATCCCGACAAGGAGCGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGTCGACa                                                                      >  1:333682/1‑94 (MQ=255)
                                                                                              aaGAACTATTTTGTCGATTGTATCCCGACAAGGAGCGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGTCGACa                                                                      >  1:1369430/1‑94 (MQ=255)
                                                                                              aaGAACTATTTTGTCGATTGTATCCCGACAAGGAGCGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGTCGACATCCTc                                                                 >  1:624363/1‑99 (MQ=255)
                                                                                              aaGAACTATTTTGTCGATTGTATCCCGACAAGGAGCGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGTCGACATCCTc                                                                 <  2:624363/99‑1 (MQ=255)
                                                                                              aaGAACTATTTTGTCGATTGTATCCCGACAAGGAGCGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGTCGACATCCTCGAAAGCAACTGTATTTGCAt                                             >  2:71672/1‑119 (MQ=255)
                                                                                              aaGAACTATTTTGTCGATTGTATCCCGACAAGGAGCGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGTCGACATCCTCGAAAGCAACTGTATTTGCAt                                             <  1:71672/119‑1 (MQ=255)
                                                                                              aaGAACTATTTTGTCGATTGTATCCCGACAAGGAGCGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGTCGACATCCTCGAAAGCAACTGTATTTGCATTTTTataa                                     >  2:171823/1‑127 (MQ=255)
                                                                                              aaGAACTATTTTGTCGATTGTATCCCGACAAGGAGCGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGTCGACATCCTCGAAAGCAACTGTATTTGCATTTTTataa                                     <  1:171823/127‑1 (MQ=255)
                                                                                              aaGAACTATTTTGTCGATTGTATCCCGACAAGGAGCGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGTCGACATCCTCGAAAGCAACTGTATTTGCATTTTTATAATATCGTGGGACa                         >  1:1541425/1‑139 (MQ=255)
                                                                                               aGAACTATTTTGTCGATTGTATCCCGACAAGGAGCGGCAGAGAGAACATTag                                                                                                               >  2:936139/1‑52 (MQ=255)
                                                                                               aGAACTATTTTGTCGATTGTATCCCGACAAGGAGCGGCAGAGAGAACATTag                                                                                                               <  1:936139/52‑1 (MQ=255)
                                                                                                gAACTATTTTGTCGATTGTATCCCGACAAGGAGCGGCa                                                                                                                            >  1:904693/1‑38 (MQ=255)
                                                                                                gAACTATTTTGTCGATTGTATCCCGACAAGGAGCGGCa                                                                                                                            <  2:904693/38‑1 (MQ=255)
                                                                                                gAACTATTTTGTCGATTGTATCCCGACAAGGAGCGGCAGAGAGAACATTAGAGCa                                                                                                           >  2:256244/1‑55 (MQ=255)
                                                                                                gAACTATTTTGTCGATTGTATCCCGACAAGGAGCGGCAGAGAGAACATTAGAGCa                                                                                                           <  1:256244/55‑1 (MQ=255)
                                                                                                gAACTATTTTGTCGATTGTATCCCGACAAGGAGCGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAGTTATTTTGTCGACATCCTCGAAAGCaa                                                         >  1:1057764/1‑105 (MQ=255)
                                                                                                gAACTATTTTGTCGATTGTATCCCGACAAGGAGCGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGTCg                                                                         >  2:1339123/1‑89 (MQ=255)
                                                                                                gAACTATTTTGTCGATTGTATCCCGACAAGGAGCGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGTCg                                                                         <  1:1339123/89‑1 (MQ=255)
                                                                                                gAACTATTTTGTCGATTGTATCCCGACAAGGAGCGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGTCGACATCCTCGAAAGCAACTGTATTTGCATTTTTataa                                     <  1:44652/125‑1 (MQ=255)
                                                                                                gAACTATTTTGTCGATTGTATCCCGACAAGGAGCGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGTCGACATCCTCGAAAGCAACTGTATTTGCATTTTTataa                                     >  2:44652/1‑125 (MQ=255)
                                                                                                gAACTATTTTGTCGATCGTATCCCGACAAGGAGCGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGTCGACATCCTCGAAAGCaa                                                         <  2:1057764/105‑1 (MQ=255)
                                                                                                     aTTTTGTCGATTGTATCCCGACAAGGAGCGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGTCGACa                                                                      >  2:319743/1‑87 (MQ=255)
                                                                                                     aTTTTGTCGATTGTATCCCGACAAGGAGCGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGTCGACa                                                                      <  1:319743/87‑1 (MQ=255)
                                                                                                      ttttGTCGATTGTATCCCGACAAGGAGCGGCAGAGAGAACATTAGAGCa                                                                                                           <  2:699277/49‑1 (MQ=255)
                                                                                                      ttttGTCGATTGTATCCCGACAAGGAGCGGCAGAGAGAACATTAGAGCa                                                                                                           >  1:699277/1‑49 (MQ=255)
                                                                                                      ttttGTCGATTGTATCCCGACAAGGAGCGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCGCCAATTATTTTGTCGACATCCTCGa                                                               >  2:1042660/1‑93 (MQ=255)
                                                                                                      ttttGTCGATTGTATCCCGACAAGGAGCGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGTCGACATCCTCgc                                                               <  1:1042660/93‑2 (MQ=255)
                                                                                                              aTTGTATCCCGACAAGGAGCGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGTCGACATCCTCGAAAGCAACTGTATTTGCATTTTTATAATATCGTGGGACATAGTCAATAATTTTAg         >  2:484360/1‑139 (MQ=255)
                                                                                                               ttGTATCCCGACAAGGAGCGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGTCGACATCCTCGAAAGCAACTGTATTTGCATTTTTATAATATCGTGGGACATAGTCAATAATTTTAGt        <  2:1541425/139‑1 (MQ=255)
                                                                                                                    tCCCGATAAGGAGCGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGTCGACATCCTCGAAAGCAACTGTATTTGCATTTTTATAATATCGTGGGACATAGTCAATAATTTTAGTTTTaa   <  2:428823/139‑1 (MQ=255)
                                                                                                                    tCCCGACAAGGAGCGGCAGAGAGAACATTAGAGGAAAGGATGAAATCGTAACCCCCAATTATTTTGTCGACATCCTCGAAAGCAACTGTATTTGCATTTTTATAATATCGTGGGACATAGTCAATAATTTTAGTTTTaa   >  1:428823/1‑139 (MQ=255)
                                                                                                                     cccGACAAGGAGCGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGTCGACATCCTCGAAAGCAACTGTATTTGCATTTTTATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAAt  <  1:1578915/139‑1 (MQ=255)
                                                                                                                       cGACAAGGAGCGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGTCGACATCCTCGAAAGCaa                                                         >  2:432959/1‑82 (MQ=255)
                                                                                                                       cGACAAGGAGCGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGTCGACATCCTCGAAAGCaa                                                         <  1:432959/82‑1 (MQ=255)
                                                                                                                         aCAAGGAGCGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGTCGACATCCTCGAAAGCAAc                                                        <  2:1431219/81‑1 (MQ=255)
                                                                                                                         aCAAGGAGCGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGTCGACATCCTCGAAAGCAAc                                                        >  1:1431219/1‑81 (MQ=255)
                                                                                                                         aCAAGGAGCGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGTCGACATCCTCGAAAGCAACTGTATTTGc                                               >  2:681868/1‑90 (MQ=255)
                                                                                                                         aCAAGGAGCGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGTCGACATCCTCGAAAGCAACTGTATTTGc                                               <  1:681868/90‑1 (MQ=255)
                                                                                                                         aCAAGGAGCGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGTCGACATCCTCGAAAGCAACTGTATTTGCATTTTTataa                                     >  2:1399867/1‑100 (MQ=255)
                                                                                                                         aCAAGGAGCGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGTCGACATCCTCGAAAGCAACTGTATTTGCATTTTTataa                                     <  1:1399867/100‑1 (MQ=255)
                                                                                                                          cAAGGAGCGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAAcc                                                                                         <  2:988566/47‑1 (MQ=255)
                                                                                                                          cAAGGAGCGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAAcc                                                                                         >  1:988566/1‑47 (MQ=255)
                                                                                                                             ggAGCGGCAGAGAGAACATTAGAGCAAAGGCTGAAATCGTAACCCCCAATTATTTTGTCGACATCCTCGAAAGCaa                                                         >  1:891928/1‑76 (MQ=255)
                                                                                                                             ggAGCGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGTCGACATCCTCGAAAGCaa                                                         <  2:891928/76‑1 (MQ=255)
                                                                                                                                |                                                                                                                               
GTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGTCGACATCCTCGAAAGCAACTGTATTTGCATTTTTATAATATCGTGGGACATAGTCAATAATTTTAGTTTTAAT  >  NC_000913/1197034‑1197289

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