Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 932,454 A→G Q3R (CAG→CGG)  ftsK → cell division protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913932,4540AG100.0% 752.9 / NA 230Q3R (CAG→CGG) ftsKcell division protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (103/127);  total (103/127)

Only 100 of 239 total aligned reads displayed.

TACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCTGATCCTTATTGTCCTGT  >  NC_000913/932364‑932543
                                                                                          |                                                                                         
tACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCggg                                                                                         <  2:330222/93‑1 (MQ=255)
tACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCggg                                                                                         >  2:5760337/1‑93 (MQ=255)
tACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCggg                                                                                         >  2:4481262/1‑93 (MQ=255)
tACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGACAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCggg                                                                                         <  2:4259683/93‑1 (MQ=255)
   aGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAAt                                                                                      >  1:741179/1‑93 (MQ=255)
   aGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAAt                                                                                      >  2:3107514/1‑93 (MQ=255)
    gCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATa                                                                                     <  1:1083319/93‑1 (MQ=255)
     cAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATAc                                                                                    <  2:5079245/93‑1 (MQ=255)
     cAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATAc                                                                                    <  1:5606395/93‑1 (MQ=255)
       aCAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACAt                                                                                  >  1:5196599/1‑93 (MQ=255)
        cAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACAtt                                                                                 >  1:5931780/1‑93 (MQ=255)
         aGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTg                                                                                <  1:1421446/93‑1 (MQ=255)
         aGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTg                                                                                <  2:2972924/93‑1 (MQ=255)
          gTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGa                                                                               <  2:4327626/93‑1 (MQ=255)
           tACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGaa                                                                              <  2:5337643/93‑1 (MQ=255)
           tACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATCGaa                                                                              <  1:2876542/93‑1 (MQ=255)
             cTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGa                                                                            >  1:4028525/1‑93 (MQ=255)
             cTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGa                                                                            >  2:622450/1‑93 (MQ=255)
              tGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGAc                                                                           >  2:5128674/1‑93 (MQ=255)
                gggTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACaa                                                                         >  1:1547237/1‑93 (MQ=255)
                   tAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAaaga                                                                      <  1:4910479/93‑1 (MQ=255)
                    aaCCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAaagaa                                                                     >  1:225004/1‑93 (MQ=255)
                    aaCCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAaagaa                                                                     <  2:4924399/93‑1 (MQ=255)
                     aCCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAaagaag                                                                    >  1:2537554/1‑93 (MQ=255)
                     aCCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAaagaag                                                                    <  2:2907058/93‑1 (MQ=255)
                       cTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTc                                                                  <  1:5105720/93‑1 (MQ=255)
                        tGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTca                                                                 >  1:3494539/1‑93 (MQ=255)
                          gTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTcaca                                                               <  2:5876992/93‑1 (MQ=255)
                          gTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTcaca                                                               >  1:652312/1‑93 (MQ=255)
                          gTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTcaca                                                               >  1:3709821/1‑93 (MQ=255)
                            aCTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATATATTGAAGACAAAGAAGTCACAtt                                                             >  1:693672/1‑93 (MQ=255)
                            aCTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACAtt                                                             <  1:2502481/93‑1 (MQ=255)
                             cTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTg                                                            >  2:1190698/1‑93 (MQ=255)
                              tgttgtCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGa                                                           >  1:2473781/1‑93 (MQ=255)
                              tgttgtCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGa                                                           >  1:1675659/1‑93 (MQ=255)
                              tgtggtCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGa                                                           <  1:5958518/93‑1 (MQ=255)
                                ttgtCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACa                                                         <  2:703802/93‑1 (MQ=255)
                                  gtCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACaaa                                                       <  2:483212/93‑1 (MQ=255)
                                  gtCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACaaa                                                       >  2:5605853/1‑93 (MQ=255)
                                       ttttATCATCGGGCAGGAAAAGCCTGTAGCCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTaa                                                  <  1:1788973/93‑1 (MQ=255)
                                       ttttAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTaa                                                  >  1:2862591/1‑93 (MQ=255)
                                        tttAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAg                                                 <  1:1139184/93‑1 (MQ=255)
                                         ttAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGt                                                >  1:5069911/1‑93 (MQ=255)
                                         ttAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGt                                                <  2:2583139/93‑1 (MQ=255)
                                         ttACCATCGGACAGGAAAAGCCGATAACCTGGAGAGCCTTTCTTGAGCCGGGAAGACATTGAAGACAAAGAAGTCACAATAACAAAGTTAAGt                                                <  2:3838535/93‑1 (MQ=255)
                                          tAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTa                                               <  1:3879607/93‑1 (MQ=255)
                                          tAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTa                                               <  2:1521246/93‑1 (MQ=255)
                                           aGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAg                                              >  1:774291/1‑93 (MQ=255)
                                           aGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAg                                              <  1:5704067/93‑1 (MQ=255)
                                           aGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAg                                              <  2:4354487/93‑1 (MQ=255)
                                            ccATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGc                                             <  2:2326858/92‑1 (MQ=255)
                                              aTCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCgg                                           <  1:5389466/93‑1 (MQ=255)
                                               tCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGgc                                          <  1:172363/93‑1 (MQ=255)
                                               tCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGgc                                          >  1:434377/1‑93 (MQ=255)
                                                cGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGgcc                                         <  1:1969461/93‑1 (MQ=255)
                                                  ggCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGgccgc                                       >  1:4500574/1‑93 (MQ=255)
                                                   gCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGgccgcc                                      <  1:2843421/93‑1 (MQ=255)
                                                    cAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGgccgccg                                     >  2:2369268/1‑93 (MQ=255)
                                                       gAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCt                                  <  1:4596900/93‑1 (MQ=255)
                                                          aaGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCt                               >  1:1201335/1‑93 (MQ=255)
                                                          aaGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCt                               <  1:1323120/93‑1 (MQ=255)
                                                          aaGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCt                               >  1:6275715/1‑93 (MQ=255)
                                                           aGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTg                              >  2:1386171/1‑93 (MQ=255)
                                                            gCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTgg                             >  1:3661012/1‑93 (MQ=255)
                                                             ccTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGa                            <  1:3202470/93‑1 (MQ=255)
                                                             ccTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGa                            <  2:418290/93‑1 (MQ=255)
                                                             ccTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGa                            <  2:1874517/93‑1 (MQ=255)
                                                             ccTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGa                            <  2:1273703/93‑1 (MQ=255)
                                                               tGTACCCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAgg                          <  1:6136052/93‑1 (MQ=255)
                                                                  aaCCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGt                       >  2:1691218/1‑93 (MQ=255)
                                                                  aaCCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGt                       >  1:4772263/1‑93 (MQ=255)
                                                                  aaCCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGt                       <  1:4815098/93‑1 (MQ=255)
                                                                  aaCCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGt                       <  1:5477458/93‑1 (MQ=255)
                                                                    ccTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTg                     >  2:3784869/1‑93 (MQ=255)
                                                                     cTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGc                    >  1:6222680/1‑93 (MQ=255)
                                                                     cTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGc                    >  2:5410064/1‑93 (MQ=255)
                                                                     cTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGc                    <  2:6038052/93‑1 (MQ=255)
                                                                     cTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGc                    <  1:887201/93‑1 (MQ=255)
                                                                      tGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCt                   >  1:296499/1‑93 (MQ=255)
                                                                       ggAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCTg                  >  1:5868966/1‑93 (MQ=255)
                                                                        gagaGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCTGa                 <  2:1213358/93‑1 (MQ=255)
                                                                        gagaGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCTGa                 <  1:1009084/93‑1 (MQ=255)
                                                                        gagaGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCTGa                 <  1:3164006/93‑1 (MQ=255)
                                                                        gagaGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCTGa                 >  2:3772079/1‑93 (MQ=255)
                                                                        gagaGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCTGa                 >  2:5022958/1‑93 (MQ=255)
                                                                         agagCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCTGAt                <  2:4021039/93‑1 (MQ=255)
                                                                          gagCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCTGCCGCCGCCTTCTGGAGGCGTTGCTGATc               >  2:2703018/1‑93 (MQ=255)
                                                                          gagCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCTGATc               >  1:2103565/1‑93 (MQ=255)
                                                                           agCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCTGATcc              <  1:2790364/93‑1 (MQ=255)
                                                                            gCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCTGATCCt             >  2:3904635/1‑93 (MQ=255)
                                                                              cTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCTGATCCTTa           <  1:4275657/93‑1 (MQ=255)
                                                                               tttCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCTGATCCTTAt          <  2:5193295/93‑1 (MQ=255)
                                                                               tttCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCTGATCCTTAt          <  2:4654702/93‑1 (MQ=255)
                                                                               tttCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCTGATCCTTAt          >  2:5430186/1‑93 (MQ=255)
                                                                                 tCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCTGATCCTTATTg        >  2:5437111/1‑93 (MQ=255)
                                                                                     gAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCTGATCCTTATTGTCCt    <  2:3337572/93‑1 (MQ=255)
                                                                                     gAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCTGATCCTTATTGTCCt    <  1:4375511/93‑1 (MQ=255)
                                                                                      aGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCTGATCCTTATTGTCCTg   >  2:5699326/1‑93 (MQ=255)
                                                                                      aGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCTGATCCTTATTGTCCTg   >  1:2378582/1‑93 (MQ=255)
                                                                                       gCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCTGATCCTTATTGTCCTGt  <  2:2240585/93‑1 (MQ=255)
                                                                                          |                                                                                         
TACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCTGATCCTTATTGTCCTGT  >  NC_000913/932364‑932543

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