Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,652,855 (G)6→7 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.G93.6% 383.0 / 20.3 110coding (356/846 nt)sseA3‑mercaptopyruvate sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base . (3/4);  new base G (48/55);  total (51/59)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

Only 100 of 121 total aligned reads displayed.

GATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑GG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAAT  >  NC_000913/2652752‑2652943
                                                                                                  |                                                                                                 
gATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑gg                                                                                                 >  2:1769192/1‑100 (MQ=255)
gATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑gg                                                                                                 >  2:2215833/1‑100 (MQ=255)
  tAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑gg‑g‑g                                                                                             <  2:1171052/100‑1 (MQ=255)
      cACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑gg‑g‑ggggt                                                                                         <  1:1851251/100‑3 (MQ=255)
        ccTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑gg‑g‑ggggtgg                                                                                       >  1:2119415/1‑96 (MQ=255)
        ccTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑gg‑g‑ggggtgg                                                                                       >  1:2850622/1‑96 (MQ=255)
         cTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑gg‑g‑ggggtgga                                                                                      >  1:2209647/1‑95 (MQ=255)
         cTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑gg‑g‑ggggtgga                                                                                      >  2:599741/1‑95 (MQ=255)
           gATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑G‑GGG‑TGGACt                                                                                    <  1:65147/100‑1 (MQ=255)
           gATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑G‑GGG‑TGGACt                                                                                    <  1:127971/100‑1 (MQ=255)
              tGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑G‑GGG‑TGGACTTGc                                                                                 >  1:3028894/1‑100 (MQ=255)
              tGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑G‑GGG‑TGGACTTGc                                                                                 <  2:2035215/100‑1 (MQ=255)
               gTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑G‑GGG‑TGGACTTGCa                                                                                <  1:599742/100‑1 (MQ=255)
               gTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑G‑GGG‑TGGACTTGCa                                                                                >  1:2392275/1‑100 (MQ=255)
                tCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑G‑GGG‑TGGATTTGCAg                                                                               >  1:2628381/1‑100 (MQ=255)
                 cTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑G‑GGG‑TGGACTTGCAgg                                                                              <  2:2062248/100‑1 (MQ=255)
                  tATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑G‑GGG‑TGGACtagatcgg                                                                             >  1:1462597/1‑93 (MQ=255)
                  tATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑G‑GGG‑TGGACTTGCAGGc                                                                             <  2:1508684/100‑1 (MQ=255)
                  tATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑G‑GGG‑TGGACTTGCAGGc                                                                             >  2:735636/1‑100 (MQ=255)
                    tGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑G‑GGG‑TGGACTTGCAGGCTg                                                                           <  1:1426587/100‑1 (MQ=255)
                    tGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑G‑GGG‑TGGACTTGCAGGCTg                                                                           <  2:761724/100‑1 (MQ=255)
                     gACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑G‑GGG‑TGGACTTGCAGGccgg                                                                          >  2:2324582/1‑97 (MQ=255)
                     gACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑G‑GGG‑TGGACTTGCAGGCTgg                                                                          >  1:227544/1‑100 (MQ=255)
                     gACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑G‑GGG‑TGGACTTGCAGGCTgg                                                                          >  1:2422802/1‑100 (MQ=255)
                     gACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑G‑GGG‑TGGACTTGCAGGCTgg                                                                          <  1:2203868/100‑1 (MQ=255)
                      aCGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑G‑GGG‑TGGACTTGCAGGCTGGc                                                                         <  1:1450762/100‑1 (MQ=255)
                         aaGGTAATCTTTTCTCAGCGCCACGAGCATGTTGGATGCTGCGCACCTTTGGTGTAGAGTAAGTGTCGATTCT‑GGTG‑GGG‑TGGACTTGCAGGCTGGc                                                                         <  2:3117052/97‑1 (MQ=255)
                         aaGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAgc                                                                      >  2:1407908/1‑100 (MQ=255)
                          aGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAgcg                                                                     <  1:571851/100‑1 (MQ=255)
                           ggTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑GG‑G‑GGGTTGGACTTGCAGGCTGGCAgcgc                                                                    <  2:1974996/100‑1 (MQ=255)
                                tCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCgatga                                                               >  2:1407204/1‑100 (MQ=255)
                                tCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑G‑GGG‑AGGACTTGCAGGCTGGCAGCGCgctga                                                               >  1:3170552/1‑100 (MQ=255)
                                 cTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCgatgat                                                              >  1:2883850/1‑100 (MQ=255)
                                 cTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCgatgat                                                              <  1:1853463/100‑1 (MQ=255)
                                      ctcAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATctgct                                                         >  2:1221490/1‑100 (MQ=255)
                                      ctcAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATctgct                                                         <  1:1236817/100‑1 (MQ=255)
                                       tcAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATctgctg                                                        >  1:564751/1‑100 (MQ=255)
                                        cAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGt                                                       <  2:1937127/100‑1 (MQ=255)
                                        cAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGt                                                       >  1:1816858/1‑100 (MQ=255)
                                          gCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTa                                                     >  1:2869995/1‑100 (MQ=255)
                                          gCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTa                                                     <  1:1221493/100‑1 (MQ=255)
                                          gCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTa                                                     <  1:1322025/100‑1 (MQ=255)
                                          gCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATGCTGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTa                                                     >  2:304375/1‑100 (MQ=255)
                                            cccACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTaga                                                   >  1:3084392/1‑100 (MQ=255)
                                            cccACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTaga                                                   >  1:2872993/1‑100 (MQ=255)
                                              cACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTagaag                                                 <  1:970396/100‑1 (MQ=255)
                                              cACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTagaag                                                 <  2:564750/100‑1 (MQ=255)
                                              cACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑G‑GGG‑TGGACGTGCAGGCTGGCAGCGCGATGATCTGCTGTTagaag                                                 <  2:844485/100‑1 (MQ=255)
                                               aCGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTagaaga                                                <  2:2034800/100‑1 (MQ=255)
                                                cGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGACAGTGTCGAGTCTGGG‑G‑GGG‑TGGCCTTGCAGGATGGCAGCGCGATGATGTGCTGTTAgaagaa                                               >  2:1568144/1‑100 (MQ=255)
                                                cGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAgaagaa                                               >  2:906496/1‑100 (MQ=255)
                                                  aGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑G‑GGG‑TGGATTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAgg                                             <  2:2628373/100‑1 (MQ=255)
                                                  aGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAgg                                             <  1:2578333/100‑1 (MQ=255)
                                                       ggtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAg                                        <  2:3084380/100‑1 (MQ=255)
                                                       ggtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAg                                        >  2:2340265/1‑100 (MQ=255)
                                                         tggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTa                                      >  1:692117/1‑100 (MQ=255)
                                                         tggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTa                                      <  1:2497779/100‑1 (MQ=255)
                                                               ctgcGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTg                                >  1:2323539/1‑100 (MQ=255)
                                                               ctgcGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑GG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGc                               >  1:2444844/1‑100 (MQ=255)
                                                                  cgcACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGcc                              <  2:1590447/99‑1 (MQ=255)
                                                                  cgcACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCg                             <  2:2422794/100‑1 (MQ=255)
                                                                    cACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGa                           >  2:810845/1‑100 (MQ=255)
                                                                      ccTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAg                         <  2:3342145/100‑1 (MQ=255)
                                                                      ccTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAg                         <  2:923928/100‑1 (MQ=255)
                                                                      ccTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAg                         >  2:112762/1‑100 (MQ=255)
                                                                      ccTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAg                         >  2:1811578/1‑100 (MQ=255)
                                                                       cTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAgg                        >  2:401294/1‑100 (MQ=255)
                                                                       cTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAgg                        <  1:1097411/100‑1 (MQ=255)
                                                                       cTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAgg                        <  1:1097443/100‑1 (MQ=255)
                                                                       cTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAgg                        <  1:3286301/100‑1 (MQ=255)
                                                                        tttGGTGTAGAGAAAGTGTCGATTCTGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGga                       >  1:1244086/1‑100 (MQ=255)
                                                                        tttGGTGTAGAGAAAGTGTCGATTCT‑GG‑G‑GGG‑TGGACTTGAAGGCTGGCAGCGCGATGCTCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGgag                      >  2:140698/1‑100 (MQ=255)
                                                                         ttGGTGTAGAGAAAGTGTCGATTCTGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGgag                      >  1:1891772/1‑100 (MQ=255)
                                                                          tGGTGTAGAGAAAGTGTCGATTCT‑GG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGagag                    <  2:410123/100‑1 (MQ=255)
                                                                           ggTGTAGAGAAAGTGTCGATTCTGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGagag                    <  1:1049674/100‑1 (MQ=255)
                                                                           ggTGTAGAGAAAGTGTCGATTCTGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGagag                    >  2:766438/1‑100 (MQ=255)
                                                                            gtgtAGAGAAAGTGTCGATTCTGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGagag                    >  2:5423/1‑99 (MQ=255)
                                                                            gtgtAGAGAAAGTGTCGATTCTGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGt                   <  1:953552/100‑1 (MQ=255)
                                                                             tgtAGAGAAAGTGTCGATTCTGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGtt                  >  2:2879199/1‑100 (MQ=255)
                                                                             tgtAGAGAAAGTGTCGATTCTGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGtt                  >  1:1033214/1‑100 (MQ=255)
                                                                             tgtAGAGAAAGTGTCGATTCTGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGtt                  <  2:391388/100‑1 (MQ=255)
                                                                             tgtAGAGAAAGTGTCGATTCT‑GG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGttt                 <  2:2444836/100‑1 (MQ=255)
                                                                               tAGAGAAAGTGTCGATTCTGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGgaagatcgg                >  2:1394749/1‑93 (MQ=255)
                                                                               tAGAGAAAGTGTCGATTCTGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTa                >  2:2027156/1‑100 (MQ=255)
                                                                                agagAAAGTGTCGATTCTGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTaa               >  2:2553000/1‑100 (MQ=255)
                                                                                agagAAAGTGTCGATTCTGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTaa               >  2:281000/1‑100 (MQ=255)
                                                                                agagAAAGTGTCGATTCTGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTaa               >  1:1466621/1‑100 (MQ=255)
                                                                                  agaAAGTGTCGATTCTGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAcg             <  1:112761/100‑1 (MQ=255)
                                                                                  agaAAGTGTCGATTCTGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAcg             <  1:1811582/100‑1 (MQ=255)
                                                                                   gaAAGTGTCGATTCTGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAcgc            <  2:1850771/100‑1 (MQ=255)
                                                                                     aaGTGTCGATTCTGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAcgccg          <  2:2978744/100‑1 (MQ=255)
                                                                                       gtgtCGATTCTGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCcgcg        >  1:2652402/1‑100 (MQ=255)
                                                                                       gtgtCGATTCTGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCcgcg        <  2:2323532/100‑1 (MQ=255)
                                                                                        tgtCGATTCTGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCcgcg        >  2:1396972/1‑99 (MQ=255)
                                                                                        tatCGATTCT‑GG‑GCGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGt       >  1:2303986/3‑100 (MQ=255)
                                                                                         gtCGATTCTGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGtt      >  1:2175042/1‑100 (MQ=255)
                                                                                         gtCGATTCTGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGtt      <  1:1407911/100‑1 (MQ=255)
                                                                                           cGATTCTGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTa    >  2:2771414/1‑100 (MQ=255)
                                                                                             aTTCTGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAAt  <  1:3320280/100‑1 (MQ=255)
                                                                                             aTTCTGGG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAAt  <  2:1932634/100‑1 (MQ=255)
                                                                                                  |                                                                                                 
GATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑GG‑G‑GGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAAT  >  NC_000913/2652752‑2652943

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