Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,652,855:1 (G)6→7 100% coding (362/846 nt) sseA → 3‑mercaptopyruvate sulfurtransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,652,8491.G93.1% 483.0 / 25.4 144coding (356/846 nt)sseA3‑mercaptopyruvate sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base . (8/2);  new base G (65/69);  total (73/71)
Fisher's exact test for biased strand distribution p-value = 9.74e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 162 total aligned reads displayed.

GGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑T‑GG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGA  >  NC_000913/2652751‑2652948
                                                                                                    |                                                                                                      
ggATTAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑T‑g                                                                                                       <  2:545116/100‑1 (MQ=255)
   tAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑T‑gg‑g‑g                                                                                                  >  2:3788279/1‑100 (MQ=255)
     aGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGGGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑T‑gg‑g‑ggg                                                                                                <  2:2048099/100‑1 (MQ=255)
     aGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑T‑gg‑g‑ggg                                                                                                >  1:444384/1‑100 (MQ=255)
     aGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑T‑gg‑g‑gg                                                                                                 <  2:1744452/99‑1 (MQ=255)
      gCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑T‑gg‑g‑gggg                                                                                               >  2:3605682/1‑99 (MQ=255)
      gCAACTGATTGTCTATGACGAAGGTAATCTTATCACAGCCTCACGAGCAGGGTGGAGGCTGGGCACATGTGGTCTAGAGAAAGTGTCGATTC‑T‑gg‑g                                                                                                    >  2:810213/1‑96 (MQ=255)
        aCCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑T‑gg‑g‑ggggtg                                                                                             <  2:1234956/100‑4 (MQ=255)
         ccTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑T‑gg‑g‑ggggtgg                                                                                            <  1:3713712/100‑5 (MQ=255)
         ccTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑T‑gg‑g‑ggggtgg                                                                                            >  2:3293589/1‑96 (MQ=255)
            gATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGG‑G‑GGG‑TGGACt                                                                                         >  1:450995/1‑100 (MQ=255)
              ttGTCTATGATGAAGGTAATCTTTTCTCAGCCTCACGAGCATGGAGGATGCTGCGCACCTTTGGCGAAGAGAAAGGGGGGAGTC‑G‑GG‑G‑GGG‑TGGACTTg                                                                                       >  2:2808738/1‑99 (MQ=255)
              ttGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑T‑gg‑g‑gggggtggactt                                                                                       <  2:576850/100‑10 (MQ=255)
                gTCTATGACTAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGG‑G‑GGG‑TGGACTTGCa                                                                                     <  2:3417730/100‑1 (MQ=255)
                gTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGG‑G‑GGG‑TGGACTTGCa                                                                                     <  2:2302358/100‑1 (MQ=255)
                gTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGG‑G‑GGG‑TGGACTTGCa                                                                                     <  1:3808986/100‑1 (MQ=255)
                gTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGG‑G‑GGG‑TGGACTTGCa                                                                                     <  2:1371576/100‑1 (MQ=255)
                 tCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGG‑G‑GGG‑AGGACTTGCAg                                                                                    >  2:2481829/1‑100 (MQ=255)
                  cTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGG‑G‑GGG‑TGGACTTGCAgg                                                                                   <  2:3413487/100‑1 (MQ=255)
                  cTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGG‑G‑GGG‑TGGACTTGCAgg                                                                                   <  1:2819962/100‑1 (MQ=255)
                  cTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGG‑G‑GGG‑TGGACTTGCAgg                                                                                   <  1:2960307/100‑1 (MQ=255)
                  cTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTCGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGG‑G‑GGG‑TGGACTTGCAg                                                                                    <  2:2672939/99‑1 (MQ=255)
                   tATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGG‑G‑GGG‑TGGACTTGCAGGc                                                                                  >  1:1815456/1‑100 (MQ=255)
                     tGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGG‑G‑GGG‑TGGACTTGCAGGCTg                                                                                <  1:1276934/100‑1 (MQ=255)
                     tGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGG‑G‑GGG‑TGGACTTGCAGGCTg                                                                                >  1:1111630/1‑100 (MQ=255)
                     tGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGG‑G‑GGG‑TGGACTTGCAGGCTg                                                                                >  1:654344/1‑100 (MQ=255)
                      gACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGG‑G‑GGG‑TGGACTTGCAGGCTgg                                                                               <  1:336200/100‑1 (MQ=255)
                       aCGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGG‑G‑GGG‑TGGACTTGCAGGCTGGc                                                                              <  1:1765728/100‑1 (MQ=255)
                        cggagggtaTCGTTTCTCTGGCCCACGAGCAGTGTGGATGCTCGGTTCCTTTGGGGTGTAGAAAATGTCTTTTC‑TGGG‑G‑GGG‑TGGACTTGCACTCTggta                                                                             >  2:3625869/9‑98 (MQ=255)
                        cGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGG‑G‑GGG‑TGGACTTGCAGGCTGGCa                                                                             >  2:1510245/1‑100 (MQ=255)
                         gAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAg                                                                            <  2:2161390/100‑1 (MQ=255)
                          aaGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAgc                                                                           <  1:3018971/100‑1 (MQ=255)
                             gTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGcgcg                                                                        >  2:3074251/1‑100 (MQ=255)
                              tAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGcgcg                                                                        >  2:1317332/1‑99 (MQ=255)
                               aaTCTTTTCTCATCCCCACGAGCATGGTGGATGATACGAACCTTTGGTGTCGAGAAAGTGTCGATTC‑T‑GGCG‑GGG‑TGGACTTGCACGCTGGCAGcacg                                                                        >  2:1660256/1‑98 (MQ=255)
                               aaTCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCgat                                                                      <  1:1000567/100‑1 (MQ=255)
                                aTCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCgatg                                                                     >  2:399639/1‑100 (MQ=255)
                                aTCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCgat                                                                      >  2:1068875/1‑99 (MQ=255)
                                 tCTTTTCGCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCgatga                                                                    <  2:580051/100‑1 (MQ=255)
                                   ttttCTCAGCCCCACGAGCATGGTGGATGCTGCTCTCCTTTGGTCTAGACACAATGCCCATTC‑A‑GG‑GCCGG‑TGAACTTGCAGACTAGCAGCCCgatgcgc                                                                  >  2:3176142/1‑97 (MQ=255)
                                   ttttCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATc                                                                  >  1:152794/1‑100 (MQ=255)
                                    tttCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATct                                                                 <  1:1069598/100‑1 (MQ=255)
                                     ttCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATctg                                                                >  2:3339643/1‑100 (MQ=255)
                                      tctcAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATctgc                                                               <  1:2144375/100‑1 (MQ=255)
                                       ctcAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATctgct                                                              >  1:410759/1‑100 (MQ=255)
                                       ctcAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATctgct                                                              >  2:2853693/1‑100 (MQ=255)
                                       ctcAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGGGTAGAGAAAGTGTCGATTC‑TGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATctgct                                                              >  1:3217744/1‑100 (MQ=255)
                                        tcAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATctgctg                                                             <  1:3625576/100‑1 (MQ=255)
                                         cAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGt                                                            <  1:1464293/100‑1 (MQ=255)
                                          aGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGtt                                                           <  1:3271425/100‑1 (MQ=255)
                                          aGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGtt                                                           <  1:1317341/100‑1 (MQ=255)
                                           gctccACGGGCATGGTTGTTGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑T‑GG‑G‑GGGTTGTACTTGCAGGCTGGCAGCGCCATGATCTGCTGTTa                                                          <  2:896497/97‑1 (MQ=255)
                                           gCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTa                                                          <  2:1186264/100‑1 (MQ=255)
                                            ccccACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTag                                                         <  2:448009/100‑1 (MQ=255)
                                             cccACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTaga                                                        >  2:1189920/1‑100 (MQ=255)
                                             cccACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTaga                                                        <  2:3814912/100‑1 (MQ=255)
                                              ccACGGGCATGGTGGATGCTGCGCACCGTTGGTGTAGAGAAAGGGTCTTTTC‑GGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGtt                                                           <  2:360228/96‑1 (MQ=255)
                                               cACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGG‑G‑GGG‑TGGACTTGCCGGCTGGCAGCGCGATGATCTGCTGTTagaag                                                      >  1:53074/1‑100 (MQ=255)
                                                  gAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGaagaag                                                   <  2:1583760/100‑1 (MQ=255)
                                                      aTGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGc                                               <  1:1488470/100‑1 (MQ=255)
                                                       tggtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCa                                              <  2:2711064/100‑1 (MQ=255)
                                                        ggtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAg                                             >  1:2249129/1‑100 (MQ=255)
                                                        ggtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAg                                             <  1:688865/100‑1 (MQ=255)
                                                        ggtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGGGGAGTCGG‑GG‑G‑GGG‑TGGACTTACAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAg                                             >  2:3747216/1‑100 (MQ=255)
                                                         gtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑T‑GG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTa                                           >  2:467031/1‑100 (MQ=255)
                                                           ggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTag                                          <  1:947640/100‑1 (MQ=255)
                                                            gATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTaga                                         <  2:527667/100‑1 (MQ=255)
                                                              tgctgcGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTagag                                        >  2:2789055/1‑99 (MQ=255)
                                                                ctgcGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTg                                     >  1:2925114/1‑100 (MQ=255)
                                                                ctgcGCACATTTGGTGTAGAGAAAGTGTCGATTC‑TGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTg                                     >  2:1552733/1‑100 (MQ=255)
                                                                cagcgcACCGTCGCTGTAGAGAAAGTGTAGATCC‑TGGG‑G‑GGG‑TGGACTTGTAGGCTGGCAGCGGGATGATCGGCGGTTAGAAGAAGGTGCAGGAGAGCTg                                     >  2:2614546/3‑100 (MQ=255)
                                                                 tgcGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGc                                    <  1:2527037/100‑1 (MQ=255)
                                                                   cgcACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCg                                  >  1:1673443/1‑100 (MQ=255)
                                                                   cgcACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCg                                  >  2:2423233/1‑100 (MQ=255)
                                                                     cACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGa                                >  1:1002708/1‑100 (MQ=255)
                                                                         tttGGTGTAGAGAAAGTGTCGATTC‑TGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAgg                             <  1:1152326/99‑1 (MQ=255)
                                                                         tttGGTGTAGAGAAAGTGTCGATTC‑TGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGga                            <  2:3655978/100‑1 (MQ=255)
                                                                           tGGTGTAGAGAAAGTGTCGATTC‑TGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGgaga                          <  1:603623/100‑1 (MQ=255)
                                                                             gtgtAGAGAAAGTGTCGATTC‑TGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGgagg                          >  2:3401435/1‑97 (MQ=255)
                                                                               gtAGAGAAAGTGTCGATTC‑TGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGGGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGt                        >  2:3663243/1‑98 (MQ=255)
                                                                               gtAGAGAAAGTGTCGATTC‑TGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGttt                      >  2:706318/1‑100 (MQ=255)
                                                                               gtAGAGAAAGTGTCGATTC‑TGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGttt                      <  2:3194921/100‑1 (MQ=255)
                                                                               gtAGAGAAAGTGTCGATTC‑TGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGttt                      <  2:79574/100‑1 (MQ=255)
                                                                                  gagaAAGTGTCGATTC‑TGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAc                   >  1:2963814/1‑100 (MQ=255)
                                                                                  gagaAAGTGTCGATTC‑TGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAc                   >  1:2005226/1‑100 (MQ=255)
                                                                                   agaAAGTGTCGATTC‑TGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAcg                  <  2:2196492/100‑1 (MQ=255)
                                                                                     aaaGTGTCGATTC‑TGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAcgcc                <  1:862651/100‑1 (MQ=255)
                                                                                       aGTGTCGATTC‑TGGG‑G‑GGG‑GGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAcgccgc              >  2:2595276/1‑100 (MQ=255)
                                                                                        gtgtCGATTC‑TGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTACCGCcgcg             <  1:1994365/100‑1 (MQ=255)
                                                                                        gtgtCGATTC‑TGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCcgcg             <  1:2277937/100‑1 (MQ=255)
                                                                                         tgtCGATTC‑TGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGt            >  1:523964/1‑100 (MQ=255)
                                                                                          gtCGGTTC‑TGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCcgcg             >  1:1596888/1‑98 (MQ=255)
                                                                                          gtCGATTC‑TGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGtt           >  2:1205411/1‑100 (MQ=255)
                                                                                          gtCGATTC‑TGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGtt           >  2:3332832/1‑100 (MQ=255)
                                                                                          gtCGATTC‑TGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGGAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGtt           >  2:2280283/1‑100 (MQ=255)
                                                                                          gtCGATTC‑TGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCCATGATCTGCTGTTAGGAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGCTTACcgccgc              >  2:2462283/1‑97 (MQ=255)
                                                                                             gATTC‑TGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTaa        >  1:1567308/1‑100 (MQ=255)
                                                                                             gATTC‑TGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTaa        <  2:1882740/100‑1 (MQ=255)
                                                                                             gATTC‑TGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTaa        <  1:3406829/100‑1 (MQ=255)
                                                                                              aTTC‑TGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAAt       >  2:2186340/1‑100 (MQ=255)
                                                                                                    |                                                                                                      
GGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑T‑GG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGA  >  NC_000913/2652751‑2652948

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑