Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,450,334 G→A A64T (GCT→ACT)  ytfT → putative transport system permease protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,450,3340GA98.9% 634.6 / NA 188A64T (GCT→ACT) ytfTputative transport system permease protein
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base A (88/97);  total (88/100)
Fisher's exact test for biased strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 193 total aligned reads displayed.

TCGATAGCCTGGTAGCCCCGCATTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAGCTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTAGGGGCGGTGATGGCTATCGCC  >  NC_000913/4450245‑4450423
                                                                                         |                                                                                         
tCGATAGCCTGGTAGCCCCGCATTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTc                                                                                        <  1:2007762/93‑1 (MQ=255)
 cGATAGCCTGGTAGCCCCGCATTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTcc                                                                                       <  1:1013664/93‑1 (MQ=255)
  gATAGCCTGGTAGCCCCGCATTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTccc                                                                                      >  1:4202193/1‑93 (MQ=255)
   aTAGCCTGGTAGCCCCGCATTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCg                                                                                     >  2:3922827/1‑93 (MQ=255)
    tAGCCTGGTAGCCCCGCATTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGt                                                                                    >  1:148151/1‑93 (MQ=255)
     aGCCTGGTAGCCCCGCATTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGtt                                                                                   <  2:3574928/93‑1 (MQ=255)
     aGCCTGGTAGCCCCGCATTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGtt                                                                                   <  1:2871127/93‑1 (MQ=255)
     aGCCTGGTAGCCCCGCATTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGtt                                                                                   <  2:4640116/93‑1 (MQ=255)
     aGCCTGGTAGCCCCGCATTTCTGGCAGGTGGTGCTTCAGGATGGGCGCTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGtt                                                                                   >  2:3631381/1‑93 (MQ=255)
      gCCTGGTAGCGCCGCATTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTg                                                                                  <  1:3827821/93‑1 (MQ=255)
      gCCTGGTAGCCCCGCATTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTg                                                                                  >  1:2380518/1‑93 (MQ=255)
       ccTGGTAGCCCCGCATTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGc                                                                                 >  1:462368/1‑93 (MQ=255)
         tGGTAGCCCCGCATTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAGCCGTGCAACTCCCGTTGCGt                                                                               <  1:349386/93‑1 (MQ=255)
         tGGTAGCCCCGCATTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGt                                                                               <  1:230283/93‑1 (MQ=255)
            tAGCCCCGCATTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAAGTCCCGTTGCGTTat                                                                            >  1:1792015/1‑92 (MQ=255)
            tAGCCCCGCATTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTGCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTAc                                                                            <  2:2818475/93‑1 (MQ=255)
             aGCCCCGCATTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACt                                                                           <  2:168205/93‑1 (MQ=255)
              gCCGCGCATTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTg                                                                          <  1:4017783/93‑1 (MQ=255)
              gCCCCGCATTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTg                                                                          <  1:1759005/93‑1 (MQ=255)
              gCCCCGCATTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTg                                                                          >  2:3428343/1‑93 (MQ=255)
              gCCCCGCATTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTg                                                                          >  1:1960132/1‑93 (MQ=255)
              gCCCCGCATTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTg                                                                          <  2:4525944/93‑1 (MQ=255)
              gCCCCGCATTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTg                                                                          <  2:4565072/93‑1 (MQ=255)
               ccccGCATTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTgg                                                                         <  1:4497368/93‑1 (MQ=255)
                cccGCATTTCTGGCAGGGGGTGCTTCAGGCTGGGCGTTTGTTCGGTAGCCCCCTAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGc                                                                        >  2:4069144/1‑93 (MQ=255)
                  cGCATTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGa                                                                      >  2:1100617/1‑93 (MQ=255)
                  cGCATTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGa                                                                      >  1:1435748/1‑93 (MQ=255)
                   gCATTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGAt                                                                     <  1:2551965/93‑1 (MQ=255)
                   gCATTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGAt                                                                     >  2:4669545/1‑93 (MQ=255)
                     aTTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTACCCGTGCAACTCCCGTTGCGTTACTGGCGATTg                                                                   <  1:2515547/93‑1 (MQ=255)
                     aTTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTg                                                                   <  1:556563/93‑1 (MQ=255)
                     aTTTCTGGCAGGTGGGGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAGCCGTGCAGCTCCCGTTGCGTTACTGGCGATTg                                                                   <  2:559098/93‑1 (MQ=255)
                           ggCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGa                                                             <  2:1554186/93‑1 (MQ=255)
                            gCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGAc                                                            >  1:2740225/1‑93 (MQ=255)
                              aGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGc                                                          >  2:2586443/1‑93 (MQ=255)
                               ggtggtGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCt                                                         <  2:592939/93‑1 (MQ=255)
                                gtggtgCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTg                                                        >  1:3166389/1‑93 (MQ=255)
                                 tggtgCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACTTTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTgg                                                       <  1:2145086/93‑1 (MQ=255)
                                 tggtgCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTgg                                                       >  2:3660028/1‑93 (MQ=255)
                                  gttgCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGt                                                      <  2:86053/91‑1 (MQ=255)
                                  ggtgCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGt                                                      <  1:710313/93‑1 (MQ=255)
                                  ggtgCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGt                                                      <  2:1805168/93‑1 (MQ=255)
                                  ggtgCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGt                                                      >  2:224972/1‑93 (MQ=255)
                                  ggtgCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGAGACTGGCATGACGCGGGt                                                      >  1:297655/1‑93 (MQ=255)
                                   gtgCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTg                                                     <  2:3145648/93‑1 (MQ=255)
                                     gCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGAt                                                   <  2:817616/93‑1 (MQ=255)
                                          aGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCa                                              <  1:1175927/93‑1 (MQ=255)
                                          aGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCa                                              >  1:3660428/1‑93 (MQ=255)
                                           ggATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCAc                                             <  1:2347635/93‑1 (MQ=255)
                                           ggATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCAc                                             >  2:3761557/1‑93 (MQ=255)
                                            gATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTACCCGTGCAACTCCCGTTGCGTTCCTGGCGATTGGCATGACGCTGGTGATCGCCAcc                                            <  2:1450152/93‑1 (MQ=255)
                                             aTGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCg                                           >  2:1385913/1‑93 (MQ=255)
                                              tGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCgg                                          >  2:2442031/1‑93 (MQ=255)
                                              tGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCgg                                          >  1:3271636/1‑93 (MQ=255)
                                              gtgtcGTTTGTTCGGTAGCCCCATAGACATTCTTATCCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCgg                                          <  1:260259/89‑1 (MQ=255)
                                               gggCGTTTGTTCGGTAGCCCCATAGACATTCTTACCCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGt                                         <  2:2530179/93‑1 (MQ=255)
                                                ggCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTg                                        <  2:328198/93‑1 (MQ=255)
                                                ggCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTg                                        >  2:128121/1‑93 (MQ=255)
                                                 gCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTgg                                       <  2:1766331/93‑1 (MQ=255)
                                                 gCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTgg                                       <  1:585524/93‑1 (MQ=255)
                                                 gCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTgg                                       <  1:3190322/93‑1 (MQ=255)
                                                  cGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTggg                                      >  2:2333529/1‑93 (MQ=255)
                                                  cGTTTGTTCGGTAGCCCCATAGACATCCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTggg                                      <  2:2476877/93‑1 (MQ=255)
                                                   gTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGa                                     <  1:3916761/93‑1 (MQ=255)
                                                   gTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGa                                     <  1:2124765/93‑1 (MQ=255)
                                                   gTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGa                                     <  2:4211414/93‑1 (MQ=255)
                                                    tttGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGAt                                    >  1:305603/1‑93 (MQ=255)
                                                    tttGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGAt                                    >  1:3340324/1‑93 (MQ=255)
                                                    tttGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGAt                                    >  2:1879721/1‑93 (MQ=255)
                                                       gTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGa                                 <  1:1230919/93‑1 (MQ=255)
                                                       gTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGa                                 >  1:3617118/1‑93 (MQ=255)
                                                          cGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGAtct                              >  1:939954/1‑93 (MQ=255)
                                                            gTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGAtctct                            >  1:4536881/1‑93 (MQ=255)
                                                             tAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGAtctctc                           <  1:2234/93‑1 (MQ=255)
                                                             tAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGAtctctc                           <  1:1515775/93‑1 (MQ=255)
                                                              aGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTcc                          <  1:4365385/93‑1 (MQ=255)
                                                               gCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCg                         >  1:4643162/1‑93 (MQ=255)
                                                                 cccATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGTATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTa                       >  1:1723304/1‑93 (MQ=255)
                                                                 cccATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTa                       <  1:987654/93‑1 (MQ=255)
                                                                  ccATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTAg                      >  2:2251774/1‑93 (MQ=255)
                                                                  ccATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTAg                      >  2:1170140/1‑93 (MQ=255)
                                                                   cATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTAgg                     <  1:3642294/93‑1 (MQ=255)
                                                                   cATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTAgg                     >  1:166613/1‑93 (MQ=255)
                                                                    aTAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTAggg                    >  1:1269659/1‑93 (MQ=255)
                                                                    aTAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTAggg                    <  1:1860953/93‑1 (MQ=255)
                                                                     tAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTAgggg                   >  1:1030389/1‑93 (MQ=255)
                                                                     tAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTAgggg                   >  1:839835/1‑93 (MQ=255)
                                                                     tAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTAgggg                   >  2:174314/1‑93 (MQ=255)
                                                                       gACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTAGGGGCg                 >  2:2061575/1‑93 (MQ=255)
                                                                        aCATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGGGGGATTGATCTCTCCGTAGGGGCgg                >  2:3729091/1‑93 (MQ=255)
                                                                          aTTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTAGGGGCGGTg              >  2:1500789/1‑93 (MQ=255)
                                                                           ttCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTAGGGGCGGTGa             >  2:2162112/1‑93 (MQ=255)
                                                                            tCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTAGGGGCGGTGAt            <  2:3953850/93‑1 (MQ=255)
                                                                             cTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTAGGGGCGGTGATg           <  2:2725649/93‑1 (MQ=255)
                                                                               tAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTAGGGGCGGTGATGGc         <  1:1609263/93‑1 (MQ=255)
                                                                                aaCCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTAGGGGCGGTGATGGCt        >  1:627055/1‑93 (MQ=255)
                                                                                 aCCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTAGGGGCGGTGATGGCTa       >  1:484509/1‑93 (MQ=255)
                                                                                   cGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTAGGGGCGGTGATGGCTATc     >  1:3529851/1‑93 (MQ=255)
                                                                                     tGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTAGGGGCGGTGATGGCTATCGc   >  1:669649/1‑93 (MQ=255)
                                                                                      gCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTAGGGGCGGTGATGGCTATCGcc  >  1:2613796/1‑93 (MQ=255)
                                                                                         |                                                                                         
TCGATAGCCTGGTAGCCCCGCATTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAGCTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTAGGGGCGGTGATGGCTATCGCC  >  NC_000913/4450245‑4450423

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