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.G96.6% 419.6 / 6.5 118coding (356/846 nt)sseA3‑mercaptopyruvate sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base . (4/0);  new base G (57/57);  total (61/57)
Fisher's exact test for biased strand distribution p-value = 1.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 135 total aligned reads displayed.

GATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑GGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAAT  >  NC_000913/2652752‑2652943
                                                                                                  |                                                                                                
 aTAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑ggg                                                                                               <  2:3119765/100‑1 (MQ=255)
  tAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑gggg                                                                                              <  1:1565173/100‑1 (MQ=255)
    aGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑gggg‑gg                                                                                           >  2:842206/1‑100 (MQ=255)
     gCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑gggg‑ggg                                                                                          >  2:3284926/1‑99 (MQ=255)
     gCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑gggg‑ggg                                                                                          <  2:2099290/100‑2 (MQ=255)
     gCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑gggg‑ggg                                                                                          >  2:1804217/1‑99 (MQ=255)
     gCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑gggg‑ggg                                                                                          >  2:2578661/1‑99 (MQ=255)
     gCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑gggg                                                                                              >  1:2642264/1‑97 (MQ=255)
        ccTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑gggg‑gggtgg                                                                                       >  1:1434935/1‑96 (MQ=255)
        ccTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑gggg‑gggtgg                                                                                       <  1:3410679/100‑5 (MQ=255)
        ccTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑gggg‑gggtgg                                                                                       <  2:2269073/100‑5 (MQ=255)
          tGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGAc                                                                                     >  2:2561903/1‑100 (MQ=255)
          tGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGAc                                                                                     >  2:3345268/1‑100 (MQ=255)
          tGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGAc                                                                                     <  2:2735824/100‑1 (MQ=255)
          tGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGAc                                                                                     >  2:1871621/1‑100 (MQ=255)
          tGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGAc                                                                                     >  1:1158521/1‑100 (MQ=255)
           gATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACt                                                                                    >  2:2898606/1‑100 (MQ=255)
           gATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACt                                                                                    <  1:2863456/100‑1 (MQ=255)
            aTTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACtt                                                                                   >  1:3372704/1‑100 (MQ=255)
             ttGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTg                                                                                  <  1:1804201/100‑1 (MQ=255)
             ttGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTg                                                                                  <  2:2067312/100‑1 (MQ=255)
               gTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCa                                                                                <  2:2662073/100‑1 (MQ=255)
               gTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCa                                                                                >  1:1603430/1‑100 (MQ=255)
                tCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAg                                                                               <  1:2070073/100‑1 (MQ=255)
                tCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAg                                                                               >  2:1774657/1‑100 (MQ=255)
                  tATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGc                                                                             <  1:1134216/100‑1 (MQ=255)
                  tATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGc                                                                             >  2:1292276/1‑100 (MQ=255)
                  tATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTGGCAGGc                                                                             >  2:1369411/1‑100 (MQ=255)
                   aTGCCGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGCAGAGAAAGTGGCGATTCTGGGGG‑GG‑TGGGCGTGCAGGCt                                                                            >  1:1498670/1‑100 (MQ=255)
                    tGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTg                                                                           <  1:467643/100‑1 (MQ=255)
                     gACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTgg                                                                          <  1:289741/100‑1 (MQ=255)
                       cGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCa                                                                        >  2:462191/1‑100 (MQ=255)
                             tAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCga                                                                  >  2:829663/1‑100 (MQ=255)
                             tAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCga                                                                  >  2:2104359/1‑100 (MQ=255)
                              aaTCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑GGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCgat                                                                 >  1:1645938/1‑99 (MQ=255)
                               aTCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCgatg                                                                <  1:412268/100‑1 (MQ=255)
                                tCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCgatga                                                               >  1:2766060/1‑100 (MQ=255)
                                tCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCgatga                                                               >  2:2448636/1‑100 (MQ=255)
                                tCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCgatga                                                               >  1:2763240/1‑100 (MQ=255)
                                 cTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCgatgat                                                              >  2:1742281/1‑100 (MQ=255)
                                 cTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCgatgat                                                              >  1:3262087/1‑100 (MQ=255)
                                       tcAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATctgctg                                                        <  2:917217/100‑1 (MQ=255)
                                        cAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGt                                                       >  1:1989915/1‑100 (MQ=255)
                                         aGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGtt                                                      <  1:344556/100‑1 (MQ=255)
                                         aGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGtt                                                      >  2:2836879/1‑100 (MQ=255)
                                         aGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGtt                                                      >  2:2737185/1‑100 (MQ=255)
                                          gCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTa                                                     <  2:811861/100‑1 (MQ=255)
                                           ccccACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGCGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCACCGCGATGATCTGCTGTTa                                                     <  1:249425/99‑1 (MQ=255)
                                            cccACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTaga                                                   <  1:2898583/100‑1 (MQ=255)
                                                cGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTagaagta                                               <  1:3329690/100‑3 (MQ=255)
                                                  aGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAgg                                             >  2:377506/1‑100 (MQ=255)
                                                  aGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAgg                                             <  2:1509284/100‑1 (MQ=255)
                                                      tggtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCa                                         >  1:1834801/1‑100 (MQ=255)
                                                       ggtggATGCTGCGCCCCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGGTGTTAGACGAAGGTGCAg                                        >  1:1011495/1‑100 (MQ=255)
                                                       ggtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAg                                        <  1:2317640/100‑1 (MQ=255)
                                                       ggtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAg                                        >  2:1745285/1‑100 (MQ=255)
                                                        gtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGt                                       <  1:1774641/100‑1 (MQ=255)
                                                        gtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGt                                       <  2:1516627/100‑1 (MQ=255)
                                                        gtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGt                                       <  1:228607/100‑1 (MQ=255)
                                                        gtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGt                                       >  2:1777973/1‑100 (MQ=255)
                                                        gtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGt                                       >  2:3082060/1‑100 (MQ=255)
                                                         tggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTa                                      <  2:234176/100‑1 (MQ=255)
                                                          ggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTag                                     <  1:2578642/100‑1 (MQ=255)
                                                          ggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTag                                     >  2:1603606/1‑100 (MQ=255)
                                                          ggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTag                                     <  2:2400844/100‑1 (MQ=255)
                                                            aTGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTagag                                   <  1:2263046/100‑1 (MQ=255)
                                                               ctgcGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTg                                >  2:935533/1‑100 (MQ=255)
                                                                tgcGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGc                               >  1:1422324/1‑100 (MQ=255)
                                                                tgcGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGc                               >  2:3395222/1‑100 (MQ=255)
                                                                   gcACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGTTTATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCgg                            >  1:1420912/1‑100 (MQ=255)
                                                                   gcACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCgg                            <  1:2517194/100‑1 (MQ=255)
                                                                    cACCTTTGGTGTAGAGAAAGTGTCGATTTTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGa                           <  2:1210409/100‑1 (MQ=255)
                                                                    cACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGTCGGa                           <  1:97109/100‑1 (MQ=255)
                                                                    cACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGa                           <  2:1227975/100‑1 (MQ=255)
                                                                    cACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGa                           <  1:1799189/100‑1 (MQ=255)
                                                                    cACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGa                           <  2:364933/100‑1 (MQ=255)
                                                                    cACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGa                           <  1:364644/100‑1 (MQ=255)
                                                                    cACCTTTGGTGTAGAGAAAGGGTCGAGTCT‑GGGGAGT‑TGGACTTGAAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGa                           >  1:1656850/1‑100 (MQ=255)
                                                                      ccTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAg                         >  2:157302/1‑100 (MQ=255)
                                                                      ccTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAg                         <  2:2410615/100‑1 (MQ=255)
                                                                       cTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAgg                        <  2:1886146/100‑1 (MQ=255)
                                                                       cTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAgg                        <  2:569785/100‑1 (MQ=255)
                                                                       cTTTGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAgg                        <  2:1478427/100‑1 (MQ=255)
                                                                         ttGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTACCGGAAGgag                      <  1:462187/100‑1 (MQ=255)
                                                                          tGGTGTAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGgaga                     >  2:2815280/1‑100 (MQ=255)
                                                                              gtAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGttt                 <  1:3284902/100‑1 (MQ=255)
                                                                              gtAGAGAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGttt                 <  1:1742265/100‑1 (MQ=255)
                                                                                agagAAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTaa               >  1:987284/1‑100 (MQ=255)
                                                                                 gagaAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAc              <  1:2543021/100‑1 (MQ=255)
                                                                                  agaAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGCAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAcg             >  1:462909/1‑100 (MQ=255)
                                                                                  agaAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAcg             <  1:2202712/100‑1 (MQ=255)
                                                                                  agaAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAcg             >  2:581591/1‑100 (MQ=255)
                                                                                   gaAAGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAcgc            >  2:302735/1‑100 (MQ=255)
                                                                                      aGTGTCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAcgccgc         <  1:2561884/100‑1 (MQ=255)
                                                                                         gtCGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGtt      >  2:1909140/1‑100 (MQ=255)
                                                                                           cGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGATCTGCCGGAAGGAGAGTTTAACGCCGCGTTTa    >  2:83855/1‑100 (MQ=255)
                                                                                           cGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTa    >  1:88921/1‑100 (MQ=255)
                                                                                           cGATTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTa    >  1:2568120/1‑100 (MQ=255)
                                                                                           cGATTCT‑GGGG‑GGTTGGACTTGCAGGCTGGCCGCGCGATGATGTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTa    >  2:2336355/1‑100 (MQ=255)
                                                                                             aTTCTGGGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAAt  <  1:196824/100‑1 (MQ=255)
                                                                                                  |                                                                                                
GATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑GGGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAAT  >  NC_000913/2652752‑2652943

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