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.4% 260.8 / 12.4 76coding (356/846 nt)sseA3‑mercaptopyruvate sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base . (2/3);  new base G (28/43);  total (30/46)
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.98e-01

AGGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑GG‑GG‑G‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAAT  >  NC_000913/2652750‑2652943
                                                                                                    |                                                                                                 
aGGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCt                                                                                                    >  2:846675/1‑100 (MQ=255)
aGGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCt                                                                                                    <  2:82035/100‑1 (MQ=255)
  gATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATTCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑gg                                                                                                 <  2:1477403/100‑1 (MQ=255)
  gATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCttg                                                                                                  <  2:1336353/100‑3 (MQ=255)
   aTAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑gg‑g                                                                                               <  1:2299066/100‑1 (MQ=255)
   aTAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑gg                                                                                                 <  2:2124504/99‑1 (MQ=255)
    tAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑gg‑gg                                                                                              >  1:1059832/1‑100 (MQ=255)
    tAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑gg‑gg                                                                                              >  1:53441/1‑100 (MQ=255)
      aGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑gg‑gg‑g‑g                                                                                          >  2:65920/1‑100 (MQ=255)
      aGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑gg‑gg‑g‑g                                                                                          >  1:1739570/1‑100 (MQ=255)
      aGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑gg‑gg‑g‑g                                                                                          >  1:843039/1‑100 (MQ=255)
       ncACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGCTTCT‑gg‑gg‑g‑gg                                                                                         >  2:1041711/2‑99 (MQ=255)
       ggaCCTGATTGTCGATGACGAAGGTAATCTTTTCTCAGCCCCGCGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑gg‑gg‑g‑gg                                                                                         >  1:1870296/3‑99 (MQ=255)
       gCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑gg‑gg‑g‑gg                                                                                         >  1:1010115/1‑99 (MQ=255)
        cACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑gg‑gg‑g‑ggt                                                                                        >  1:289473/1‑98 (MQ=255)
              aTTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GG‑G‑GTGGACtt                                                                                   >  2:173636/1‑100 (MQ=255)
              aTTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GG‑G‑GTGGACtt                                                                                   <  1:2381681/100‑1 (MQ=255)
                tGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GG‑G‑GTGGACTTGc                                                                                 <  2:144121/100‑1 (MQ=255)
                  tCTATGACGAAGGTAATCTTTTCTCAGCCCTACGAGCATTGTCGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑GGTGG‑G‑GTGGACTTGCAg                                                                               >  2:2272995/1‑100 (MQ=255)
                     aTGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GG‑G‑GTGGACTTGCAGGCt                                                                            >  2:717995/1‑100 (MQ=255)
                      tGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GG‑G‑GTGGACTTGCAGGCTg                                                                           <  2:657248/100‑1 (MQ=255)
                      tGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑GG‑GG‑G‑GTGGACTTGCAGGCTgg                                                                          <  2:718225/100‑1 (MQ=255)
                       gACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GG‑G‑GTGGACTTGCAGGCTgg                                                                          <  1:2674070/100‑1 (MQ=255)
                        tcGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GG‑G‑GTGGACTTGCAGGCTGGc                                                                         <  1:91778/99‑1 (MQ=255)
                        aCGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GG‑G‑GTGGACTTGCAGGCTGGc                                                                         >  1:247536/1‑100 (MQ=255)
                        aCGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GG‑G‑GTGGACTTGCAGGCTGGc                                                                         <  2:843033/100‑1 (MQ=255)
                        aCGAAGGGAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GG‑G‑GTGGACTTGCAGGCTGGc                                                                         <  2:560886/100‑1 (MQ=255)
                         cGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GG‑G‑GTGGACTTGCAGGCTGGCa                                                                        >  1:1391000/1‑100 (MQ=255)
                         cGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GG‑G‑GTGGACTTGCAGGCTGGCa                                                                        >  2:91776/1‑100 (MQ=255)
                            aGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GG‑G‑GTGGACTTGCAGGCTGGCAgcg                                                                     >  1:450272/1‑100 (MQ=255)
                             ggTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GG‑G‑GTGGACTTGCAGGCTGGCAgcgc                                                                    <  2:1328255/100‑1 (MQ=255)
                              gTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GG‑G‑GTGGACTTGCAGGCTGGCAgcgc                                                                    <  2:1327970/99‑1 (MQ=255)
                               tAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GG‑G‑GTGGACTTGCAGGCTGGCAGCGCga                                                                  >  1:1639795/1‑100 (MQ=255)
                                  tCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GG‑G‑GTGGACTTGCAGGCTGGCAGCGCgatga                                                               <  2:1301161/100‑1 (MQ=255)
                                  tCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCGGGG‑GG‑G‑GTGGACTTGCAGGCTGGCAGCGCgatga                                                               <  1:2327325/100‑1 (MQ=255)
                                     tttCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GG‑G‑GTGGACTTGCAGGCTGGCAGCGCGATGATct                                                            <  1:708275/100‑1 (MQ=255)
                                     tttCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GG‑G‑GTGGACTTGCAGGCTGGCAGCGCGATGATct                                                            <  1:754043/100‑1 (MQ=255)
                                        ctcAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GG‑G‑GTGGGCTTGCAGGCTGGCAGCGCGATGATctgct                                                         >  1:1835389/1‑100 (MQ=255)
                                        ctcAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GG‑G‑GTGGACTTGCAGGCTGGGAGCGCGATGATctgct                                                         >  2:1613614/1‑100 (MQ=255)
                                        ctcAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GG‑G‑GTGGACTTGCAGGCTGGCAGCGCGATGATctgct                                                         >  2:2452855/1‑100 (MQ=255)
                                        ctcAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GG‑G‑GTGGACTTGCAGGCTGGCAGCGCGATGATctgct                                                         >  1:492499/1‑100 (MQ=255)
                                         tcAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GG‑G‑GTGGACTTGCAGGCTGGCAGCGCGATGATctgcag                                                        >  2:1678619/1‑98 (MQ=255)
                                          cAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GG‑G‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGt                                                       <  1:847741/100‑1 (MQ=255)
                                           aGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GG‑G‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGtt                                                      <  2:823659/100‑1 (MQ=255)
                                           aGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑GG‑GGCG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGtt                                                      <  2:1699307/100‑1 (MQ=255)
                                             ccccACTAGCATGGTTGTGGCTGTTCACCGTTGGTGTAGAGAAAGTGTCGTTTCTGGG‑GG‑G‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTag                                                    <  1:195701/100‑1 (MQ=255)
                                             ccccACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GG‑G‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTag                                                    <  2:2661971/100‑1 (MQ=255)
                                              cccACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GG‑G‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTaga                                                   <  2:1756045/100‑1 (MQ=255)
                                                cACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GG‑G‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTagaag                                                 >  2:1840863/1‑100 (MQ=255)
                                                cACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GG‑G‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTagaag                                                 <  2:910031/100‑1 (MQ=255)
                                                    aGCATGGTGGGTGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GG‑G‑GTGGACTTGCAGACTGGCAGCGCGATGATCTGCTGTTAGAAGAAgg                                             <  2:1809456/100‑1 (MQ=255)
                                                    aGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GG‑G‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAgg                                             <  1:2638470/100‑1 (MQ=255)
                                                     gCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GG‑G‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGt                                            >  2:2715667/1‑100 (MQ=255)
                                                       ntggtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GG‑G‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGc                                          >  2:965797/2‑100 (MQ=255)
                                                       aTGGTTGATGCTGCGCACACGTGGTGTAGAGAGAGTGTCGATTCTGGG‑GG‑G‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGc                                          <  2:247331/100‑1 (MQ=255)
                                                       aTGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GG‑G‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGc                                          <  1:1442792/100‑1 (MQ=255)
                                                       aTGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GG‑G‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGc                                          <  2:21533/100‑1 (MQ=255)
                                                        tggtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GG‑G‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCa                                         <  2:884033/100‑1 (MQ=255)
                                                         ggtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GG‑G‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAg                                        >  2:558703/1‑100 (MQ=255)
                                                          gtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑GG‑GG‑G‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTa                                      <  1:2098821/100‑1 (MQ=255)
                                                              aTGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GG‑G‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTagag                                   <  2:2269695/100‑1 (MQ=255)
                                                               tgctgcGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GG‑G‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTagag                                   <  1:883062/99‑1 (MQ=255)
                                                               tgctgcGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GG‑G‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGc                                  <  1:1774821/100‑1 (MQ=255)
                                                               tgctgcGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GG‑G‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGc                                  >  1:824329/1‑100 (MQ=255)
                                                               tgctgcGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GG‑G‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGc                                  <  2:1568118/100‑1 (MQ=255)
                                                                 ctgcGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GG‑G‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTg                                <  1:965803/100‑1 (MQ=255)
                                                                  tgcGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GG‑G‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGc                               >  2:1146985/1‑100 (MQ=255)
                                                                   gcgcACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GG‑G‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGcc                              >  2:875883/1‑100 (MQ=255)
                                                                   gcgcACCTTTGGTGTAGAGAAAGTGTCGATTCT‑GG‑GG‑GTGTGGACTTGCAGGCTGGTAGCGCGATGATCTGCGGTTAGAAGAAGGTGCAGTAGAGCt                                 >  2:2627619/1‑97 (MQ=255)
                                                                      cACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GG‑G‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGa                           >  1:156281/1‑100 (MQ=255)
                                                                         cTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GG‑G‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAgg                        >  1:769993/1‑100 (MQ=255)
                                                                         cTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GG‑G‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAgg                        >  2:843808/1‑100 (MQ=255)
                                                                          tttGGTGTAGAGAAAGTGTCGATTCTGGG‑GG‑G‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGga                       <  1:718000/100‑1 (MQ=255)
                                                                            tGGTGTAGAGAAAGTGTCGATTCTGGG‑GG‑G‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGGAGTAGAGCTGCCGGAAGgaga                     >  2:2654127/1‑100 (MQ=255)
                                                                            tGGTGTAGAGAAAGTGTCGATTCTGGG‑GG‑G‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGgaga                     <  1:2654147/100‑1 (MQ=255)
                                                                             ggTGTAGAGAAAGTGTCGATTCTGGG‑GG‑G‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGagag                    <  1:846681/100‑1 (MQ=255)
                                                                              gtgtAGAGAAAGTGTCGATTCTGGG‑GG‑G‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGagag                    <  1:2482344/99‑1 (MQ=255)
                                                                              gtgtAGAGAAAGTGTCGATTCTGGG‑GG‑G‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGt                   <  2:1870281/100‑1 (MQ=255)
                                                                              gtgtAGAGAAAGTGTCGATTCTGGG‑GG‑G‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGt                   <  2:1010108/100‑1 (MQ=255)
                                                                                 tAGAGAAAGTGTCGATTCTGGG‑GG‑G‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTa                <  1:173639/100‑1 (MQ=255)
                                                                                 tAGAGAAAGTGTCGATTCTGGG‑GG‑G‑GTGGACTTGCAGGCTGGAAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTa                <  2:684953/100‑1 (MQ=255)
                                                                                   gagaAAGTGTCGATTCTGGG‑GG‑G‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAc              <  2:492495/100‑1 (MQ=255)
                                                                                     gaAAGTGTCGATTCTGGG‑GG‑G‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAcgc            >  1:920869/1‑100 (MQ=255)
                                                                                          tgtCGATTCTGGG‑GG‑G‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGt       >  2:384659/1‑100 (MQ=255)
                                                                                            tCGATTCTGGG‑GG‑G‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGttt     <  1:2715686/100‑1 (MQ=255)
                                                                                            tCGATTCTGGG‑GG‑G‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGttt     >  1:495910/1‑100 (MQ=255)
                                                                                            tCGATTCTGGG‑GG‑G‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGttt     <  1:523802/100‑1 (MQ=255)
                                                                                            tCGATTCTGGG‑GG‑G‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGttt     <  1:2010427/100‑1 (MQ=255)
                                                                                               aTTCTGGG‑GG‑G‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAAt  >  1:1544181/1‑100 (MQ=255)
                                                                                               aTTCTGGG‑GG‑G‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAAt  >  1:2154871/1‑100 (MQ=255)
                                                                                               aTTCTGGG‑GG‑G‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAAt  <  1:558707/100‑1 (MQ=255)
                                                                                                    |                                                                                                 
AGGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑GG‑GG‑G‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAAT  >  NC_000913/2652750‑2652943

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