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.G95.7% 405.3 / 11.1 115coding (356/846 nt)sseA3‑mercaptopyruvate sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base . (4/1);  new base G (51/59);  total (55/60)
Fisher's exact test for biased strand distribution p-value = 1.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.92e-01

Only 100 of 123 total aligned reads displayed.

AGGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CT‑‑GG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAAT  >  NC_000913/2652750‑2652943
                                                                                                     |                                                                                                
aGGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑Ct                                                                                                   >  1:1082716/1‑100 (MQ=255)
    tAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CT‑‑gg‑gg                                                                                            <  2:2349878/100‑1 (MQ=255)
      aGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CT‑‑gg‑gggg                                                                                          >  2:1574803/1‑100 (MQ=255)
       gCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCACACCTTTGGTGTAGAGAAAGTGTCGATT‑CT‑‑gg‑ggggg                                                                                         <  2:271568/100‑2 (MQ=255)
          ccTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CT‑‑gg‑gggggtgg                                                                                      <  1:2956305/100‑5 (MQ=255)
          ccTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CT‑‑gg‑gggggtgg                                                                                      >  2:319875/1‑96 (MQ=255)
           cTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CT‑‑gg‑gggggtgga                                                                                     >  2:748585/1‑95 (MQ=255)
            tGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTG‑GG‑GGGGTGGAc                                                                                     >  2:213927/1‑100 (MQ=255)
             gATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTG‑GG‑GGGGTGGACt                                                                                    <  2:800216/100‑1 (MQ=255)
             gATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTG‑GG‑GGGGTGGACt                                                                                    >  2:2199306/1‑100 (MQ=255)
             gATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGCAGAGAAAGTGTCGATT‑CT‑‑gg‑gggggtggatt                                                                                   <  2:2062614/100‑8 (MQ=255)
             gATTGTATATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTG‑GG‑GGGGTGGACt                                                                                    <  2:1409822/100‑1 (MQ=255)
              aTTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTG‑GG‑GGGGTGGAcct                                                                                   >  1:1854519/1‑98 (MQ=255)
               ttGTCTACGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTG‑GG‑GGGGTGGACTTg                                                                                  >  2:1606768/1‑100 (MQ=255)
                 gTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATTCTGCGCAGCTTTGGTGTAGAGAAAGTGTCGATT‑CTG‑GG‑GGGGTGGACTTGc                                                                                 <  1:676120/99‑1 (MQ=255)
                 gTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTG‑GG‑GGGGTGGACTTGc                                                                                 <  2:2513171/99‑1 (MQ=255)
                  tCTATGACGAAGGTAATCTTTTCTCAGCCCCACGATCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTG‑GG‑GGGGTGGACTTGc                                                                                 >  1:2358151/1‑98 (MQ=255)
                  tCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTG‑GG‑GGGGTGGACTTGCAg                                                                               <  1:2972418/100‑1 (MQ=255)
                  tCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTG‑GG‑GGGGTGGACTTGCAg                                                                               <  1:2784214/100‑1 (MQ=255)
                   cTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTCGAGCAAGTGTCGATT‑CTG‑GG‑GGGGTGGACTTGCAgg                                                                              >  1:936167/1‑100 (MQ=255)
                    tATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTG‑GG‑GGGGTGGACTTGCAGGc                                                                             <  2:1115483/100‑1 (MQ=255)
                    tATGACGAAGGTAATCTTTTCTCAGCCCCACGACCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CT‑‑GG‑GGGGTGGACTTGCAGGCt                                                                            >  1:2003650/1‑100 (MQ=255)
                      tGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTG‑GG‑GGGGTGGACTTGCAGGCTg                                                                           <  1:2199293/100‑1 (MQ=255)
                      tGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTG‑GG‑GGGGTGGACTTGCAGGCTg                                                                           <  2:1376617/100‑1 (MQ=255)
                      tGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTG‑GG‑GGGGTGGACTTGCAGGCTg                                                                           <  1:610614/100‑1 (MQ=255)
                       gACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTG‑GG‑GGGGTGGACTTGCAGGCTgg                                                                          <  2:1502523/100‑1 (MQ=255)
                        aCGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTG‑GG‑GGGGTGGACTTGCAGGCTGGc                                                                         <  1:821458/100‑1 (MQ=255)
                         cGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTGGGG‑GGGGTGGACTTGCAGGCTGGc                                                                         <  1:94684/100‑1 (MQ=255)
                         cGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTG‑GG‑GGGGTGGACTTGCAGGCTGGCa                                                                        >  2:2476203/1‑100 (MQ=255)
                         cGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTG‑GG‑GGGGTGGACTTGCAGGCTGGCa                                                                        >  2:144320/1‑100 (MQ=255)
                               tAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTG‑GG‑GGGGTGGACTTGCAGGCTGGCAGCGCga                                                                  >  2:1072784/1‑100 (MQ=255)
                               tAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTG‑GG‑GGGGTGGACTTGCAGGCTGGCAGCGCga                                                                  >  2:1727608/1‑100 (MQ=255)
                                  tCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTG‑GG‑GGGGTGGACTTGCAGGCTGGCAGCGCgatga                                                               <  2:1082722/100‑1 (MQ=255)
                                   cTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTG‑GG‑GGGGCGGACTTGCAGGCTGGCAGCGCgatga                                                               >  1:472747/1‑99 (MQ=255)
                                    ttttCTCAGCGCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTTGATT‑CTG‑GG‑GGGGTGGACTTGCAGGATGGCAGCGCGATGATc                                                             <  1:2846622/100‑1 (MQ=255)
                                    ttttCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTG‑GG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATc                                                             >  1:2696271/1‑100 (MQ=255)
                                    ttttCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTG‑GG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATc                                                             >  2:3027957/1‑100 (MQ=255)
                                    ttttCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTG‑GG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATc                                                             >  2:783914/1‑100 (MQ=255)
                                    ttttCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTG‑GG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATc                                                             <  1:1796609/100‑1 (MQ=255)
                                    ttttCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTG‑GG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATc                                                             <  1:1991125/100‑1 (MQ=255)
                                    ttttCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTG‑GG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATc                                                             <  1:2021397/100‑1 (MQ=255)
                                     tttCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTG‑GG‑GGGGGGGACTTGCAGGCTGGCAGCGCGATGATct                                                            >  2:1302631/1‑100 (MQ=255)
                                      ttCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTG‑GG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATctg                                                           >  1:43521/1‑100 (MQ=255)
                                       tctcAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTG‑GG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATctgc                                                          <  2:1396647/100‑1 (MQ=255)
                                       tctcAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTG‑GG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATctgc                                                          <  2:159605/100‑1 (MQ=255)
                                        ctcAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTG‑GG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATctgct                                                         >  2:2712771/1‑100 (MQ=255)
                                         tcAGCCCCACGGGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTG‑GG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATctgctg                                                        <  2:496013/100‑1 (MQ=255)
                                         tcAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTG‑GG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATctgctg                                                        >  1:148236/1‑100 (MQ=255)
                                         tcAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTG‑GG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATctgctg                                                        <  1:2264605/100‑1 (MQ=255)
                                          cAGCCCCACGAGCATGGTGGATGCTGCGCGCCTTTGGTGTAGAGAAAGTGTCGATT‑CTG‑GG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATctgctg                                                        <  1:1004903/99‑1 (MQ=255)
                                           aGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTG‑GG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGtt                                                      >  2:125357/1‑100 (MQ=255)
                                           aGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTG‑GG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGtt                                                      <  1:2944641/100‑1 (MQ=255)
                                           aGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTG‑GG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGtt                                                      <  2:1549819/100‑1 (MQ=255)
                                            gCCCCACGAGCATGGTGGATGCTGCGCACCGTTGGTGTAGAGAAAGTGTCGATT‑CTG‑GG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTa                                                     <  1:2654064/100‑1 (MQ=255)
                                             ccccACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTG‑GG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTag                                                    <  2:1007764/100‑1 (MQ=255)
                                               ccACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTG‑GG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTagaa                                                  <  2:1854530/100‑1 (MQ=255)
                                                 aCGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTG‑GG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTagaaga                                                <  1:1399209/100‑1 (MQ=255)
                                                  cGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTG‑GG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAgaagaa                                               <  1:1431327/100‑1 (MQ=255)
                                                  cGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTG‑GG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAgaagaa                                               >  1:1746342/1‑100 (MQ=255)
                                                   gAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CT‑‑GG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAgg                                             <  2:1990167/100‑1 (MQ=255)
                                                    aGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTG‑GG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAgg                                             >  1:703127/1‑100 (MQ=255)
                                                    aGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTG‑GG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAgg                                             <  2:2823788/100‑1 (MQ=255)
                                                    aGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTG‑GG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAgg                                             <  2:1408835/100‑1 (MQ=255)
                                                       aTGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTG‑GG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGc                                          <  1:463912/100‑1 (MQ=255)
                                                         ggtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTG‑GG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAg                                        <  2:1587939/100‑1 (MQ=255)
                                                         ggtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTG‑GG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAg                                        <  1:1075616/100‑1 (MQ=255)
                                                         ggtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTG‑GG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAg                                        <  2:471739/100‑1 (MQ=255)
                                                          gtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTG‑GG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGt                                       <  2:2884906/100‑1 (MQ=255)
                                                          gtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTG‑GG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGt                                       <  2:2883635/100‑1 (MQ=255)
                                                           tggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTG‑GG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTa                                      <  2:426530/100‑1 (MQ=255)
                                                             gATGCTGCGCACCGTTGGTGTAGAGAAAGTGACGATT‑CTG‑GG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATGTGCTGTTAGAAGAAGGTGCAGTaga                                    >  1:1520814/1‑100 (MQ=255)
                                                               tgctgcGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTG‑GG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGc                                  >  2:3060345/1‑100 (MQ=255)
                                                               tgctgcGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTG‑GG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGc                                  >  2:3060261/1‑100 (MQ=255)
                                                               tgctgcGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTG‑GG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGc                                  <  2:43522/100‑1 (MQ=255)
                                                               tgctgcGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTG‑GG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGc                                  >  2:1942825/1‑100 (MQ=255)
                                                               tgctgcGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTG‑GG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGc                                  <  2:1454686/100‑1 (MQ=255)
                                                                gctgcGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTG‑GG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCt                                 >  2:2750043/1‑100 (MQ=255)
                                                                 ctgcGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CTG‑GG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTg                                <  2:2936928/100‑1 (MQ=255)
                                                                   gcgcACCTTTGGTGTAGAGAAAGTGTCGATT‑CTG‑GG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGcc                              <  1:367248/100‑1 (MQ=255)
                                                                   gcgcACCTTTGGTGTAGAGAAAGTGTCGATT‑CTG‑GG‑GGGGTGGACGTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGcc                              <  2:1237644/100‑1 (MQ=255)
                                                                       aCCTTTGGTGTAGAGAAAGTGTCGATT‑CTG‑GG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGa                           <  1:597661/99‑1 (MQ=255)
                                                                         cTTTGGTGTAGAGAAAGTGTCGATT‑CTG‑GG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAgg                        <  1:2298804/100‑1 (MQ=255)
                                                                          tttGGTGTAGAGAAAGTGTCGATT‑CTGGGG‑GGGGTGTACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAgg                        <  2:1855935/100‑1 (MQ=255)
                                                                               tgtAGAGAAAGTGTCGATT‑CTG‑GG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGtt                  >  1:1177560/1‑100 (MQ=255)
                                                                                 ccgatctaGTGTCGATT‑CTG‑GG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTa                <  1:558254/93‑1 (MQ=255)
                                                                                  agagAAAGTGTCGATT‑CTG‑GG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTaa               >  1:62709/1‑100 (MQ=255)
                                                                                    agaAAGTGTCGATT‑CTG‑GG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAcg             >  2:149543/1‑100 (MQ=255)
                                                                                     gaAAGTGTCGATTCCG‑‑GG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAcgc            >  1:2472716/1‑100 (MQ=255)
                                                                                     gaAAGTGTCGATT‑CTG‑GG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAcgc            >  1:1660151/1‑100 (MQ=255)
                                                                                      aaaGTGTCGATT‑CTG‑GG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAcgcc           <  2:936172/100‑1 (MQ=255)
                                                                                       aaGTGTCGATT‑CTG‑GG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAcgccg          >  2:2716636/1‑100 (MQ=255)
                                                                                         gagtCGATT‑CT‑‑GGAGGGGTGGACTTGCAGGCAGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCcgcg        >  2:2057578/3‑100 (MQ=255)
                                                                                           gtCGATT‑CTG‑GG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGtt      >  2:3025604/1‑100 (MQ=255)
                                                                                           gtCGATT‑CTG‑GG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGtt      >  1:746839/1‑100 (MQ=255)
                                                                                           gtCGATT‑CTG‑GG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGtt      <  1:1003339/100‑1 (MQ=255)
                                                                                             cGATT‑CTG‑GG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGttt     >  1:1053945/1‑99 (MQ=255)
                                                                                             cGATT‑CTG‑GG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTa    >  1:734393/1‑100 (MQ=255)
                                                                                             cGATT‑CTG‑GG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTa    >  1:1603378/1‑100 (MQ=255)
                                                                                              gATT‑CTG‑GG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTaa   >  1:581002/1‑100 (MQ=255)
                                                                                               aTT‑CTG‑GG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAAt  >  2:1852101/1‑100 (MQ=255)
                                                                                                     |                                                                                                
AGGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATT‑CT‑‑GG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAAT  >  NC_000913/2652750‑2652943

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