Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,009,732 C→G 16.7% intergenic (+20/‑89) yedF → / → yedK putative TusA family sulfurtransferase/DUF159 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,009,7320CG16.7% 301.7 / 52.8 119intergenic (+20/‑89)yedF/yedKputative TusA family sulfurtransferase/DUF159 family protein
Reads supporting (aligned to +/- strand):  ref base C (47/52);  new base G (7/13);  total (54/65)
Fisher's exact test for biased strand distribution p-value = 3.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.51e-01

Only 100 of 124 total aligned reads displayed.

GGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATC‑CCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACG  >  NC_000913/2009634‑2009831
                                                                                                   |                                                                                                   
ggATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATC‑CCTCATTCCCGCTGTGAcc                                                                                                    >  1:672855/1‑100 (MQ=255)
ggATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATC‑CCTCATTCCCGCTGTGAcc                                                                                                    >  1:3915447/1‑100 (MQ=255)
ggATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATC‑CCTCATTCCCGCTGTGAcc                                                                                                    >  1:1745949/1‑100 (MQ=255)
ggATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACAGTCCTACCATTCGTTATTTAATTCAGAAGTAATC‑CCTCATTCCCGCTGTGAcc                                                                                                    <  2:619777/100‑1 (MQ=255)
     cgcgTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATC‑CCTCATTCCCGCTGTGACCCACAg                                                                                               >  1:3299169/1‑100 (MQ=255)
     cgcgTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATC‑CCTCATTCCCGCTGTGACCCACAg                                                                                               >  2:3247936/1‑100 (MQ=255)
       cgTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATC‑CCTCATTCCCGCTGTGACCCACAGCg                                                                                             <  1:1273315/100‑1 (MQ=255)
       cgTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATC‑CCTCATTCCCGCTGTGACCCACAGCg                                                                                             >  1:2159274/1‑100 (MQ=255)
         tAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATC‑CCTCATTCCCGCTGTGGGTCACAGCggg                                                                                           <  1:1123869/100‑1 (MQ=255)
         tAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATC‑CCTCATTCCCGCTGTGGGTCACAGCggg                                                                                           <  2:1354151/100‑1 (MQ=255)
         tAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATC‑CCTCATTCCCGCTGTGGGTCACAGCggg                                                                                           <  1:241969/100‑1 (MQ=255)
         tAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATC‑CCTCATTCCCGCTGTGGGTCACAGCggg                                                                                           <  1:241849/100‑1 (MQ=255)
         tAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATC‑CCTCATTCCCGCTGTGGGTCACAGCggg                                                                                           <  1:2026830/100‑1 (MQ=255)
         tAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATC‑CCTCATTCCCGCTGTGGGTCACAGCggg                                                                                           <  1:128691/100‑1 (MQ=255)
         tAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATC‑CCTCATTCCCGCTGTGACCCACAGCggg                                                                                           >  1:3094325/1‑100 (MQ=255)
          aaTCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATC‑CCTCATTCCCGCTGTGGGTCACAGCGGGa                                                                                          <  2:653320/100‑1 (MQ=255)
          aaTCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATC‑CCTCATTCCCGCTGTGGGTCACAGCGGGa                                                                                          <  1:2106238/100‑1 (MQ=255)
          aaTCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATC‑CCTCATTCCCGCTGTGGGTCACAGCGGGa                                                                                          <  1:1577062/100‑1 (MQ=255)
          aaTCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATC‑CCTCATTCCCGCTGTGACCCACAGCGGGa                                                                                          >  1:900246/1‑100 (MQ=255)
           aTCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATC‑CCTCATTCCCGCTGTGGGTCACAGCGGGaa                                                                                         <  2:1327044/100‑1 (MQ=255)
           aTCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATC‑CCTCATTCCCGCTGTGACCCACAGCGGGaa                                                                                         >  1:1834169/1‑100 (MQ=255)
            tCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATC‑CCTCATTCCCGCTGTGACCCACAGCGGGAAc                                                                                        <  1:796765/100‑1 (MQ=255)
             cACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATC‑CCTCATTCCCGCTGTGACCCACAGCGGGAACa                                                                                       >  2:796761/1‑100 (MQ=255)
                ggCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATC‑CCTCATTCCCGCTGTGACCCACAGCGGGAACATCg                                                                                    <  2:389789/100‑1 (MQ=255)
                 gCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATC‑CCTCATTCCCGCTGTGACCCACAGCGGGAACATCGt                                                                                   <  1:2861705/100‑1 (MQ=255)
                 gCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATC‑CCTCATTCCCGCTGTGACCCACAGCGGGAACATCGt                                                                                   <  1:1343376/100‑1 (MQ=255)
                  cTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATC‑GCTCATTCCCGCTGTGACCCACAGCGGGAACATCGtt                                                                                  >  2:2227025/1‑100 (MQ=255)
                   tataCGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATC‑CCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTa                                                                                 >  2:3235665/1‑100 (MQ=255)
                   tataCGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATC‑CCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTa                                                                                 <  1:3530695/100‑1 (MQ=255)
                    ataCGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATC‑CCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTAc                                                                                >  1:1197248/1‑100 (MQ=255)
                     cgtcttGCTGAATATCCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATC‑CCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTAcc                                                                               >  1:2718099/6‑100 (MQ=255)
                      gtcttGCTGAATATCCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATC‑CCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCt                                                                              >  2:1048976/5‑100 (MQ=255)
                       cGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATC‑CCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCtt                                                                             >  2:539725/1‑100 (MQ=255)
                         gTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATC‑CCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTAACttt                                                                            >  2:3636105/1‑99 (MQ=255)
                           gCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATC‑CCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTa                                                                         <  1:629449/100‑1 (MQ=255)
                           gCTGGATATTCAGCAAGACGGGCCGACCATTCGATATTTAATTAAGAAGTAATC‑CCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTa                                                                         >  2:1686502/1‑100 (MQ=255)
                            cTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATC‑CCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTAt                                                                        <  2:2901088/100‑1 (MQ=255)
                             tGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATC‑CCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTa                                                                         <  2:3697557/98‑1 (MQ=255)
                             tGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATC‑CCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATc                                                                       <  2:2569755/100‑1 (MQ=255)
                               gATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATC‑CCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATcac                                                                     <  1:1877683/100‑1 (MQ=255)
                               gATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATC‑CCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATcac                                                                     <  1:2442726/100‑1 (MQ=255)
                               gATAGTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATC‑CCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATcac                                                                     <  2:1110290/100‑1 (MQ=255)
                                 tatTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATC‑CCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACAt                                                                   <  2:3419983/100‑1 (MQ=255)
                                  atTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATC‑CCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATa                                                                  <  1:1623284/100‑1 (MQ=255)
                                   ttCAGCAAGAGGGGCCGACCATTCGTTATTTACTTCAGGAGAAAGCGCCTCA‑TAACGGTGTGAGCCACAGCGGGGACAGCGGTACCCTTGGATcagagag                                                                 >  2:3326855/1‑95 (MQ=255)
                                   ttCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATC‑CCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAg                                                                 >  2:2958180/1‑100 (MQ=255)
                                   ttCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATC‑CCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAg                                                                 >  2:1299964/1‑100 (MQ=255)
                                    tCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATC‑CCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGa                                                                <  2:900241/100‑1 (MQ=255)
                                      aGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATC‑CCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGa                                                                >  2:1274534/1‑98 (MQ=255)
                                      aGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATC‑CCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGAtt                                                              >  1:2191944/1‑100 (MQ=255)
                                       gCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATC‑CCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACAGAGAttt                                                             <  1:3272445/100‑1 (MQ=255)
                                         aaGACGGGCCGACACTTCGTTATTTAATTCAGCAGTAGTC‑GCTCCTTCCCGCGGTGACCCACAGCGGGAACATCGGTACCTTTGTATCACATAGAtttt                                                            >  2:1599088/1‑99 (MQ=255)
                                         aaGACGCGCCGACCATTCGTTATTTAATTCAGAAGTAATC‑CCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTa                                                           <  2:3148066/100‑1 (MQ=255)
                                          aGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATC‑CCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTAc                                                          <  1:1188744/100‑1 (MQ=255)
                                               ggCCGACCATTCGTTATTTAATTCAGAAGTAATC‑CCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGcc                                                     <  1:3925098/100‑1 (MQ=255)
                                                gCCGACCATTCGTTATTTAATTCAGAAGTAATC‑CCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGccc                                                    <  2:1997125/100‑1 (MQ=255)
                                                  cGACCATTCGTTATTTAATTCAGAAGTAATC‑CCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGccccc                                                  >  2:2693054/1‑100 (MQ=255)
                                                   gACCATTCGTTATTTAATTCAGAAGTAATC‑CCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCa                                                 >  1:192537/1‑100 (MQ=255)
                                                      cATTCGTTATTTAATTCAGAAGTAATC‑CCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTc                                              >  2:1855176/1‑100 (MQ=255)
                                                       aTTCGTTATTTAATTCAGAAGTAATC‑CCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTcc                                             <  1:115125/100‑1 (MQ=255)
                                                          cGTTATTTAATTCAGAAGTAATC‑CCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTAGCTTTGTATCACATAGATTTTACCCGCCCCCACTCCt                                            >  2:2144575/1‑98 (MQ=255)
                                                          cGTTATTTAATTCAGAAGTAATC‑CCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTAt                                          >  1:195787/1‑100 (MQ=255)
                                                          cGTTATTTAATTCAGAAGTAATC‑CCCCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTa                                           <  1:3141830/99‑1 (MQ=255)
                                                            ttATTTAATTCAGAAGTAATC‑CCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTc                                        >  1:2184063/1‑100 (MQ=255)
                                                            ttATTTAATTCAGAAGTAATC‑CCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTc                                        <  2:514230/100‑1 (MQ=255)
                                                             tATTTAATTCAGAAGTAATC‑CCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCt                                       <  2:3473225/100‑1 (MQ=255)
                                                               tttAATTCAGAAGTAATC‑CCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCttt                                     <  1:581068/100‑1 (MQ=255)
                                                                 tAATTCAGAAGTAATC‑CCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTAGTCTTTTc                                   <  1:2063459/100‑1 (MQ=255)
                                                                   aTTCAGAAGTATTT‑CCTCATTCTCGCTGTGACTCTCAGTGGGAAGATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAg                                 <  2:18853/100‑1 (MQ=255)
                                                                      cAGAAGTAATC‑CCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTAt                              >  1:3176214/1‑100 (MQ=255)
                                                                       aGAAGTAATC‑CCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCTCTCCTATTCTTTTCAGTATg                             <  2:1868695/100‑1 (MQ=255)
                                                                       aGAAGTAATC‑CCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTAt                              >  1:2465044/1‑99 (MQ=255)
                                                                        gAAGTAATC‑CCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGa                            <  1:3671469/100‑1 (MQ=255)
                                                                        gAAGTAATC‑CCGCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTGCCCGCCCCCACTCCTATTCTTTTCAGTATGa                            <  2:1071184/100‑1 (MQ=255)
                                                                        gAAGTAAGC‑CCCCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATGTTACCCGCCCCCACTCCTATTCTTTTCAGTATGa                            >  2:2074833/1‑100 (MQ=255)
                                                                         aaGTAATC‑CCTCATTCCCGCTGTGACCCACAGCTGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGAt                           <  2:1298749/100‑1 (MQ=255)
                                                                          aGTAATC‑CCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGAtc                          <  2:425968/100‑1 (MQ=255)
                                                                          aGTAATC‑CCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGAtc                          >  1:1178015/1‑100 (MQ=255)
                                                                               taacgatgttCCCGCTGTGGGTCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCag                      >  2:142082/9‑100 (MQ=37)
                                                                                c‑ccTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCagga                    >  1:938865/1‑100 (MQ=255)
                                                                                c‑ccTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCagga                    >  1:1955123/1‑100 (MQ=255)
                                                                                  ccTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCTGCCCCCACTCCTATTCTTTTCAGTATGATCTCaggag                   >  1:1564381/1‑100 (MQ=255)
                                                                                  ccTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCaggag                   <  1:3098032/100‑1 (MQ=255)
                                                                                  ccTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCaggag                   <  2:2138298/100‑1 (MQ=255)
                                                                                   gatgttCCCGCTGTGGGTCACAGCGGGAACATCGTTACGTTTGTATCACATAGATTTTAAACGCCCCCACTCCTATTCTTTTCAGTATGATATCagga                    >  2:644831/5‑98 (MQ=37)
                                                                                        tCCCGCTGTGGGTCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCt             >  2:3063121/1‑100 (MQ=255)
                                                                                        tCCCGCTGTGACCCACAGCGGGAACATCGTTCCCTTTATATCACATAGATTTTACGCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCt             <  2:2654111/100‑1 (MQ=255)
                                                                                        tCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTc              <  1:2551977/99‑1 (MQ=255)
                                                                                        tCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCt             <  2:672851/100‑1 (MQ=255)
                                                                                        tCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCt             >  1:1760280/1‑100 (MQ=255)
                                                                                         cccGCTGTGGGTCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTa            >  2:3739152/1‑100 (MQ=255)
                                                                                         cccGCTGTGGGTCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGCTTTCt             >  1:893696/1‑99 (MQ=255)
                                                                                         cccGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTa            >  1:638422/1‑100 (MQ=255)
                                                                                         cccGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTa            >  2:107069/1‑100 (MQ=255)
                                                                                          ccgatctgACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTAt           <  2:2693780/94‑1 (MQ=255)
                                                                                          ccGCTGTGGGTCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTAt           >  1:233115/1‑100 (MQ=255)
                                                                                          acGCTGTGACCCACAGCGGCAACATCGTTACCTTTGTATGACATAGATTTTACCCGCGCCCACTCCTATGCTGTTCAGTATGATCTCaggagggtact             >  2:90160/2‑93 (MQ=255)
                                                                                                tgACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTgg     <  2:338274/100‑1 (MQ=255)
                                                                                                 gACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGa    <  2:832112/100‑1 (MQ=255)
                                                                                                   cccACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACg  <  2:1078713/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
GGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATC‑CCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACG  >  NC_000913/2009634‑2009831

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: ‑