Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,781,020 +G intergenic (‑144/+65) ypjA ← / ← b2648 putative ATP‑binding component of a transport system/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,781,0191.G100.0% 705.2 / NA 213intergenic (‑143/+66)ypjA/b2648putative ATP‑binding component of a transport system/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (101/112);  total (101/112)

Only 100 of 239 total aligned reads displayed.

GGCATTTTTTAAACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGA‑GCCGTTGCCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTGAGATCG  >  NC_000913/2780927‑2781106
                                                                                             |                                                                                       
ggCATTTTTTAAACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGa                                                                                          <  2:6050189/93‑1 (MQ=255)
ggCATTTTTTAAACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGa                                                                                          <  2:2389392/93‑1 (MQ=255)
ggCATTTTTTAAACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGa                                                                                          <  2:2315625/93‑1 (MQ=255)
ggCATTTTTTAAACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGa                                                                                          <  2:1024022/93‑1 (MQ=255)
ggCATTTTTTAAACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGa                                                                                          <  2:4412352/93‑1 (MQ=255)
  cATTTTTTAAACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGA‑gg                                                                                       >  1:2416792/1‑92 (MQ=255)
   aTTTTTTAAACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGA‑ggc                                                                                      <  1:5474443/93‑3 (MQ=255)
   aTTTTTTAAACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGA‑ggc                                                                                      <  1:1945630/93‑3 (MQ=255)
    ttttttAAACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGA‑ggcc                                                                                     >  2:5451094/1‑90 (MQ=255)
       tttAAACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGtt                                                                                   <  1:5591205/93‑1 (MQ=255)
        ttAAACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTg                                                                                  <  1:4952672/93‑1 (MQ=255)
          aaaCGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTggc                                                                                >  2:6317498/1‑91 (MQ=255)
            aCGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCaa                                                                              >  2:2170292/1‑93 (MQ=255)
            aCGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCaa                                                                              <  2:4210844/93‑1 (MQ=255)
             cGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCaaa                                                                             <  1:2664702/93‑1 (MQ=255)
             cGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAAGCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCaaa                                                                             <  2:5313906/93‑1 (MQ=25)
              gCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCaaaa                                                                            <  2:3793752/93‑1 (MQ=25)
               cAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCaaaat                                                                           <  2:4536614/93‑1 (MQ=25)
               cAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCaaaat                                                                           >  1:3728566/1‑93 (MQ=255)
                aGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCaaaata                                                                          <  2:5584390/93‑1 (MQ=25)
                 gTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCaaaataa                                                                         >  1:1526387/1‑93 (MQ=255)
                     gcTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATg                                                                     >  2:710217/1‑93 (MQ=255)
                     gcTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATg                                                                     >  1:3084443/1‑93 (MQ=255)
                        gAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGAtt                                                                  >  1:6430584/1‑93 (MQ=255)
                         aaGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGAttt                                                                 >  1:3923721/1‑93 (MQ=255)
                          aGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTa                                                                <  2:5155167/93‑1 (MQ=37)
                           gtgtGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTAc                                                               >  2:1491195/1‑93 (MQ=255)
                           gtgtGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTAc                                                               >  1:851884/1‑93 (MQ=255)
                           gtgtGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTAc                                                               <  2:6145179/93‑1 (MQ=37)
                            tgtgGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACa                                                              <  2:3259086/93‑1 (MQ=37)
                            tgtgGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACa                                                              <  2:390070/93‑1 (MQ=37)
                            tgtgGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACa                                                              <  2:1246026/93‑1 (MQ=37)
                            tgtgGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACa                                                              <  1:6282385/93‑1 (MQ=255)
                             gtgGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACAt                                                             <  2:3982423/93‑1 (MQ=37)
                             gtgGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACAt                                                             <  1:3532217/93‑1 (MQ=37)
                             gtgGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATCATTTACAt                                                             <  1:4169814/93‑1 (MQ=37)
                               ggTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATaa                                                           <  2:978847/93‑1 (MQ=37)
                               ggTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATaa                                                           >  1:3720268/1‑93 (MQ=255)
                                gTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAAt                                                          <  1:1229657/93‑1 (MQ=37)
                                gTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAAt                                                          <  1:1426599/93‑1 (MQ=37)
                                gTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAAt                                                          <  2:3072196/93‑1 (MQ=37)
                                 ttGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATc                                                         >  1:6472415/1‑93 (MQ=255)
                                  tGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCg                                                        >  1:2409568/1‑93 (MQ=255)
                                   ggATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGt                                                       <  2:3473000/93‑1 (MQ=38)
                                    gATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGtt                                                      >  1:1046504/1‑93 (MQ=255)
                                    gATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGtt                                                      >  2:1169610/1‑93 (MQ=255)
                                    gATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGtt                                                      >  1:3123426/1‑93 (MQ=255)
                                    gATAAAAAAGTCAATCCATTCAGGAAATACGGGCGGATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGtt                                                      >  1:3880487/1‑93 (MQ=37)
                                     aTAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGttt                                                     >  1:5797157/1‑93 (MQ=255)
                                     aTAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGttt                                                     >  1:4553155/1‑93 (MQ=255)
                                       aaaaaaGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCt                                                   >  1:4083400/1‑93 (MQ=255)
                                          aaaGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCtgat                                                >  2:1648376/1‑93 (MQ=255)
                                           aaGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCtgatg                                               >  2:1051758/1‑93 (MQ=255)
                                            aGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCtgatga                                              >  2:4674974/1‑93 (MQ=255)
                                             gTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGaa                                             >  2:5835394/1‑93 (MQ=255)
                                               cAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAata                                           <  1:4288873/93‑1 (MQ=37)
                                                aaTCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAatat                                          >  1:6046958/1‑93 (MQ=255)
                                                aaTCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAatat                                          >  1:2774942/1‑93 (MQ=255)
                                                  tCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATAtct                                        <  2:5359944/93‑1 (MQ=255)
                                                  tCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATAtct                                        <  2:2929669/93‑1 (MQ=255)
                                                  tCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATAtct                                        <  2:5972230/93‑1 (MQ=25)
                                                    cATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATAtcttc                                      <  2:192848/93‑1 (MQ=255)
                                                     aTTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATAtcttct                                     >  2:3831148/1‑93 (MQ=255)
                                                      ttCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTg                                    >  1:3395247/1‑93 (MQ=255)
                                                        cAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCt                                  <  1:525758/93‑1 (MQ=21)
                                                         aGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTc                                 >  2:5016201/1‑93 (MQ=255)
                                                             aaTACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACAt                             <  2:4595998/93‑1 (MQ=21)
                                                              aTACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATa                            <  2:4516742/93‑1 (MQ=255)
                                                              aTACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATa                            >  2:4957829/1‑93 (MQ=18)
                                                              aTACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATa                            <  2:3758000/93‑1 (MQ=21)
                                                                  gggCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATaaaaa                        >  2:3960843/1‑93 (MQ=255)
                                                                   ggCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATACAAAt                       >  1:796965/1‑93 (MQ=17)
                                                                   ggCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAAt                       >  2:4321425/1‑93 (MQ=18)
                                                                      gTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATcac                    >  2:2067627/1‑93 (MQ=18)
                                                                       tATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATcaca                   >  2:4034058/1‑93 (MQ=255)
                                                                         ttcttttcttTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATcacaca                 >  1:4053186/1‑93 (MQ=17)
                                                                            ttttctttCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATa              <  2:5675575/93‑1 (MQ=255)
                                                                            ttttctttCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATa              <  1:2524364/93‑1 (MQ=255)
                                                                             tttctttcGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATaa             <  1:1472799/93‑1 (MQ=255)
                                                                              ttctttcGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAAc            <  1:5350268/93‑1 (MQ=255)
                                                                               tctttcGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACt           <  1:5720033/93‑1 (MQ=255)
                                                                               tctttcGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACt           >  1:4900475/1‑93 (MQ=17)
                                                                               tctttcGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACt           <  2:1352348/93‑1 (MQ=255)
                                                                               tctttcGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACt           <  1:451542/93‑1 (MQ=255)
                                                                               tctttcGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACt           <  1:3240296/93‑1 (MQ=14)
                                                                               tctttcGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACt           <  2:2621871/93‑1 (MQ=17)
                                                                               tctttcGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACt           <  2:4563329/93‑1 (MQ=17)
                                                                               tctttcGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACt           <  2:2074371/93‑1 (MQ=255)
                                                                                ctttcGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACtt          <  1:165150/93‑1 (MQ=17)
                                                                                ctttcGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACtt          <  2:3725456/93‑1 (MQ=17)
                                                                                 tttcGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACttt         <  1:4748144/93‑1 (MQ=255)
                                                                                 tttcGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACttt         >  1:639033/1‑93 (MQ=255)
                                                                                  ttcGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTg        >  1:4676114/1‑93 (MQ=17)
                                                                                  ttcGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTg        <  2:6228752/93‑1 (MQ=255)
                                                                                  ttcGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTg        <  1:1000824/93‑1 (MQ=255)
                                                                                    cGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTgag      <  2:886343/93‑1 (MQ=255)
                                                                                      aCAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTGAGAt    >  1:4389595/1‑93 (MQ=18)
                                                                                      aCAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTGAGAt    <  1:6145103/93‑1 (MQ=255)
                                                                                       cAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTGAGATc   >  2:2415675/1‑93 (MQ=255)
                                                                                       cAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTGAGATc   <  2:3619364/93‑1 (MQ=255)
                                                                                             |                                                                                       
GGCATTTTTTAAACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGA‑GCCGTTGCCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTGAGATCG  >  NC_000913/2780927‑2781106

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑