Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,295,417 G→A P510L (CCG→CTG)  fdhF ← selenopolypeptide subunit of formate dehydrogenase H

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,295,4170GA99.7% 1060.8 / ‑5.6 285P510L (CCG→CTG) fdhFselenopolypeptide subunit of formate dehydrogenase H
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (160/124);  total (161/124)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.19e-01

Only 100 of 294 total aligned reads displayed.

GTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATT  >  NC_000913/4295327‑4295506
                                                                                          |                                                                                         
gTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCagg                                                                                         >  1:1586936/1‑93 (MQ=255)
gTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCagg                                                                                         >  2:6498839/1‑93 (MQ=255)
   cccTGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGcac                                                                                      >  1:4093699/1‑93 (MQ=255)
   cccTGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGcac                                                                                      >  2:930227/1‑93 (MQ=255)
   cccTGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGcac                                                                                      >  2:6015184/1‑93 (MQ=255)
    ccTGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGcaca                                                                                     <  2:1083616/93‑1 (MQ=255)
     cTGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAg                                                                                    >  2:5562938/1‑93 (MQ=255)
      tGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGa                                                                                   >  2:2920253/1‑93 (MQ=255)
      tGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGa                                                                                   >  2:2062633/1‑93 (MQ=255)
       gATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGAt                                                                                  <  2:6344880/93‑1 (MQ=255)
       gATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGAt                                                                                  >  2:2165360/1‑93 (MQ=255)
       gATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGAt                                                                                  <  2:3235480/93‑1 (MQ=255)
        aTCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATg                                                                                 <  2:6293903/93‑1 (MQ=255)
         tCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGa                                                                                <  1:646699/93‑1 (MQ=255)
         tCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGa                                                                                >  2:4772429/1‑93 (MQ=255)
          cGGCATCTGAAGTATCGCGGCAAAGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGAc                                                                               >  2:4348354/1‑93 (MQ=255)
           ggCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACg                                                                              >  2:5211475/1‑93 (MQ=255)
           ggCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACg                                                                              >  2:3939041/1‑93 (MQ=255)
           ggCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACg                                                                              >  1:5869851/1‑93 (MQ=255)
            gCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGc                                                                             >  1:6267140/1‑93 (MQ=255)
              aTCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCaa                                                                           >  1:5048647/1‑93 (MQ=255)
              aTCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCaa                                                                           >  1:1018891/1‑93 (MQ=255)
               tCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAAc                                                                          >  1:6107544/1‑93 (MQ=255)
               tCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAAc                                                                          <  2:67385/93‑1 (MQ=255)
                cTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACt                                                                         <  2:2505583/93‑1 (MQ=255)
                 tGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTc                                                                        >  1:4324005/1‑93 (MQ=255)
                 tGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTc                                                                        >  2:5397453/1‑93 (MQ=255)
                 tGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTc                                                                        <  1:976890/93‑1 (MQ=255)
                  gAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACGGTAGAAATCCAGGCACAGATGACGCAACTCa                                                                       <  1:926863/93‑1 (MQ=255)
                  gAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCa                                                                       <  2:2337352/93‑1 (MQ=255)
                  gAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCa                                                                       >  1:5696476/1‑93 (MQ=255)
                      tATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATccc                                                                   <  1:4953549/93‑1 (MQ=255)
                       aTCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCa                                                                  <  2:6085834/93‑1 (MQ=255)
                         cgcgGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGa                                                                <  2:3349722/93‑1 (MQ=255)
                           cgGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGAtc                                                              <  2:4259879/93‑1 (MQ=255)
                           cgGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGAtc                                                              <  1:5134477/93‑1 (MQ=255)
                           cgGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGAtc                                                              >  2:2803032/1‑93 (MQ=255)
                            ggCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGAtct                                                             >  2:6441895/1‑93 (MQ=255)
                             gCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGAtctc                                                            >  2:1221583/1‑93 (MQ=255)
                             gCAAGCCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTcc                                                           <  1:2841393/93‑1 (MQ=255)
                              cAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTcc                                                           <  1:1141948/93‑1 (MQ=255)
                              cAAGGCCACTGAATGAAGCCCAGGTCGCCCATTTTCTCGTAAGTCGTACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTcc                                                           >  1:5391863/1‑93 (MQ=255)
                                aGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTg                                                         >  1:2439000/1‑93 (MQ=255)
                                  gCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTggg                                                       <  2:5756337/93‑1 (MQ=255)
                                   ccACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGgt                                                      >  1:3999631/1‑93 (MQ=255)
                                   ccACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGgt                                                      >  2:5861501/1‑93 (MQ=255)
                                     aCTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGgtgt                                                    >  1:479927/1‑93 (MQ=255)
                                     aCTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGgtgt                                                    <  1:4569652/93‑1 (MQ=255)
                                     aCTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGgtgt                                                    >  1:1973765/1‑93 (MQ=255)
                                       tgaatgaaGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGtgttg                                                  >  2:3585115/1‑93 (MQ=255)
                                        gaatgaaGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGtgttgt                                                 <  1:5296822/93‑1 (MQ=255)
                                        gaatgaaGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGtgttgt                                                 >  2:2584614/1‑93 (MQ=255)
                                         aatgaaGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGtgttgtt                                                <  1:1719766/93‑1 (MQ=255)
                                          atgaaGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGtgttgttg                                               >  1:2399663/1‑93 (MQ=255)
                                           tgaaGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGtgttgttgt                                              >  1:3120251/1‑93 (MQ=255)
                                             aaGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAg                                            <  1:2313362/93‑1 (MQ=255)
                                             aaGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAg                                            >  1:5513381/1‑93 (MQ=255)
                                             aaGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAg                                            <  2:1583056/93‑1 (MQ=255)
                                              aGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGt                                           >  2:6160943/1‑93 (MQ=255)
                                               gCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTg                                          >  2:3198934/1‑93 (MQ=255)
                                                cccAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGc                                         <  2:3426925/93‑1 (MQ=255)
                                                cccAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGc                                         <  2:4564418/93‑1 (MQ=255)
                                                 ccAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCa                                        <  1:6697883/93‑1 (MQ=255)
                                                  cAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCAt                                       >  1:6202917/1‑93 (MQ=255)
                                                   aGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATc                                      >  2:119354/1‑93 (MQ=255)
                                                   aGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATc                                      <  1:259421/93‑1 (MQ=255)
                                                   aGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATc                                      <  2:545606/93‑1 (MQ=255)
                                                    gTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCg                                     >  1:776263/1‑93 (MQ=255)
                                                     ttCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCgg                                    >  1:1783779/1‑93 (MQ=255)
                                                       cGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGAt                                  <  2:2709486/93‑1 (MQ=255)
                                                       cGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGAt                                  <  2:5936313/93‑1 (MQ=255)
                                                        gCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATa                                 >  2:1816425/1‑93 (MQ=255)
                                                        gCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATa                                 >  2:252903/1‑93 (MQ=255)
                                                           cATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAAcc                              >  1:1260260/1‑93 (MQ=255)
                                                           cATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAAcc                              >  1:2184694/1‑93 (MQ=255)
                                                           cATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAAcc                              >  1:5347835/1‑93 (MQ=255)
                                                            aTTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCAGAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAAccc                             <  1:1413915/93‑1 (MQ=255)
                                                            aTTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAAccc                             <  1:5978337/93‑1 (MQ=255)
                                                            aTTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAAccc                             >  1:3265359/1‑93 (MQ=255)
                                                              tttCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCAt                           >  1:4125184/1‑93 (MQ=255)
                                                               ttCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATa                          <  2:4990973/93‑1 (MQ=255)
                                                                tctcGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATAc                         <  1:3155904/93‑1 (MQ=255)
                                                                tctcGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATAc                         <  1:3473697/93‑1 (MQ=255)
                                                                    gTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACgggg                     >  2:1662698/1‑92 (MQ=255)
                                                                    gTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGt                     <  1:2699214/93‑1 (MQ=255)
                                                                    gTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGt                     >  1:442944/1‑93 (MQ=255)
                                                                      aaGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTgg                   >  2:4599607/1‑93 (MQ=255)
                                                                        gTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCg                 >  2:2256279/1‑93 (MQ=255)
                                                                          cGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGAt               <  2:1431159/93‑1 (MQ=255)
                                                                           gCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGAtt              >  1:4523617/1‑93 (MQ=255)
                                                                           gCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGAtt              >  2:1540929/1‑93 (MQ=255)
                                                                           gCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGAtt              <  1:3761910/93‑1 (MQ=255)
                                                                            cACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGAttt             >  2:5253462/1‑93 (MQ=255)
                                                                             aCCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTc            >  2:2191915/1‑93 (MQ=255)
                                                                             aCCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTc            <  2:417114/93‑1 (MQ=255)
                                                                                  aGAAATCCAGGCACAGATGACGCACCTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACtga       <  2:3116598/93‑1 (MQ=255)
                                                                                  aGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACtga       >  1:3570645/1‑93 (MQ=255)
                                                                                  aGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACtga       <  2:6588644/93‑1 (MQ=255)
                                                                                   gAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACtgat      >  1:1930345/1‑93 (MQ=255)
                                                                                   gAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACtgat      <  1:3059246/93‑1 (MQ=255)
                                                                                          |                                                                                         
GTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATT  >  NC_000913/4295327‑4295506

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