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.

GGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACG  >  NC_000913/2009634‑2009831
                                                                                                   |                                                                                                   
ggATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGAcc                                                                                                    >  1:672855/1‑100 (MQ=255)
ggATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGAcc                                                                                                    >  1:3915447/1‑100 (MQ=255)
ggATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGAcc                                                                                                    >  1:1745949/1‑100 (MQ=255)
ggATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACAGTCCTACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGAcc                                                                                                    <  2:619777/100‑1 (MQ=255)
     cgcgTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAg                                                                                               >  1:3299169/1‑100 (MQ=255)
     cgcgTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAg                                                                                               >  2:3247936/1‑100 (MQ=255)
       cgTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCg                                                                                             <  1:1273315/100‑1 (MQ=255)
       cgTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCg                                                                                             >  1:2159274/1‑100 (MQ=255)
         tAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGGGTCACAGCggg                                                                                           <  1:1123869/100‑1 (MQ=255)
         tAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGGGTCACAGCggg                                                                                           <  2:1354151/100‑1 (MQ=255)
         tAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGGGTCACAGCggg                                                                                           <  1:241969/100‑1 (MQ=255)
         tAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGGGTCACAGCggg                                                                                           <  1:241849/100‑1 (MQ=255)
         tAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGGGTCACAGCggg                                                                                           <  1:2026830/100‑1 (MQ=255)
         tAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGGGTCACAGCggg                                                                                           <  1:128691/100‑1 (MQ=255)
         tAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCggg                                                                                           >  1:3094325/1‑100 (MQ=255)
          aaTCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGGGTCACAGCGGGa                                                                                          <  2:653320/100‑1 (MQ=255)
          aaTCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGGGTCACAGCGGGa                                                                                          <  1:2106238/100‑1 (MQ=255)
          aaTCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGGGTCACAGCGGGa                                                                                          <  1:1577062/100‑1 (MQ=255)
          aaTCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGa                                                                                          >  1:900246/1‑100 (MQ=255)
           aTCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGGGTCACAGCGGGaa                                                                                         <  2:1327044/100‑1 (MQ=255)
           aTCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGaa                                                                                         >  1:1834169/1‑100 (MQ=255)
            tCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAAc                                                                                        <  1:796765/100‑1 (MQ=255)
             cACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACa                                                                                       >  2:796761/1‑100 (MQ=255)
                ggCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCg                                                                                    <  2:389789/100‑1 (MQ=255)
                 gCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGt                                                                                   <  1:2861705/100‑1 (MQ=255)
                 gCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGt                                                                                   <  1:1343376/100‑1 (MQ=255)
                  cTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCGCTCATTCCCGCTGTGACCCACAGCGGGAACATCGtt                                                                                  >  2:2227025/1‑100 (MQ=255)
                   tataCGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTa                                                                                 >  2:3235665/1‑100 (MQ=255)
                   tataCGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTa                                                                                 <  1:3530695/100‑1 (MQ=255)
                    ataCGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTAc                                                                                >  1:1197248/1‑100 (MQ=255)
                     cgtcttGCTGAATATCCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTAcc                                                                               >  1:2718099/6‑100 (MQ=255)
                      gtcttGCTGAATATCCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCt                                                                              >  2:1048976/5‑100 (MQ=255)
                       cGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCtt                                                                             >  2:539725/1‑100 (MQ=255)
                         gTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTAACttt                                                                            >  2:3636105/1‑99 (MQ=255)
                           gCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTa                                                                         <  1:629449/100‑1 (MQ=255)
                           gCTGGATATTCAGCAAGACGGGCCGACCATTCGATATTTAATTAAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTa                                                                         >  2:1686502/1‑100 (MQ=255)
                            cTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTAt                                                                        <  2:2901088/100‑1 (MQ=255)
                             tGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTa                                                                         <  2:3697557/98‑1 (MQ=255)
                             tGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATc                                                                       <  2:2569755/100‑1 (MQ=255)
                               gATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATcac                                                                     <  1:1877683/100‑1 (MQ=255)
                               gATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATcac                                                                     <  1:2442726/100‑1 (MQ=255)
                               gATAGTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATcac                                                                     <  2:1110290/100‑1 (MQ=255)
                                 tatTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACAt                                                                   <  2:3419983/100‑1 (MQ=255)
                                  atTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATa                                                                  <  1:1623284/100‑1 (MQ=255)
                                   ttCAGCAAGAGGGGCCGACCATTCGTTATTTACTTCAGGAGAAAGCGCCTCATAACGGTGTGAGCCACAGCGGGGACAGCGGTACCCTTGGATcagagag                                                                 >  2:3326855/1‑95 (MQ=255)
                                   ttCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAg                                                                 >  2:2958180/1‑100 (MQ=255)
                                   ttCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAg                                                                 >  2:1299964/1‑100 (MQ=255)
                                    tCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGa                                                                <  2:900241/100‑1 (MQ=255)
                                      aGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGa                                                                >  2:1274534/1‑98 (MQ=255)
                                      aGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGAtt                                                              >  1:2191944/1‑100 (MQ=255)
                                       gCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACAGAGAttt                                                             <  1:3272445/100‑1 (MQ=255)
                                         aaGACGGGCCGACACTTCGTTATTTAATTCAGCAGTAGTCGCTCCTTCCCGCGGTGACCCACAGCGGGAACATCGGTACCTTTGTATCACATAGAtttt                                                            >  2:1599088/1‑99 (MQ=255)
                                         aaGACGCGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTa                                                           <  2:3148066/100‑1 (MQ=255)
                                          aGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTAc                                                          <  1:1188744/100‑1 (MQ=255)
                                               ggCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGcc                                                     <  1:3925098/100‑1 (MQ=255)
                                                gCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGccc                                                    <  2:1997125/100‑1 (MQ=255)
                                                  cGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGccccc                                                  >  2:2693054/1‑100 (MQ=255)
                                                   gACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCa                                                 >  1:192537/1‑100 (MQ=255)
                                                      cATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTc                                              >  2:1855176/1‑100 (MQ=255)
                                                       aTTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTcc                                             <  1:115125/100‑1 (MQ=255)
                                                          cGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTAGCTTTGTATCACATAGATTTTACCCGCCCCCACTCCt                                            >  2:2144575/1‑98 (MQ=255)
                                                          cGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTAt                                          >  1:195787/1‑100 (MQ=255)
                                                          cGTTATTTAATTCAGAAGTAATCCCCCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTa                                           <  1:3141830/99‑1 (MQ=255)
                                                            ttATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTc                                        >  1:2184063/1‑100 (MQ=255)
                                                            ttATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTc                                        <  2:514230/100‑1 (MQ=255)
                                                             tATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCt                                       <  2:3473225/100‑1 (MQ=255)
                                                               tttAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCttt                                     <  1:581068/100‑1 (MQ=255)
                                                                 tAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTAGTCTTTTc                                   <  1:2063459/100‑1 (MQ=255)
                                                                   aTTCAGAAGTATTTCCTCATTCTCGCTGTGACTCTCAGTGGGAAGATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAg                                 <  2:18853/100‑1 (MQ=255)
                                                                      cAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTAt                              >  1:3176214/1‑100 (MQ=255)
                                                                       aGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCTCTCCTATTCTTTTCAGTATg                             <  2:1868695/100‑1 (MQ=255)
                                                                       aGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTAt                              >  1:2465044/1‑99 (MQ=255)
                                                                        gAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGa                            <  1:3671469/100‑1 (MQ=255)
                                                                        gAAGTAATCCCGCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTGCCCGCCCCCACTCCTATTCTTTTCAGTATGa                            <  2:1071184/100‑1 (MQ=255)
                                                                        gAAGTAAGCCCCCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATGTTACCCGCCCCCACTCCTATTCTTTTCAGTATGa                            >  2:2074833/1‑100 (MQ=255)
                                                                         aaGTAATCCCTCATTCCCGCTGTGACCCACAGCTGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGAt                           <  2:1298749/100‑1 (MQ=255)
                                                                          aGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGAtc                          <  2:425968/100‑1 (MQ=255)
                                                                          aGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGAtc                          >  1:1178015/1‑100 (MQ=255)
                                                                               taacgatgttCCCGCTGTGGGTCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCag                      >  2:142082/9‑100 (MQ=37)
                                                                                cccTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCagga                    >  1:938865/1‑100 (MQ=255)
                                                                                cccTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCagga                    >  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)
                                                                                                   |                                                                                                   
GGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACG  >  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: