Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,564,525 G→A 100% M1M (ATG→ATA) † yffS → CPZ‑55 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,564,5250GA100.0% 312.5 / NA 92M1M (ATG→ATA) †yffSCPZ‑55 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (46/46);  total (46/46)

TGATCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGC  >  NC_000913/2564429‑2564620
                                                                                                |                                                                                               
tgATCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGt                                                                                              >  1:499736/1‑100 (MQ=255)
 gATCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGtt                                                                                             <  1:795048/100‑1 (MQ=255)
 gATCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGtt                                                                                             >  2:2510800/1‑100 (MQ=255)
   tCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTAc                                                                                           <  2:582209/100‑1 (MQ=255)
   tCTGCGACAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGTGATGGGAGATAACAATCATAAGTTa                                                                                            >  2:2342384/1‑99 (MQ=255)
      gCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGaa                                                                                        <  1:2155442/100‑1 (MQ=255)
       cGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGaaa                                                                                       >  1:2007549/1‑100 (MQ=255)
       cGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGaaa                                                                                       >  1:2007467/1‑100 (MQ=255)
        gTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAAt                                                                                      <  1:665250/100‑1 (MQ=255)
         tCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATc                                                                                     >  2:2485635/1‑100 (MQ=255)
          cAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCa                                                                                    >  1:728740/1‑100 (MQ=255)
              atTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAAt                                                                                <  2:1456552/100‑1 (MQ=255)
                 cccTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTg                                                                             <  1:2473777/100‑1 (MQ=255)
                  ccTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGc                                                                            >  1:143604/1‑100 (MQ=255)
                  ccTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGc                                                                            >  2:284154/1‑100 (MQ=255)
                   cTTCCGTTACGCTTACCCTCCACCAGTACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCg                                                                           <  1:2010937/100‑1 (MQ=255)
                   cTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCg                                                                           <  2:1830331/100‑1 (MQ=255)
                   cTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCg                                                                           <  2:342589/100‑1 (MQ=255)
                   cTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCg                                                                           <  1:875349/100‑1 (MQ=255)
                   cTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAAAAAAGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCg                                                                           <  2:119212/100‑1 (MQ=255)
                    ttCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGa                                                                          <  2:2007453/100‑1 (MQ=255)
                    ttCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGa                                                                          <  2:722850/100‑1 (MQ=255)
                    ttCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGa                                                                          <  2:2007535/100‑1 (MQ=255)
                     tCCGTTACGCTTACCCTCCACCAATCCGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAAGTTGCGa                                                                          >  2:487565/1‑99 (MQ=255)
                       cGTTACGCTTTCCCTCCACCATTTCGCGGGTTAATCTCTCAATGGATTCTATGAGATGTGATATAAAAATCATAATTTTCGAAATCAAAATTTGCGa                                                                          <  2:1115150/97‑1 (MQ=255)
                       cGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACAt                                                                       >  1:1646729/1‑100 (MQ=255)
                            cGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCa                                                                  >  1:206483/1‑100 (MQ=255)
                            cGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCa                                                                  >  2:1353258/1‑100 (MQ=255)
                              cTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCaaa                                                                <  1:862375/100‑1 (MQ=255)
                              cTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCaaa                                                                >  2:861322/1‑100 (MQ=255)
                                 aCCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGgc                                                             <  1:200400/100‑1 (MQ=255)
                                 aCCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCAATGAGATGGGAGATAAAAATCAAAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGgc                                                             >  1:1812725/1‑100 (MQ=255)
                                  cccTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGgcg                                                            >  1:765444/1‑100 (MQ=255)
                                  cccTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGgcg                                                            >  1:2248486/1‑100 (MQ=255)
                                   ccTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGgcgc                                                           >  1:50948/1‑100 (MQ=255)
                                    cTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGgtgc                                                           >  1:1659936/1‑99 (MQ=255)
                                    cTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCt                                                          >  1:1635747/1‑100 (MQ=255)
                                     tCCACCAATACGCGTGTTAACCTCCCAATGTATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTg                                                         <  2:1699555/100‑1 (MQ=255)
                                     tCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTg                                                         <  1:829499/100‑1 (MQ=255)
                                      ccaccaATACGCGTGTTAACCTCCCAACGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGc                                                        <  1:2432482/100‑1 (MQ=255)
                                       caccaATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAACAATCATAAGTTACGAAATCAAAATTTGCGAAAGCCTCAAAGGCGCTGCg                                                       >  2:1239873/1‑100 (MQ=255)
                                       caccaATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCg                                                       <  2:631668/100‑1 (MQ=255)
                                         ccaATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGAt                                                     >  1:1772920/1‑100 (MQ=255)
                                          caATGCGCGTGTTAACCTCCCAATGGAGTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATg                                                    <  2:1270895/100‑1 (MQ=255)
                                          caATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATg                                                    <  1:2268913/100‑1 (MQ=255)
                                          caATACGCGTGTTAACCTCCCAATCGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATg                                                    <  2:1402059/100‑1 (MQ=255)
                                            aTACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGa                                                  >  1:598770/1‑100 (MQ=255)
                                            aTACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGa                                                  <  2:871609/100‑1 (MQ=255)
                                             tACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGaa                                                 >  1:1910307/1‑100 (MQ=255)
                                              aggcgTGTTAACCTGGCAATGGATTCTATGAGATGGGAGATAAAAAGCATAAGTTACGAAAACAAAATTTGCGACATCCTCAAAGGCGCGGCGATgg                                                   <  2:2331742/95‑1 (MQ=255)
                                               cgcgTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAgg                                               >  1:1971582/1‑100 (MQ=255)
                                                 cgTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGAc                                             <  2:760157/100‑1 (MQ=255)
                                                 cgTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGAc                                             <  1:1407882/100‑1 (MQ=255)
                                                  gtgtTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGAc                                             <  2:223146/99‑1 (MQ=255)
                                                  gtgtTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACa                                            <  1:516671/100‑1 (MQ=255)
                                                   tgtTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACaa                                           >  2:1015144/1‑100 (MQ=255)
                                                   tgtTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACaa                                           >  1:2145134/1‑100 (MQ=255)
                                                   tgtTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATATGCGACATCCTCAAAGGCGCTGCGATGGAAGGACaa                                           >  1:1593148/1‑100 (MQ=255)
                                                     ttAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAata                                         >  1:1638505/1‑100 (MQ=255)
                                                      tAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAatat                                        >  1:1404902/1‑100 (MQ=255)
                                                        aCCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATaa                                      <  2:206481/100‑1 (MQ=255)
                                                         ccTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATaan                                     <  2:992050/100‑2 (MQ=255)
                                                         ccTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAg                                     <  2:1971568/100‑1 (MQ=255)
                                                           tCCCAATGGATTCTATGAGATGGGCGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGAAAGGACAATATAAggg                                   >  2:1495125/1‑100 (MQ=255)
                                                           tCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAggg                                   <  2:1623610/100‑1 (MQ=255)
                                                           tCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAggg                                   <  1:2155603/100‑1 (MQ=255)
                                                            cccAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGgc                                  <  1:1579895/100‑1 (MQ=255)
                                                            cccAATGGATTCGATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGgc                                  >  1:1681906/1‑100 (MQ=255)
                                                                  ggATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCaa                            <  1:2047268/100‑1 (MQ=255)
                                                                   gATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAAc                           >  1:1025350/1‑100 (MQ=255)
                                                                        tATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTggg                      <  1:17632/100‑1 (MQ=255)
                                                                         aTGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTgggg                     <  2:728733/100‑1 (MQ=255)
                                                                         aTGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTgggg                     <  1:1899307/100‑1 (MQ=255)
                                                                          tGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGc                    >  1:1015190/1‑100 (MQ=255)
                                                                          tGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGc                    >  1:1301830/1‑100 (MQ=255)
                                                                            agaTGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCa                  >  2:748123/1‑100 (MQ=255)
                                                                               tGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAAt               >  1:2258463/1‑100 (MQ=255)
                                                                                gggAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATg              >  2:689780/1‑100 (MQ=255)
                                                                                 ggAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGc             >  1:1920221/1‑100 (MQ=255)
                                                                                 ggAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGc             >  2:2474368/1‑100 (MQ=255)
                                                                                  gagaTAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCg            <  2:1868494/100‑1 (MQ=255)
                                                                                    gaTAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCgag          <  1:668542/100‑1 (MQ=255)
                                                                                     aTAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCgagg         >  1:1222164/1‑100 (MQ=255)
                                                                                     aTAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCgagg         <  1:10685/100‑1 (MQ=255)
                                                                                      tAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCgagga        >  2:236988/1‑100 (MQ=255)
                                                                                          aaTCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATc    <  2:257666/100‑1 (MQ=255)
                                                                                           aTCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGCTTGAAGGACAATATAAGGGCGCTCAAGGTCGGGCCAAATGCGCGgaga      >  2:2269075/1‑97 (MQ=255)
                                                                                           aTCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGATAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCg   >  1:244390/1‑100 (MQ=255)
                                                                                           aTCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCg   >  2:1707871/1‑100 (MQ=255)
                                                                                            tCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGc  <  2:1638493/100‑1 (MQ=255)
                                                                                            tCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGc  <  2:1466211/100‑1 (MQ=255)
                                                                                            tCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGc  >  1:955259/1‑100 (MQ=255)
                                                                                                |                                                                                               
TGATCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGC  >  NC_000913/2564429‑2564620

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