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.

GGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAAT  >  NC_000913/2652751‑2652943
                                                                                                     |                                                                                                  
ggATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTg                                                                                                   >  1:1552059/1‑100 (MQ=255)
   tAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCAGGGTGCATGCTGGGCGCCTTTGGTGTCGAGAAAGTGTCGATTCTgg                                                                                                 >  2:2745557/1‑98 (MQ=255)
    aaGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTgggga                                                                                             <  1:1207111/100‑2 (MQ=255)
    aaGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTggggg                                                                                            <  2:3661036/100‑1 (MQ=255)
      gCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTggggagg                                                                                           >  1:3226861/1‑97 (MQ=255)
      gCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGATGTAGAGAAAGTGTCGATTCTggggggg                                                                                          >  1:3429082/1‑99 (MQ=255)
       cACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTgggggggt                                                                                         <  2:2108155/100‑3 (MQ=255)
       cACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTgggggggt                                                                                         >  2:3340501/1‑98 (MQ=255)
        aCCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTggggagatc                                                                                         >  2:3735316/1‑95 (MQ=255)
         ccTGATTGTCTATGACGATGGTAATCTTTTCTAAGCCCCACGTGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTgggggggtgg                                                                                       <  1:574959/100‑5 (MQ=255)
         ccTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTgggggggtgg                                                                                       <  2:466954/100‑5 (MQ=255)
         ccTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTgggggggtgg                                                                                       >  2:908502/1‑96 (MQ=255)
         ccTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTgggggggtgg                                                                                       >  1:1178518/1‑96 (MQ=255)
         ccTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGAACACCTTTGGTGTAGAGAAAGTGTCGATTCTgggggggtgg                                                                                       >  1:267788/1‑96 (MQ=255)
          cTGATTGTCTATGACGAAGGTAATCCTTTCTCAGCCCCACGACCATGGTGGATGCTGCGCACCTTTGGTGTAGAGACAGTGTCGCTTCTgggggggtg                                                                                        >  2:601506/1‑95 (MQ=255)
           tGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTGGATTCTtggg                                                                                              >  2:1658376/1‑92 (MQ=255)
           tGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCTATTCTGGGGGGGTGGAc                                                                                     >  2:1143030/1‑100 (MQ=255)
           tGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGAc                                                                                     <  2:3481774/100‑1 (MQ=255)
            gATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACt                                                                                    <  1:2204928/100‑1 (MQ=255)
            gATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACt                                                                                    >  1:2646693/1‑100 (MQ=255)
            gATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACt                                                                                    >  2:159931/1‑100 (MQ=255)
              ttGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTg                                                                                  <  2:1976750/100‑1 (MQ=255)
                 tCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAg                                                                               <  1:2856610/100‑1 (MQ=255)
                  cTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAgg                                                                              >  2:2760431/1‑100 (MQ=255)
                  cTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGAAACTTTGATGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAgg                                                                              >  1:276991/1‑100 (MQ=255)
                   tATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGTTGGATGCTGCGCACCTTCGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGc                                                                             <  1:3081979/100‑1 (MQ=255)
                   tATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGc                                                                             >  1:3808268/1‑100 (MQ=255)
                   tATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGc                                                                             <  2:2063488/100‑1 (MQ=255)
                   tATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGc                                                                             <  1:3411921/100‑1 (MQ=255)
                   tATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGTGGGGTGGACTTGCAGGc                                                                             <  2:616970/100‑1 (MQ=255)
                     tGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTg                                                                           <  2:3081658/100‑1 (MQ=255)
                     tGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTg                                                                           <  2:2582365/100‑1 (MQ=255)
                     tGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTg                                                                           <  2:2605601/100‑1 (MQ=255)
                      gACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTgg                                                                          <  2:2932904/100‑1 (MQ=255)
                      gACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTgg                                                                          <  1:140550/100‑1 (MQ=255)
                       aCGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGc                                                                         <  2:1057450/100‑1 (MQ=255)
                        cGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCa                                                                        <  2:1822195/100‑1 (MQ=255)
                        cGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTATGGGGGGGTGGACTTGCAGGCTGGCa                                                                        >  1:618169/1‑100 (MQ=255)
                           aGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAgcg                                                                     <  1:1884517/100‑1 (MQ=255)
                           aGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAgcg                                                                     <  1:1004557/100‑1 (MQ=255)
                            ggTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAgcgc                                                                    <  1:324447/100‑1 (MQ=255)
                            ggTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAgcgc                                                                    <  2:3625090/100‑1 (MQ=255)
                             gcaaTCTTTGCTCAGCCGCACGAGCATGGTGGATGCTGCGCACCTGTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGcgcg                                                                   <  2:3728803/98‑1 (MQ=255)
                              tAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGACAAATAAGGGGGGGTGGACTTGCAGGCTGGCAGCGCga                                                                  >  1:1221832/1‑100 (MQ=255)
                              tAATCTTTTCTCAGCCACACGAGCATGGTGGATGCTGCGCACGTTTGGTGTAAAGAAAGTGTGGAGTCTGGGGGGGAGGACTTGCATGCTGGCAgcg                                                                     >  2:3413644/1‑97 (MQ=255)
                               aaTCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCgat                                                                 >  2:2375158/1‑100 (MQ=255)
                               aaTCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCgat                                                                 >  1:1093938/1‑100 (MQ=255)
                               aaTCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCgat                                                                 <  1:3653606/100‑1 (MQ=255)
                                aTCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCgatg                                                                <  2:688183/100‑1 (MQ=255)
                                aTCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCATTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCgatg                                                                <  2:2458365/100‑1 (MQ=255)
                                 tCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCgatga                                                               >  1:938887/1‑100 (MQ=255)
                                 tCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCgatga                                                               >  1:3875189/1‑100 (MQ=255)
                                     ttCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATctg                                                           <  1:3780389/100‑1 (MQ=255)
                                      tctcAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATctgc                                                          <  1:436984/100‑1 (MQ=255)
                                       ctcAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTCTAGAGAAAGTGTCGAATATAAGGGAATGGACTAGCAGGCTGGCAGCGCGATGATctgct                                                         >  1:444739/1‑100 (MQ=255)
                                          aGCCCCACTAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGtt                                                      <  1:656347/100‑1 (MQ=255)
                                          aGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGtt                                                      <  1:159936/100‑1 (MQ=255)
                                          aGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGtt                                                      <  1:3007484/100‑1 (MQ=255)
                                             cccACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTaga                                                   <  1:3814971/100‑1 (MQ=255)
                                             cccACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTaga                                                   >  2:2494550/1‑100 (MQ=255)
                                              ccACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTagaa                                                  <  1:812889/100‑1 (MQ=255)
                                              ccACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTagaa                                                  <  2:3359779/100‑1 (MQ=255)
                                              ccACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTaga                                                   <  2:1957792/99‑1 (MQ=255)
                                                 cGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAgaagaa                                               >  2:121030/1‑100 (MQ=255)
                                                 cGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAgaagaa                                               >  1:1544722/1‑100 (MQ=255)
                                                   aGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAgg                                             <  2:396352/100‑1 (MQ=255)
                                                    gCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCAGACAGCACGAAGAACTGCTATTAGAAGAAGGt                                            >  1:1444297/1‑100 (MQ=255)
                                                       tggtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGAGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCa                                         >  2:2034310/1‑100 (MQ=255)
                                                        ggtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAg                                        >  1:1487043/1‑100 (MQ=255)
                                                         gtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGt                                       <  1:1739498/100‑1 (MQ=255)
                                                          tggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTa                                      <  1:2001043/100‑1 (MQ=255)
                                                           ggATTCTGCGCACCTTTTGTGTAGTGAAAGCGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTag                                     <  1:243766/100‑1 (MQ=255)
                                                            gATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGAGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTaga                                    <  2:233417/100‑1 (MQ=255)
                                                              tgctgcGCACCTTTGGTGTAGAGAAATATTCGTTTCTTTGGGGCTGGACGTGCAGGTTGGCAGCGTGATGATCTTCTTTTATAAGATGGTGCAGTAGAGc                                  <  2:293067/100‑1 (MQ=255)
                                                              tgctgcGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGc                                  >  1:3044588/1‑100 (MQ=255)
                                                                ctgcGCACCTTTGGTGTAGAGAAAGAGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCAGATAGAAGAAGGAGCAGTaga                                    >  2:2861242/1‑96 (MQ=255)
                                                                 tgcGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGc                               >  2:942459/1‑100 (MQ=255)
                                                                  gcgtACCTTTTTTGTAGAGATAGTGTCGATTCTGTTTGGGGAGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGcc                              <  1:2276162/100‑1 (MQ=255)
                                                                    gcACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCgg                            <  2:3449299/100‑1 (MQ=255)
                                                                    gcACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCgg                            >  1:1945783/1‑100 (MQ=255)
                                                                     cACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGa                           <  1:1189984/100‑1 (MQ=255)
                                                                        cTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAgg                        >  1:2490258/1‑100 (MQ=255)
                                                                            ggTGTAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGagag                    >  1:3170725/1‑100 (MQ=255)
                                                                            ggTGTAGAGAAAGTGGCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGagag                    <  2:3022646/100‑1 (MQ=255)
                                                                             gtgtAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGt                   >  1:562917/1‑100 (MQ=255)
                                                                              tgtAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGtt                  <  2:749887/100‑1 (MQ=255)
                                                                               gtAGAGAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGttt                 <  1:3685817/100‑1 (MQ=255)
                                                                                 agagAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTaa               >  2:1999300/1‑100 (MQ=255)
                                                                                 agagAAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGGAGAGCTGCCGGAAGGACACTTTaa               >  2:3017162/1‑100 (MQ=255)
                                                                                 agagAAAGTGTCGATTCTGGGGGGGTGGAATTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCCGAAGgaga                     >  2:3110966/1‑94 (MQ=255)
                                                                                 agagAAAGTGTCGATTCTGGGGGGATGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGGAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTaa               >  2:1149899/1‑100 (MQ=255)
                                                                                  gagaAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAc              >  1:112348/1‑100 (MQ=255)
                                                                                  gagaAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAc              <  2:1301369/100‑1 (MQ=255)
                                                                                  gagaAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAc              >  2:3338684/1‑100 (MQ=255)
                                                                                   agaAAGTGTCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAcg             <  2:3836821/100‑1 (MQ=255)
                                                                                     aaaGTGTCGAGTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAcgc            <  1:957479/99‑1 (MQ=255)
                                                                                       aGTGTCGATTCTGGTGGGGTGGACTTGCAGGCTGGCAGAGCGATGATCTGCTGTTAGAAGAAGTTGCAGTAGAGCTGCCGGAAGGAGAGGTTAAcgc            >  2:2072469/1‑97 (MQ=255)
                                                                                          gtCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGtt      >  2:1330288/1‑100 (MQ=255)
                                                                                           tCGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGttt     <  2:2914093/100‑1 (MQ=255)
                                                                                            cGATTCTGGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTa    <  2:430565/100‑1 (MQ=255)
                                                                                                     |                                                                                                  
GGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAAT  >  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: