Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,297,884 A→G 6.1% Y495H (TAC→CAC)  fdhF ← formate dehydrogenase‑H, selenopolypeptide subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,297,8840AG 6.1% 192.3 / 2.4 65Y495H (TAC→CAC) fdhFformate dehydrogenase‑H, selenopolypeptide subunit
Reads supporting (aligned to +/- strand):  ref base A (24/37);  new base G (2/2);  total (26/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTCGGTTCAACCGCCTTGAAGAAACGCTGGAAGCCACGGTCAGCCGCAGTAAACACGCCTTCATGCTCGCCCCACGA  >  NC_000913/4297736‑4298028
                                                                                                                                                    |                                                                                                                                                
aaCAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGCTGTTGTTGTa                                                                                                                                                  >  2:744116/1‑149 (MQ=255)
 aCAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATAGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATACGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTat                                                                                                                                                 >  2:564875/1‑148 (MQ=255)
   aGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTg                                                                                                                                               >  1:535661/1‑149 (MQ=255)
     aTAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGCTGTTGTTGTAGGGCa                                                                                                                                             >  2:921580/1‑149 (MQ=255)
       aagaagTCCCCTGATCGGCATCTGAAGTATAGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATc                                                                                                                                           <  2:1091903/149‑1 (MQ=255)
        cgaagTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTCGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCg                                                                                                                                          <  2:1121952/148‑1 (MQ=255)
         gaagTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCgg                                                                                                                                         >  2:255142/1‑149 (MQ=255)
             tCCCCTGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATaa                                                                                                                                     <  2:214896/149‑1 (MQ=255)
             tCCCCTGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGGTGTTGTAGTGCATCGGATaa                                                                                                                                     >  1:560569/1‑149 (MQ=255)
              ccccTGCTCGGCATCTGAAGTATCGCGGCAAGGCCCCTGAATGAAGCCCAGTTCGCCCATGTTCTCGTAAGTCGCACCGTAGAAATCCGGGCCCAGATGACGCAACTCATCCCAGATGTCCTGGGGGTTGTTGTAGTGCAGCGGATAAc                                                                                                                                    >  1:738831/1‑149 (MQ=255)
                  tGATCGGCATCTGAAGTATCGCCGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTACAAATCCGTGCACAGATGACGCAACTAATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATATCCCAt                                                                                                                                <  1:105214/149‑1 (MQ=255)
                      cGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACgg                                                                                                                            <  1:919956/149‑1 (MQ=255)
                          ttCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCTCAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTgg                                                                                                                        <  2:1042821/148‑1 (MQ=255)
                          aTCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTgg                                                                                                                        <  1:329743/149‑1 (MQ=255)
                          aTCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTGGTTGTAGTGCATCGGATAACCCATACGGGTgg                                                                                                                        <  2:125971/149‑1 (MQ=255)
                           tCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGc                                                                                                                       >  1:1072340/1‑149 (MQ=255)
                             tGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGa                                                                                                                     <  1:107807/149‑1 (MQ=255)
                               aaGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGCTCTCCTGGGGGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGagt                                                                                                                   >  1:21901/1‑147 (MQ=255)
                                    tCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCCCAGCTGACGCCACTCATCCGAGCTCTCCTTCTTTTCTTTGTAGTGCATCGGAGACACGAACCGGGGGGCGATTTcccg                                                                                                              >  2:1186564/1‑146 (MQ=255)
                                      gcgGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACtga                                                                                                            <  2:551454/149‑1 (MQ=255)
                                       cgGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATATCCCATACGGGTGGCGATTTCACtgat                                                                                                           <  2:135189/149‑1 (MQ=255)
                                         ccAAGGCCACTGTATGAAGCCCAGATCGCCCTTTTTCTCGTAAGTCGCACCGTAGAAATGCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACtgatga                                                                                                         <  2:299338/148‑1 (MQ=255)
                                            tggCCACTGAATGAAGCCCAGTTCGACCATTTTCTCTTAAGTCGCACCGTAGAAATACGGGTACAGATGACGCAACTCATCCCAGATCTCCTGTGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGAttt                                                                                                      <  2:244534/148‑1 (MQ=255)
                                            aGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGAttt                                                                                                      <  1:858542/149‑1 (MQ=255)
                                              tccACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCTTATAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGc                                                                                                    <  2:1064910/148‑1 (MQ=255)
                                              gCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGc                                                                                                    >  2:617668/1‑149 (MQ=255)
                                                   tgaatgaaGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTc                                                                                               >  1:589401/1‑149 (MQ=255)
                                                      atgaaGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTGGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGt                                                                                            >  1:681336/1‑149 (MQ=255)
                                                      atgaaGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGGTGTTGTGGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGt                                                                                            >  1:50539/1‑149 (MQ=255)
                                                         caGCCCAGTTCGCCCATTTTCTCGTCAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGgttt                                                                                         <  2:583038/148‑1 (MQ=255)
                                                         aaGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACACATGACGCAACTCATCCCAGATCTCCAGGGTGTTGTTCTAGTGCATCGGAAAACCCATACGGGTGCCGATTTCACTGATCAATTGCCAGTCCGtatt                                                                                         >  1:7236/1‑149 (MQ=255)
                                                          aGCCCAGTTCGCCCATTTTCACGTAAGTCGCACCGGAGAAAACCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTc                                                                                        <  2:778153/149‑1 (MQ=255)
                                                            cccAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCCACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGTGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAg                                                                                      >  2:1065462/1‑149 (MQ=255)
                                                            cccAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGCTCTCCTGGGTGTTGGTGTAGTCCATCGGATAACCCATACGGGTGGCGATTTCACTGATGTTTTGCCAGTCCGGTTTCat                                                                                      >  1:240091/1‑148 (MQ=255)
                                                                gTTCGCCCATTATCTCGACAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTTCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATcc                                                                                  <  2:61191/149‑1 (MQ=255)
                                                                 ttcggacATTTTCTCGTAAGTCTCACCGTAGTTATACGGGCACAGATGACGCATCTCATCCCAGATCTCTTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATccc                                                                                 <  1:921580/143‑1 (MQ=255)
                                                                    gCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCAtt                                                                              <  1:164387/149‑1 (MQ=255)
                                                                       cATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTCg                                                                           <  2:457248/149‑1 (MQ=255)
                                                                        aTTATCTCGTAAGTCGGACCGTAGAAAGCCGGGCACAGATGACCCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCTGATAACCCATACGGGTGGCGATTTCACTGATGATTAGCCAGTCCGTTTTCAGATCCCATTTCgg                                                                          <  1:699671/149‑1 (MQ=255)
                                                                           ttCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTCGGTTc                                                                       <  1:578980/149‑1 (MQ=255)
                                                                                 tAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTCGGTTCAACCGc                                                                 >  2:389694/1‑149 (MQ=255)
                                                                                    ggcGCACCGTATAAATCCGGGCACAGATTACGCAAATCATCCCAGATCTCCTGGGTGTTGTTGTAGGGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTCGGTTCAACCGCCtt                                                              <  2:417213/147‑1 (MQ=255)
                                                                                    gTCGCACCGTCGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTCGGTTCAACCGCCtt                                                              <  2:422264/149‑1 (MQ=255)
                                                                                       gCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGAGCTCCTGGGGGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTCGGTTCAACCGCCTTgaa                                                           <  2:574392/149‑1 (MQ=255)
                                                                                           cGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTCGGTTCAACCGCCTTGAAGaaa                                                       >  1:157031/1‑149 (MQ=255)
                                                                                            gTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGGGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTCGGTTCAACCGCCTTGAAGAAAc                                                      <  1:528479/149‑1 (MQ=255)
                                                                                             tAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTCTGTTCAACCGCCTTGAAGAAACg                                                     >  1:404512/1‑149 (MQ=255)
                                                                                             tAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGGGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTCGGTTCAACCGCCTTGAAGAAACg                                                     <  2:854377/149‑1 (MQ=255)
                                                                                              aGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGGGGCGATTTCACTGATGAGTTGCCAGTCCGTTTTCAGAACCCATTTCGGTTCAACCGCCTTGAAGAAACGc                                                    >  2:972534/1‑149 (MQ=255)
                                                                                               gAAATCCGTGCACAGATGACGCAACTCATCCAAGATCTCCTGGTTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTCGGTTCAACCGCCTTGAAGAAACGCt                                                   <  1:564875/149‑1 (MQ=255)
                                                                                                aaaTCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTCGGTTCAACCGCCTTGAAGAAACGCTg                                                  <  1:342285/149‑1 (MQ=255)
                                                                                                    cctggCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTTGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTCGGTTCAACCGCCTTGAAGAAACGCTGGAAg                                              <  2:307813/146‑1 (MQ=255)
                                                                                                        gCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTCGGTTCAACCGCCTTGAAGAAACGCTGGAAGCCAc                                          >  1:870381/1‑149 (MQ=255)
                                                                                                         cacaGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTCGGTTCAACCGCCTTGAAGAAACGCTGGAAGCCACg                                         <  1:953163/149‑1 (MQ=255)
                                                                                                         cacaGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTCGGTTCAACCGCCTTGAAGAAACGCTGGAAGCCACg                                         <  2:871820/149‑1 (MQ=255)
                                                                                                            aGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTCGGTTCAACCGCCTTGAAGAACCGCTGGAAGCGACGGTc                                      >  1:258314/1‑149 (MQ=255)
                                                                                                                 aCGCATCTCATCCCAGATCTCCTGTGTGTTGTTTTAGTGCATCGGATAACGCATACGGGTGGCGATTTCACTGATTATTTGCCAGTCCGTTTTCAGATCCCATTTCGGTTCAACCGCCTTGAAGAAACGCTGGAAGCCACTGTCAGCCg                                 <  1:744116/149‑1 (MQ=255)
                                                                                                                   gCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTCGGTTCAACCGCCTTGAAGAAACGCTGGAAGCCACGGTCAGCCGCa                               >  2:934691/1‑149 (MQ=255)
                                                                                                                     aaCTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGCTCCCATTTCGGTTCAACCGCCTTGAAGAAACGCTGGAAGCCACGGTCAGCCGCAGt                             >  2:1145306/1‑149 (MQ=255)
                                                                                                                             ccAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGGCCGTTTTCAGATCCCATTTCGGTTCAACCGCCTTGAAGAAACGCTGGAAGCCACGGTCAGCCGCAGTAAACACGc                     <  1:534491/149‑1 (MQ=255)
                                                                                                                                gATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTCGGTTCAACCGCCTTGAAGAAACGCTGGAAGCCACGGTCAGCCGCAGTAAACACGCCtt                  >  2:682477/1‑149 (MQ=255)
                                                                                                                                gATCTCCTGGGTGTTGGTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTCGGTTCAACCGCCTTGACGAAACGCTGGAAGCCACGGTCAGCCGCAGTAAACACGCCtt                  <  1:255142/149‑1 (MQ=255)
                                                                                                                                 aTCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCGATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTCGGTTCAACCGCCTTGAAGAAACGCTGGAAGCCACGGTCAGCCGCAGTAAACACGCCTTc                 <  1:1114088/149‑1 (MQ=255)
                                                                                                                                           tgttgttgtGGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTCGGTTCAACCGCCTTGAAGAAACGCTGGAAGCCACGGTCAGCCGCAGTAAACACGCCTTCATGCTCGccc       <  2:50539/149‑1 (MQ=255)
                                                                                                                                           tgttgttggGGTGCAGCGGATAACCCCTACGGGTGGCGATTTCACTGATGATTTGCCAGGCCGTTTTCAGATCCCATTGCGGTTCAACCGCCTTGAAGAAACGCTGGAAGCCACGGTCAGCCGCAGTAAACACGCCTTCATGCTCGccc       <  2:681336/149‑1 (MQ=255)
                                                                                                                                              tgttgtAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTCGGTTCAACCTCCTTGAAGAAACGCGGGAAGCCACGGTCAGCCGCAGTAAACACGCCGTCATTCGCGCCCCac    >  2:289401/1‑149 (MQ=255)
                                                                                                                                                ttgtAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTCGGTTCAACCGCCTTGAAGAAACGCTGGAAGCCACGGTCAGCCGCAGTAAACACGCCTTCATGCTCGCCCCacga  >  1:715288/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
AACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTCGGTTCAACCGCCTTGAAGAAACGCTGGAAGCCACGGTCAGCCGCAGTAAACACGCCTTCATGCTCGCCCCACGA  >  NC_000913/4297736‑4298028

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑