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.G91.5% 411.8 / 31.8 129coding (356/846 nt)sseA3‑mercaptopyruvate sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base . (7/4);  new base G (47/70);  total (55/74)
Fisher's exact test for biased strand distribution p-value = 2.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 149 total aligned reads displayed.

GGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTC‑GATTC‑T‑G‑G‑GG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAAT  >  NC_000913/2652751‑2652943
                                                                                                     |                                                                                                  
ggATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTC‑GATTC‑T‑g                                                                                                   >  1:1552059/1‑100 (MQ=255)
   tAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCAGGGTGCATGCTGGGCGCCTTTGGTGTCGAGAAAGTGTC‑GATTC‑T‑g‑g                                                                                                 >  2:2745557/1‑98 (MQ=255)
    aaGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTC‑GATTC‑T‑g‑g‑gga                                                                                             <  1:1207111/100‑2 (MQ=255)
    aaGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTC‑GATTC‑T‑g‑g‑gg‑g                                                                                            <  2:3661036/100‑1 (MQ=255)
      gCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTC‑GATTC‑T‑g‑g‑ggagg                                                                                           >  1:3226861/1‑97 (MQ=255)
      gCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGATGTAGAGAAAGTGTC‑GATTC‑T‑g‑g‑gg‑ggg                                                                                          >  1:3429082/1‑99 (MQ=255)
       cACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTC‑GATTC‑T‑g‑g‑gg‑gggt                                                                                         <  2:2108155/100‑3 (MQ=255)
       cACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTC‑GATTC‑T‑g‑g‑gg‑gggt                                                                                         >  2:3340501/1‑98 (MQ=255)
        aCCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTC‑GATTC‑T‑g‑g‑ggagatc                                                                                         >  2:3735316/1‑95 (MQ=255)
         ccTGATTGTCTATGACGATGGTAATCTTTTCTAAGCCCCACGTGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTC‑GATTC‑T‑g‑g‑gg‑gggtgg                                                                                       <  1:574959/100‑5 (MQ=255)
         ccTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTC‑GATTC‑T‑g‑g‑gg‑gggtgg                                                                                       <  2:466954/100‑5 (MQ=255)
         ccTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTC‑GATTC‑T‑g‑g‑gg‑gggtgg                                                                                       >  2:908502/1‑96 (MQ=255)
         ccTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTC‑GATTC‑T‑g‑g‑gg‑gggtgg                                                                                       >  1:1178518/1‑96 (MQ=255)
         ccTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGAACACCTTTGGTGTAGAGAAAGTGTC‑GATTC‑T‑g‑g‑gg‑gggtgg                                                                                       >  1:267788/1‑96 (MQ=255)
          cTGATTGTCTATGACGAAGGTAATCCTTTCTCAGCCCCACGACCATGGTGGATGCTGCGCACCTTTGGTGTAGAGACAGTGTC‑GCTTC‑T‑g‑g‑gg‑gggtg                                                                                        >  2:601506/1‑95 (MQ=255)
           tGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTG‑GATTC‑T‑t‑g‑gg                                                                                              >  2:1658376/1‑92 (MQ=255)
           tGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTC‑TATTC‑TGG‑G‑GG‑GG‑TGGAc                                                                                     >  2:1143030/1‑100 (MQ=255)
           tGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTC‑GATTC‑TGG‑G‑GG‑GG‑TGGAc                                                                                     <  2:3481774/100‑1 (MQ=255)
            gATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTC‑GATTC‑TGG‑G‑GG‑GG‑TGGACt                                                                                    <  1:2204928/100‑1 (MQ=255)
            gATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTC‑GATTC‑TGG‑G‑GG‑GG‑TGGACt                                                                                    >  1:2646693/1‑100 (MQ=255)
            gATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTC‑GATTC‑TGG‑G‑GG‑GG‑TGGACt                                                                                    >  2:159931/1‑100 (MQ=255)
              ttGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTC‑GATTC‑TGG‑G‑GG‑GG‑TGGACTTg                                                                                  <  2:1976750/100‑1 (MQ=255)
                 tCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTC‑GATTC‑TGG‑G‑GG‑GG‑TGGACTTGCAg                                                                               <  1:2856610/100‑1 (MQ=255)
                  cTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTC‑GATTC‑TGG‑G‑GG‑GG‑TGGACTTGCAgg                                                                              >  2:2760431/1‑100 (MQ=255)
                  cTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGAAACTTTGATGTAGAGAAAGTGTC‑GATTC‑TGG‑G‑GG‑GG‑TGGACTTGCAgg                                                                              >  1:276991/1‑100 (MQ=255)
                   tATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGTTGGATGCTGCGCACCTTCGGTGTAGAGAAAGTGTC‑GATTC‑TGG‑G‑GG‑GG‑TGGACTTGCAGGc                                                                             <  1:3081979/100‑1 (MQ=255)
                   tATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTC‑GATTC‑TGG‑G‑GG‑GG‑TGGACTTGCAGGc                                                                             >  1:3808268/1‑100 (MQ=255)
                   tATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTC‑GATTC‑TGG‑G‑GG‑GG‑TGGACTTGCAGGc                                                                             <  2:2063488/100‑1 (MQ=255)
                   tATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTC‑GATTC‑TGG‑G‑GG‑GG‑TGGACTTGCAGGc                                                                             <  1:3411921/100‑1 (MQ=255)
                   tATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTC‑GATTC‑T‑G‑GTGG‑GG‑TGGACTTGCAGGc                                                                             <  2:616970/100‑1 (MQ=255)
                     tGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTC‑GATTC‑TGG‑G‑GG‑GG‑TGGACTTGCAGGCTg                                                                           <  2:3081658/100‑1 (MQ=255)
                     tGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTC‑GATTC‑TGG‑G‑GG‑GG‑TGGACTTGCAGGCTg                                                                           <  2:2582365/100‑1 (MQ=255)
                     tGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTC‑GATTC‑TGG‑G‑GG‑GG‑TGGACTTGCAGGCTg                                                                           <  2:2605601/100‑1 (MQ=255)
                      gACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTC‑GATTC‑TGG‑G‑GG‑GG‑TGGACTTGCAGGCTgg                                                                          <  2:2932904/100‑1 (MQ=255)
                      gACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTC‑GATTC‑TGG‑G‑GG‑GG‑TGGACTTGCAGGCTgg                                                                          <  1:140550/100‑1 (MQ=255)
                       aCGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTC‑GATTC‑TGG‑G‑GG‑GG‑TGGACTTGCAGGCTGGc                                                                         <  2:1057450/100‑1 (MQ=255)
                        cGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTC‑GATTC‑TGG‑G‑GG‑GG‑TGGACTTGCAGGCTGGCa                                                                        <  2:1822195/100‑1 (MQ=255)
                        cGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTC‑GATTA‑TGG‑G‑GG‑GG‑TGGACTTGCAGGCTGGCa                                                                        >  1:618169/1‑100 (MQ=255)
                           aGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTC‑GATTC‑TGG‑G‑GG‑GG‑TGGACTTGCAGGCTGGCAgcg                                                                     <  1:1884517/100‑1 (MQ=255)
                           aGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTC‑GATTC‑TGG‑G‑GG‑GG‑TGGACTTGCAGGCTGGCAgcg                                                                     <  1:1004557/100‑1 (MQ=255)
                            ggTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTC‑GATTC‑TGG‑G‑GG‑GG‑TGGACTTGCAGGCTGGCAgcgc                                                                    <  1:324447/100‑1 (MQ=255)
                            ggTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTC‑GATTC‑TGG‑G‑GG‑GG‑TGGACTTGCAGGCTGGCAgcgc                                                                    <  2:3625090/100‑1 (MQ=255)
                             gcaaTCTTTGCTCAGCCGCACGAGCATGGTGGATGCTGCGCACCTGTGGTGTAGAGAAAGTGTC‑GATTC‑TGG‑G‑GG‑GG‑TGGACTTGCAGGCTGGCAGcgcg                                                                   <  2:3728803/98‑1 (MQ=255)
                              tAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGACAAATAA‑G‑G‑G‑GG‑GG‑TGGACTTGCAGGCTGGCAGCGCga                                                                  >  1:1221832/1‑100 (MQ=255)
                              tAATCTTTTCTCAGCCACACGAGCATGGTGGATGCTGCGCACGTTTGGTGTAAAGAAAGTGTG‑GAGTC‑TGG‑G‑GG‑GG‑AGGACTTGCATGCTGGCAgcg                                                                     >  2:3413644/1‑97 (MQ=255)
                               aaTCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTC‑GATTC‑TGG‑G‑GG‑GG‑TGGACTTGCAGGCTGGCAGCGCgat                                                                 >  2:2375158/1‑100 (MQ=255)
                               aaTCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTC‑GATTC‑TGG‑G‑GG‑GG‑TGGACTTGCAGGCTGGCAGCGCgat                                                                 >  1:1093938/1‑100 (MQ=255)
                               aaTCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTC‑GATTC‑TGG‑G‑GG‑GG‑TGGACTTGCAGGCTGGCAGCGCgat                                                                 <  1:3653606/100‑1 (MQ=255)
                                aTCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTC‑GATTC‑TGG‑G‑GG‑GG‑TGGACTTGCAGGCTGGCAGCGCgatg                                                                <  2:688183/100‑1 (MQ=255)
                                aTCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCATTGGTGTAGAGAAAGTGTC‑GATTC‑TGG‑G‑GG‑GG‑TGGACTTGCAGGCTGGCAGCGCgatg                                                                <  2:2458365/100‑1 (MQ=255)
                                 tCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTC‑GATTC‑TGG‑G‑GG‑GG‑TGGACTTGCAGGCTGGCAGCGCgatga                                                               >  1:938887/1‑100 (MQ=255)
                                 tCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTC‑GATTC‑TGG‑G‑GG‑GG‑TGGACTTGCAGGCTGGCAGCGCgatga                                                               >  1:3875189/1‑100 (MQ=255)
                                     ttCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTC‑GATTC‑TGG‑G‑GG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATctg                                                           <  1:3780389/100‑1 (MQ=255)
                                      tctcAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTC‑GATTC‑TGG‑G‑GG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATctgc                                                          <  1:436984/100‑1 (MQ=255)
                                       ctcAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTCTAGAGAAAGTGTC‑GAATA‑TAA‑G‑GG‑AA‑TGGACTAGCAGGCTGGCAGCGCGATGATctgct                                                         >  1:444739/1‑100 (MQ=255)
                                          aGCCCCACTAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTC‑GATTC‑TGG‑G‑GG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGtt                                                      <  1:656347/100‑1 (MQ=255)
                                          aGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTC‑GATTC‑TGG‑G‑GG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGtt                                                      <  1:159936/100‑1 (MQ=255)
                                          aGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTC‑GATTC‑TGG‑G‑GG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGtt                                                      <  1:3007484/100‑1 (MQ=255)
                                             cccACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTC‑GATTC‑TGG‑G‑GG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTaga                                                   <  1:3814971/100‑1 (MQ=255)
                                             cccACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTC‑GATTC‑TGG‑G‑GG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTaga                                                   >  2:2494550/1‑100 (MQ=255)
                                              ccACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTC‑GATTC‑TGG‑G‑GG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTagaa                                                  <  1:812889/100‑1 (MQ=255)
                                              ccACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTC‑GATTC‑TGG‑G‑GG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTagaa                                                  <  2:3359779/100‑1 (MQ=255)
                                              ccACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTC‑GATTC‑TGG‑G‑GG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTaga                                                   <  2:1957792/99‑1 (MQ=255)
                                                 cGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTC‑GATTC‑TGG‑G‑GG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAgaagaa                                               >  2:121030/1‑100 (MQ=255)
                                                 cGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTC‑GATTC‑TGG‑G‑GG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAgaagaa                                               >  1:1544722/1‑100 (MQ=255)
                                                   aGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTC‑GATTC‑TGG‑G‑GG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAgg                                             <  2:396352/100‑1 (MQ=255)
                                                    gCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTC‑GATTC‑TGG‑G‑GG‑GG‑TGGACTTGCAGGCAGACAGCACGAAGAACTGCTATTAGAAGAAGGt                                            >  1:1444297/1‑100 (MQ=255)
                                                       tggtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTC‑GATTC‑T‑G‑G‑GGAGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCa                                         >  2:2034310/1‑100 (MQ=255)
                                                        ggtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTC‑GATTC‑TGG‑G‑GG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAg                                        >  1:1487043/1‑100 (MQ=255)
                                                         gtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTC‑GATTC‑TGG‑G‑GG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGt                                       <  1:1739498/100‑1 (MQ=255)
                                                          tggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTC‑GATTC‑TGG‑G‑GG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTa                                      <  1:2001043/100‑1 (MQ=255)
                                                           ggATTCTGCGCACCTTTTGTGTAGTGAAAGCGTC‑GATTC‑TGG‑G‑GG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTag                                     <  1:243766/100‑1 (MQ=255)
                                                            gATGCTGCGCACCTTTGGTGTAGAGAAAGTGTC‑GATTC‑T‑G‑GAGG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTaga                                    <  2:233417/100‑1 (MQ=255)
                                                              tgctgcGCACCTTTGGTGTAGAGAAATATTC‑GTTTCTT‑T‑G‑GG‑GC‑TGGACGTGCAGGTTGGCAGCGTGATGATCTTCTTTTATAAGATGGTGCAGTAGAGc                                  <  2:293067/100‑1 (MQ=255)
                                                              tgctgcGCACCTTTGGTGTAGAGAAAGTGTC‑GATTC‑TGG‑G‑GG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGc                                  >  1:3044588/1‑100 (MQ=255)
                                                                ctgcGCACCTTTGGTGTAGAGAAAGAGTC‑GATTC‑TGG‑G‑GG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCAGATAGAAGAAGGAGCAGTaga                                    >  2:2861242/1‑96 (MQ=255)
                                                                 tgcGCACCTTTGGTGTAGAGAAAGTGTC‑GATTC‑TGG‑G‑GG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGc                               >  2:942459/1‑100 (MQ=255)
                                                                  gcgtACCTTTTTTGTAGAGATAGTGTC‑GATTC‑T‑GTT‑TG‑GG‑GAGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGcc                              <  1:2276162/100‑1 (MQ=255)
                                                                    gcACCTTTGGTGTAGAGAAAGTGTC‑GATTC‑TGG‑G‑GG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCgg                            <  2:3449299/100‑1 (MQ=255)
                                                                    gcACCTTTGGTGTAGAGAAAGTGTC‑GATTC‑TGG‑G‑GG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCgg                            >  1:1945783/1‑100 (MQ=255)
                                                                     cACCTTTGGTGTAGAGAAAGTGTC‑GATTC‑TGG‑G‑GG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGa                           <  1:1189984/100‑1 (MQ=255)
                                                                        cTTTGGTGTAGAGAAAGTGTC‑GATTC‑TGG‑G‑GG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAgg                        >  1:2490258/1‑100 (MQ=255)
                                                                            ggTGTAGAGAAAGTGTC‑GATTC‑TGG‑G‑GG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGagag                    >  1:3170725/1‑100 (MQ=255)
                                                                            ggTGTAGAGAAAGTGGC‑GATTC‑TGG‑G‑GG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGagag                    <  2:3022646/100‑1 (MQ=255)
                                                                             gtgtAGAGAAAGTGTC‑GATTC‑TGG‑G‑GG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGt                   >  1:562917/1‑100 (MQ=255)
                                                                              tgtAGAGAAAGTGTC‑GATTC‑TGG‑G‑GG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGtt                  <  2:749887/100‑1 (MQ=255)
                                                                               gtAGAGAAAGTGTC‑GATTC‑TGG‑G‑GG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGttt                 <  1:3685817/100‑1 (MQ=255)
                                                                                 agagAAAGTGTC‑GATTC‑TGG‑G‑GG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTaa               >  2:1999300/1‑100 (MQ=255)
                                                                                 agagAAAGTGTC‑GATTC‑TGG‑G‑GG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGGAGAGCTGCCGGAAGGACACTTTaa               >  2:3017162/1‑100 (MQ=255)
                                                                                 agagAAAGTGTC‑GATTC‑TGG‑G‑GG‑GG‑TGGAATTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCCGAAGgaga                     >  2:3110966/1‑94 (MQ=255)
                                                                                 agagAAAGTGTC‑GATTC‑T‑G‑G‑GG‑GGATGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGGAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTaa               >  2:1149899/1‑100 (MQ=255)
                                                                                  gagaAAGTGTC‑GATTC‑TGG‑G‑GG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAc              >  1:112348/1‑100 (MQ=255)
                                                                                  gagaAAGTGTC‑GATTC‑TGG‑G‑GG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAc              <  2:1301369/100‑1 (MQ=255)
                                                                                  gagaAAGTGTC‑GATTC‑TGG‑G‑GG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAc              >  2:3338684/1‑100 (MQ=255)
                                                                                   agaAAGTGTC‑GATTC‑TGG‑G‑GG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAcg             <  2:3836821/100‑1 (MQ=255)
                                                                                     aaaGTGTC‑GAGTC‑TGG‑G‑GG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAcgc            <  1:957479/99‑1 (MQ=255)
                                                                                       aGTGTC‑GATTC‑T‑G‑GTGG‑GG‑TGGACTTGCAGGCTGGCAGAGCGATGATCTGCTGTTAGAAGAAGTTGCAGTAGAGCTGCCGGAAGGAGAGGTTAAcgc            >  2:2072469/1‑97 (MQ=255)
                                                                                          gtC‑GATTC‑TGG‑G‑GG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGtt      >  2:1330288/1‑100 (MQ=255)
                                                                                           tC‑GATTC‑TGG‑G‑GG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGttt     <  2:2914093/100‑1 (MQ=255)
                                                                                            c‑GATTC‑TGG‑G‑GG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTa    <  2:430565/100‑1 (MQ=255)
                                                                                                     |                                                                                                  
GGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTC‑GATTC‑T‑G‑G‑GG‑GG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAAT  >  NC_000913/2652751‑2652943

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