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,4170GA97.6% 574.9 / ‑0.5 165P510L (CCG→CTG) fdhFselenopolypeptide subunit of formate dehydrogenase H
Reads supporting (aligned to +/- strand):  ref base G (1/3);  new base A (79/82);  total (80/85)
Fisher's exact test for biased strand distribution p-value = 6.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 171 total aligned reads displayed.

GTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATT  >  NC_000913/4295327‑4295506
                                                                                          |                                                                                         
gTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCagg                                                                                         >  2:3266443/1‑93 (MQ=255)
gTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCACCCATTTTCTCGTAAGTCGCACCGTAGAAATCCagg                                                                                         <  2:1366718/93‑1 (MQ=255)
 gccccTGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGc                                                                                        <  1:4076828/92‑1 (MQ=255)
  ccccTGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGca                                                                                       <  2:700452/93‑1 (MQ=255)
   cccTGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGcac                                                                                      >  1:901377/1‑93 (MQ=255)
     cTGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCCCAg                                                                                    >  1:4404241/1‑93 (MQ=255)
      tGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGa                                                                                   >  1:478729/1‑93 (MQ=255)
      tGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGa                                                                                   >  2:3172686/1‑93 (MQ=255)
      tGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGa                                                                                   >  1:566897/1‑93 (MQ=255)
       gATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGAt                                                                                  >  2:203873/1‑93 (MQ=255)
       gATCGCCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCCCCGTAGAAATCCGGGCACAGAt                                                                                  <  2:3087655/93‑1 (MQ=255)
        aTCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATg                                                                                 >  2:3798025/1‑93 (MQ=255)
        aTCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATg                                                                                 <  1:4381826/93‑1 (MQ=255)
         tCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGa                                                                                <  1:283614/93‑1 (MQ=255)
           ggCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACg                                                                              >  1:1872686/1‑93 (MQ=255)
            gCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGc                                                                             <  1:3000245/93‑1 (MQ=255)
               tCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAAc                                                                          <  2:1202799/93‑1 (MQ=255)
               tCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGAAGCAAc                                                                          >  1:726878/1‑93 (MQ=255)
                 tGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTc                                                                        >  2:1117053/1‑93 (MQ=255)
                 tGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTc                                                                        >  2:1500449/1‑93 (MQ=255)
                  gAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCa                                                                       >  2:3546814/1‑93 (MQ=255)
                   aaGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCAt                                                                      >  1:4324061/1‑93 (MQ=255)
                   aaGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCAt                                                                      <  2:2531955/93‑1 (MQ=255)
                   aaGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCAt                                                                      >  1:4741614/1‑93 (MQ=255)
                   aaGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCAt                                                                      >  2:3440155/1‑93 (MQ=255)
                    aGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATc                                                                     <  1:1202058/93‑1 (MQ=255)
                     gTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATcc                                                                    <  2:213874/93‑1 (MQ=255)
                     gTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATcc                                                                    <  2:2740954/93‑1 (MQ=255)
                      gaTCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATccc                                                                   <  1:1242251/92‑1 (MQ=255)
                         cgcgGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGa                                                                <  2:2324391/93‑1 (MQ=255)
                         cgcgGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGa                                                                <  1:1436477/93‑1 (MQ=255)
                          gcgGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGAt                                                               <  2:1078236/93‑1 (MQ=255)
                           cgGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGAtc                                                              <  1:199862/93‑1 (MQ=255)
                             gCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGAtctc                                                            <  1:2710972/93‑1 (MQ=255)
                              cccggCCCCTGGCTGAAGCCCCGTTCGCCCATTTTCTCGTCAGTCGCTCCGTCGCACTCCAGGCCCCGCTGCCGCAACTCATCCCAGATCTcc                                                           <  1:4197507/90‑1 (MQ=255)
                               aaGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCt                                                          >  1:2704336/1‑93 (MQ=255)
                               aaGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCt                                                          <  1:4340132/93‑1 (MQ=255)
                                aGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTg                                                         <  1:1619445/93‑1 (MQ=255)
                                aGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTg                                                         <  1:4181727/93‑1 (MQ=255)
                                 ggCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTgg                                                        >  2:4703312/1‑93 (MQ=255)
                                  gCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTggg                                                       >  2:1573714/1‑93 (MQ=255)
                                  gCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTggg                                                       >  1:4530751/1‑93 (MQ=255)
                                  gCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTggg                                                       >  1:2167256/1‑93 (MQ=255)
                                    cACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGgtg                                                     <  1:58260/93‑1 (MQ=255)
                                     aCTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGgtgt                                                    >  1:2941258/1‑93 (MQ=255)
                                     aCTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGgtgt                                                    <  1:4702080/93‑1 (MQ=255)
                                     aCTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGgtgt                                                    <  1:1843118/93‑1 (MQ=255)
                                      cTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGtgtt                                                   >  2:3582758/1‑93 (MQ=255)
                                      cTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGtgtt                                                   >  2:4562614/1‑93 (MQ=255)
                                      cTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGtgtt                                                   >  1:3595899/1‑93 (MQ=255)
                                       tgaatgaaGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGtgttg                                                  >  2:4623968/1‑93 (MQ=255)
                                            gaaGTCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTa                                             <  1:1494184/93‑1 (MQ=255)
                                               gCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTg                                          >  1:3207283/1‑93 (MQ=255)
                                                 ccAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCa                                        >  2:2965252/1‑93 (MQ=255)
                                                 ccAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCa                                        <  2:4460953/93‑1 (MQ=255)
                                                 ccAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCa                                        <  2:1279110/93‑1 (MQ=255)
                                                  cAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCAt                                       <  1:2264842/93‑1 (MQ=255)
                                                   aGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATc                                      >  2:607869/1‑93 (MQ=255)
                                                    gTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCg                                     <  1:1523357/93‑1 (MQ=255)
                                                    gTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCg                                     <  1:538384/93‑1 (MQ=255)
                                                     ttCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCgg                                    >  1:3554409/1‑93 (MQ=255)
                                                     ttCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCgg                                    <  2:2233820/93‑1 (MQ=255)
                                                     ttCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCgg                                    <  1:965703/93‑1 (MQ=255)
                                                      tCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGa                                   <  1:737251/93‑1 (MQ=255)
                                                      tCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGa                                   <  2:4364740/93‑1 (MQ=255)
                                                      tCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGa                                   <  2:4228028/93‑1 (MQ=255)
                                                      tCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGa                                   <  2:4437850/93‑1 (MQ=255)
                                                      tCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGa                                   <  1:3651083/93‑1 (MQ=255)
                                                       cGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGAt                                  <  1:3442445/93‑1 (MQ=255)
                                                          ccATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAAc                               >  2:3030021/1‑93 (MQ=255)
                                                              tttCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCAt                           >  1:1047705/1‑93 (MQ=255)
                                                               ttCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATa                          <  1:389377/93‑1 (MQ=255)
                                                                tctcGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATAc                         >  1:3214902/1‑93 (MQ=255)
                                                                tctcGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATAc                         <  1:2362251/93‑1 (MQ=255)
                                                                tctcGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATAc                         <  2:41477/93‑1 (MQ=255)
                                                                 ctcGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACg                        <  2:457606/93‑1 (MQ=255)
                                                                  tcGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACgg                       >  2:2753095/1‑93 (MQ=255)
                                                                   cGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACggg                      <  1:3864991/93‑1 (MQ=255)
                                                                    gTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGt                     <  1:2520408/93‑1 (MQ=255)
                                                                    gTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGt                     >  2:4390488/1‑93 (MQ=255)
                                                                      aaGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTgg                   >  2:2462681/1‑93 (MQ=255)
                                                                      aaGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTgg                   >  1:166542/1‑93 (MQ=255)
                                                                         tCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGa                <  1:503401/93‑1 (MQ=255)
                                                                             cccGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTc            <  1:1281302/92‑1 (MQ=255)
                                                                             aCCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTc            <  1:431856/93‑1 (MQ=255)
                                                                              ccGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCa           >  1:3085120/1‑93 (MQ=255)
                                                                              ccGTAAAAAGCCAGGCCCAGAGGAACCAGCTCATCCCAGATCTCCTGGGTGTTGTTGTGGTGCATCGGATAACCCATACGGGTGGCGATTTCa           <  1:4501640/93‑1 (MQ=255)
                                                                               cGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCAc          >  1:25399/1‑93 (MQ=255)
                                                                               cGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCAc          >  1:25400/1‑93 (MQ=255)
                                                                               cGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCAc          <  1:3545282/93‑1 (MQ=255)
                                                                               cGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCAc          <  1:3910892/93‑1 (MQ=255)
                                                                                 tAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGCTTTCACtg        >  2:915919/1‑93 (MQ=255)
                                                                                 tAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACtg        <  1:2954788/93‑1 (MQ=255)
                                                                                  aGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACtga       <  2:4060850/93‑1 (MQ=255)
                                                                                  aGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGAGGCGATTTCACtga       >  1:1302569/1‑93 (MQ=255)
                                                                                   gAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACtgat      >  1:2612143/1‑93 (MQ=255)
                                                                                   gAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACtgat      <  1:4466185/93‑1 (MQ=255)
                                                                                    aaaTCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACtgatg     <  1:3266443/93‑1 (MQ=255)
                                                                                       tCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGAtt  <  1:3363101/93‑1 (MQ=255)
                                                                                       accGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGAtt  <  1:3995117/92‑1 (MQ=255)
                                                                                          |                                                                                         
GTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATT  >  NC_000913/4295327‑4295506

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