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.G97.5% 283.1 / 0.2 80coding (356/846 nt)sseA3‑mercaptopyruvate sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base . (2/0);  new base G (40/38);  total (42/38)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.66e-01

GGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑GGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAA  >  NC_000913/2652751‑2652942
                                                                                                   |                                                                                              
ggATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑g                                                                                               >  2:338884/1‑100 (MQ=255)
ggATAAGCACCTAATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGCTTCT‑g                                                                                               >  1:1570318/1‑100 (MQ=255)
    aaGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑ggggg                                                                                           <  1:878133/100‑1 (MQ=255)
     aGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑gggggg                                                                                          <  2:1107045/100‑1 (MQ=255)
     aGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑gggggg                                                                                          >  1:2315679/1‑100 (MQ=255)
     aGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑gggggg                                                                                          <  1:1511341/100‑1 (MQ=255)
          cTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑gggggggtggn                                                                                     <  2:986082/100‑6 (MQ=255)
           tGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCGCCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGAc                                                                                    <  2:567482/100‑1 (MQ=255)
                gTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCa                                                                               <  1:1552901/100‑1 (MQ=255)
                 tCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAg                                                                              <  1:804253/100‑1 (MQ=255)
                  cTATGACGAAGGTTATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAgg                                                                             <  1:1383423/100‑1 (MQ=255)
                  cTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAgg                                                                             >  2:405135/1‑100 (MQ=255)
                   tATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGGCGATTCTGGGGGGG‑TGGACTTGCAgg                                                                             >  1:813498/1‑99 (MQ=255)
                     tGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTg                                                                          >  1:216828/1‑100 (MQ=255)
                     tGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTg                                                                          >  1:207102/1‑100 (MQ=255)
                              tAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCga                                                                 <  2:868353/100‑1 (MQ=255)
                               aaTCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCgat                                                                <  2:302054/100‑1 (MQ=255)
                                aTCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCgatg                                                               >  2:2285354/1‑100 (MQ=255)
                                 tCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCgatga                                                              >  1:636089/1‑100 (MQ=255)
                                  cTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCgatgat                                                             >  2:1786433/1‑100 (MQ=255)
                                  cTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCgatgat                                                             >  1:1294769/1‑100 (MQ=255)
                                   ttttCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATc                                                            <  1:842619/100‑1 (MQ=255)
                                   ttttCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATc                                                            >  1:29653/1‑100 (MQ=255)
                                    tttCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATct                                                           <  1:1412377/100‑1 (MQ=255)
                                      tctcAGCCCCACGAGCATGGTGGATGCTGGGCACCTTTGGTGTAGAGAAGGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATctg                                                          <  2:864571/99‑1 (MQ=255)
                                       ctcAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑GGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATctgctg                                                       >  2:679868/1‑100 (MQ=255)
                                        tcAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATctgctg                                                       <  2:248459/100‑1 (MQ=255)
                                        tcAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATctgctg                                                       <  2:1209930/100‑1 (MQ=255)
                                        tcAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATctgctg                                                       <  1:628990/100‑1 (MQ=255)
                                        tcAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATctgctg                                                       <  1:588528/100‑1 (MQ=255)
                                         cAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGt                                                      >  1:1896870/1‑100 (MQ=255)
                                         cAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGt                                                      <  1:2150674/100‑1 (MQ=255)
                                         cAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGt                                                      <  1:405137/100‑1 (MQ=255)
                                          aGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGtt                                                     <  2:1663733/100‑1 (MQ=255)
                                          aGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGtt                                                     <  1:600514/100‑1 (MQ=255)
                                          aGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGtt                                                     >  1:1145984/1‑100 (MQ=255)
                                          aGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGtt                                                     >  1:1711597/1‑100 (MQ=255)
                                          aGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGtt                                                     <  1:1633036/100‑1 (MQ=255)
                                             cccACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTaga                                                  >  2:1194948/1‑100 (MQ=255)
                                             cccACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTaga                                                  >  2:2313988/1‑100 (MQ=255)
                                                   aGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAgg                                            <  2:1390676/100‑1 (MQ=255)
                                                   aGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAgg                                            <  1:688317/100‑1 (MQ=255)
                                                   aGCATGGTCGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCAATGATCTGCTGTTAGAAGAAgg                                            <  2:1488469/100‑1 (MQ=255)
                                                    gCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGt                                           >  2:2343513/1‑100 (MQ=255)
                                                    gCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGt                                           <  1:466177/100‑1 (MQ=255)
                                                      aTGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGc                                         <  2:1363652/100‑1 (MQ=255)
                                                      aTGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGc                                         <  1:641994/100‑1 (MQ=255)
                                                        ggtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAg                                       >  1:265100/1‑100 (MQ=255)
                                                           ggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTag                                    <  1:1090664/100‑1 (MQ=255)
                                                           ggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTag                                    <  1:983416/100‑1 (MQ=255)
                                                            gATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTaga                                   <  1:142069/100‑1 (MQ=255)
                                                             aTGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTagag                                  <  1:953552/100‑1 (MQ=255)
                                                             aTGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTagag                                  >  1:1873839/1‑100 (MQ=255)
                                                               gctgcGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCt                                <  1:716706/100‑1 (MQ=255)
                                                                ctgcGCACCTTTGGTGTAGGGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTg                               <  1:1941001/100‑1 (MQ=255)
                                                                 tgcGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGc                              >  2:1112985/1‑100 (MQ=255)
                                                                 tgcGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGc                              >  1:1208844/1‑100 (MQ=255)
                                                                  gcgcACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGcc                             <  2:2159421/100‑1 (MQ=255)
                                                                  gcgcACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGcc                             >  2:752418/1‑100 (MQ=255)
                                                                  gcgcACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGcc                             <  1:1112990/100‑1 (MQ=255)
                                                                   cgcACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TTGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCg                            <  2:549413/100‑1 (MQ=255)
                                                                       ccTTTGGTGTAGACAAAGTGCAGATTCTGGGGGGG‑CGGAATTGCAGCCTGGCAGACCGATGATATGCTGTTAGAAGTAGGTGCAGTAGAGCTGcag                            >  2:443919/1‑94 (MQ=255)
                                                                           tggtggAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGgaga                    <  2:210319/100‑1 (MQ=255)
                                                                           tGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGgaga                    >  1:352747/1‑100 (MQ=255)
                                                                            ggTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGagag                   >  2:2184127/1‑100 (MQ=255)
                                                                             gtgtAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGagag                   >  1:716564/1‑99 (MQ=255)
                                                                              tgtAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGGCGGAAGGagag                   >  2:2105596/1‑98 (MQ=255)
                                                                              tgtAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGtt                 >  2:89389/1‑100 (MQ=255)
                                                                              tgtAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGtt                 >  1:1762654/1‑100 (MQ=255)
                                                                               gtAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGttt                >  2:2055406/1‑100 (MQ=255)
                                                                                tAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTa               >  1:1471687/1‑100 (MQ=255)
                                                                                tAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTa               >  2:1117548/1‑100 (MQ=255)
                                                                                 agagAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTaa              <  1:2314009/100‑1 (MQ=255)
                                                                                 agagAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTa               >  1:1196776/1‑99 (MQ=255)
                                                                                 agagAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGGTTaa              >  2:1404590/1‑100 (MQ=255)
                                                                                  gagaAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAc             >  2:1019440/1‑100 (MQ=255)
                                                                                   agaAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAcg            >  2:1503180/1‑100 (MQ=255)
                                                                                       aGTGTCGATTCTGGGGGGG‑TGGACTTGCTGGCTGGCAGCGAGATGATCTGCGGTTAGAAGAAGGTGCAGAAGAGCGGCCGGAAGGAGAGTTTAAcgccgc        >  2:2270227/1‑100 (MQ=255)
                                                                                       aGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAcgccgc        <  2:207101/100‑1 (MQ=255)
                                                                                          ntCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGtt     >  1:305206/2‑100 (MQ=255)
                                                                                          gtCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGtt     <  1:1137353/100‑1 (MQ=255)
                                                                                          gtCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCCGGTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTGAACGCCGCGt      >  2:913039/1‑99 (MQ=255)
                                                                                           tCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGttt    <  1:1116208/100‑1 (MQ=255)
                                                                                            caaTTCT‑GGGGGGATGGACTTGCAGGCTGGCAGCGCGATGATTTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTa   >  1:1615675/3‑100 (MQ=255)
                                                                                            cGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTa   >  1:846923/1‑100 (MQ=255)
                                                                                            cGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAAGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTa   >  2:1669593/1‑100 (MQ=255)
                                                                                             gATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTaa  >  1:2337784/1‑100 (MQ=255)
                                                                                                   |                                                                                              
GGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑GGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAA  >  NC_000913/2652751‑2652942

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