Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,021,684 Δ1 bp coding (542/609 nt)|coding (1/243 nt) b4556 →
tatD →
hypothetical protein
Mg‑dependent DNase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,021,6830C.100.0% 866.5 / NA 190coding (541/609 nt)b4556hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (97/93);  total (97/93)

Only 100 of 205 total aligned reads displayed.

TCTTCATTGCTTTACCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACCGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTCCCT  >  NC_000913/4021592‑4021770
                                                                                            |                                                                                       
tcttcATTGCTTTACCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAAcg                                                                                        >  2:172103/1‑92 (MQ=255)
 cttcATTGCTTTACCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAAcga                                                                                       <  2:4460882/93‑3 (MQ=255)
 cttcATTGCTTTACCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAAcga                                                                                       <  1:470758/93‑3 (MQ=255)
   tcATTGCTTTACCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAAcgacg                                                                                     >  2:1517654/1‑89 (MQ=255)
    cATTGCTTTACCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGAcgc                                                                                   <  1:2292627/93‑1 (MQ=255)
     aTTGCTTTACCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGTACGAcgcg                                                                                  >  2:2342379/1‑93 (MQ=255)
        gCTTTACCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGACGAACGACGCGGAc                                                                               >  1:1956815/1‑93 (MQ=255)
         cTTTACCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACt                                                                              >  2:2812012/1‑93 (MQ=255)
          tttACCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTg                                                                             >  1:1277919/1‑93 (MQ=255)
          tttACCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTg                                                                             >  1:1448913/1‑93 (MQ=255)
          tttACCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTg                                                                             <  2:33163/93‑1 (MQ=255)
             aCCGGCACACGCGAAGAGATGCAGGCGAGCGTGGCGCATGGAATTTATATAGGCATTACCGGTTAGGTTTGCGATGCACGACGCGGACTggtg                                                                          >  2:463558/1‑91 (MQ=255)
              ccGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGc                                                                         <  1:1947128/93‑1 (MQ=255)
              ccGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGc                                                                         <  1:7462/93‑1 (MQ=255)
              ccGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGc                                                                         <  2:1068139/93‑1 (MQ=255)
              ccGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGc                                                                         <  2:2583704/93‑1 (MQ=255)
              ccGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGc                                                                         >  2:3444043/1‑93 (MQ=255)
                ggCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTg                                                                       >  1:1715933/1‑93 (MQ=255)
                ggCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTg                                                                       >  1:2630531/1‑93 (MQ=255)
                 gCAGACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGc                                                                      <  1:964301/93‑1 (MQ=255)
                 gCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGc                                                                      <  2:2666323/93‑1 (MQ=255)
                 gCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGc                                                                      <  1:1818370/93‑1 (MQ=255)
                 gCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGc                                                                      >  2:359683/1‑93 (MQ=255)
                 gCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGc                                                                      >  1:2558695/1‑93 (MQ=255)
                  cacaCGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCg                                                                     >  1:1224803/1‑93 (MQ=255)
                   acacGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCgg                                                                    >  2:4393169/1‑93 (MQ=255)
                   acacGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCgg                                                                    >  1:4191563/1‑93 (MQ=255)
                    cacGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCggg                                                                   >  2:3138147/1‑93 (MQ=255)
                        cgAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACt                                                               <  2:917588/93‑1 (MQ=255)
                         gAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACtt                                                              <  1:2030221/93‑1 (MQ=255)
                          aaGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACttt                                                             <  1:2730055/93‑1 (MQ=255)
                           agagATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACtttt                                                            <  2:1295947/93‑1 (MQ=255)
                            gagaTGCAGGCGTGCGTGGCGCATGGAATTTATGTCGGCATTACGGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTg                                                           <  2:1935775/93‑1 (MQ=255)
                            gagaTGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTAACCGGTTGGGTTGGGGTGGACGACAGGGACTGGGGGCTCGGGAACttttc                                                           >  2:3073939/1‑92 (MQ=255)
                              gaTGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGGACTTTTGcc                                                         >  2:539402/1‑93 (MQ=255)
                              gaTGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGcc                                                         <  1:161252/93‑1 (MQ=255)
                                tGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGt                                                       <  2:609231/93‑1 (MQ=255)
                                  cAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTg                                                     >  1:2589909/1‑93 (MQ=255)
                                   aGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGa                                                    >  1:4486560/1‑93 (MQ=255)
                                    ggCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGAt                                                   <  1:4600785/93‑1 (MQ=255)
                                      cgtgcgtgGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTc                                                 <  2:2353712/93‑1 (MQ=255)
                                       gtgcgtgGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTcc                                                <  1:2246157/93‑1 (MQ=255)
                                       gtgcgtgGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTcc                                                >  2:3055302/1‑93 (MQ=255)
                                         gcgtgGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCcgg                                              >  2:3382601/1‑93 (MQ=255)
                                            tgGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCcggcgg                                           >  2:1769750/1‑93 (MQ=255)
                                             ggCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGTCGTTGATTCCGGCGGa                                          >  1:2176352/1‑93 (MQ=255)
                                             ggCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGa                                          >  1:2286069/1‑93 (MQ=255)
                                                gcATGGAATTTATATCGGCATTACGGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGaaaa                                       <  2:3516245/93‑1 (MQ=255)
                                                gcATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGaaaa                                       >  1:100884/1‑93 (MQ=255)
                                                gcATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGaaaa                                       >  1:1999522/1‑93 (MQ=255)
                                                gcATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGaaaa                                       >  1:3138506/1‑93 (MQ=255)
                                                 cATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGaaaaa                                      >  1:1303757/1‑93 (MQ=255)
                                                  aTGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAAt                                     >  2:3488014/1‑93 (MQ=255)
                                                  aTGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAAt                                     >  1:1732321/1‑93 (MQ=255)
                                                  aTGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAAt                                     >  1:3865234/1‑93 (MQ=255)
                                                  aTGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAAt                                     >  1:2375077/1‑93 (MQ=255)
                                                  aTGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAAt                                     <  2:2053396/93‑1 (MQ=255)
                                                   tGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAAtt                                    >  2:1396403/1‑93 (MQ=255)
                                                    ggAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGTAAAATTa                                   >  2:4417417/1‑93 (MQ=255)
                                                    ggAATTTAGATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTa                                   <  2:2274706/93‑1 (MQ=255)
                                                     gAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTAc                                  <  1:834736/93‑1 (MQ=255)
                                                     gAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTAc                                  >  1:1067769/1‑93 (MQ=255)
                                                     gAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTAc                                  >  2:263415/1‑93 (MQ=255)
                                                       aTTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTg                                >  2:3208578/1‑93 (MQ=255)
                                                       aTTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTg                                >  1:4322995/1‑93 (MQ=255)
                                                       aTTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTg                                <  1:3453739/93‑1 (MQ=255)
                                                        tttATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGa                               >  2:122746/1‑93 (MQ=255)
                                                          tataTCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATc                             <  2:1739366/93‑1 (MQ=255)
                                                          tataTCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATc                             <  1:4256245/93‑1 (MQ=255)
                                                            tatCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGa                           >  1:3161007/1‑93 (MQ=255)
                                                            tatCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGa                           <  1:1531124/93‑1 (MQ=255)
                                                              tCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGaaa                         <  2:767717/93‑1 (MQ=255)
                                                               cGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAAc                        <  1:3395077/93‑1 (MQ=255)
                                                                  cATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGa                     <  1:3149247/93‑1 (MQ=255)
                                                                   aTTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGAt                    <  2:1476670/93‑1 (MQ=255)
                                                                   aTTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGAt                    >  1:3787526/1‑93 (MQ=255)
                                                                   aTTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGAt                    <  2:1111308/93‑1 (MQ=255)
                                                                   aTTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGAt                    >  2:1898921/1‑93 (MQ=255)
                                                                     tACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATgg                  >  1:3409525/1‑92 (MQ=255)
                                                                     tACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATgc                  <  2:3759878/93‑1 (MQ=255)
                                                                     tACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATgc                  <  1:1123824/93‑1 (MQ=255)
                                                                         cGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGcc               <  2:2105174/93‑1 (MQ=255)
                                                                          ggTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCg              >  1:3948232/1‑93 (MQ=255)
                                                                          ggTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCg              >  1:1745720/1‑93 (MQ=255)
                                                                           gTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGt             >  2:819180/1‑93 (MQ=255)
                                                                            ttGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTa            <  1:561938/93‑1 (MQ=255)
                                                                            ttGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTa            <  1:2310526/93‑1 (MQ=255)
                                                                              gggTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATc          <  1:3709984/93‑1 (MQ=255)
                                                                              gggTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATc          <  1:2137141/93‑1 (MQ=255)
                                                                              gggTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATc          <  2:3959174/93‑1 (MQ=255)
                                                                              gggTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATc          >  1:606541/1‑93 (MQ=255)
                                                                               ggTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCt         >  1:2604195/1‑93 (MQ=255)
                                                                                gTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTg        <  1:3078794/93‑1 (MQ=255)
                                                                                 tttGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGc       <  1:3006764/93‑1 (MQ=255)
                                                                                 tttGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGc       <  1:1801897/93‑1 (MQ=255)
                                                                                  ttGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCt      <  2:2126547/93‑1 (MQ=255)
                                                                                   tGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTc     <  2:4550623/93‑1 (MQ=255)
                                                                                    gCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTcc    <  2:4195980/93‑1 (MQ=255)
                                                                                     cGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTccc   >  2:2152800/1‑93 (MQ=255)
                                                                                      gATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTCCCt  >  2:4366058/1‑93 (MQ=255)
                                                                                            |                                                                                       
TCTTCATTGCTTTACCGGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACCGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGCGCCGTATCTGCTCCCT  >  NC_000913/4021592‑4021770

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