Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,009,731 A→G 15.5% intergenic (+19/‑90) 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,7310AG15.5% 276.6 / 43.1 123intergenic (+19/‑90)yedF/yedKputative TusA family sulfurtransferase/DUF159 family protein
Reads supporting (aligned to +/- strand):  ref base A (51/53);  new base G (6/13);  total (57/66)
Fisher's exact test for biased strand distribution p-value = 2.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.52e-01

Only 100 of 123 total aligned reads displayed.

GGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGA  >  NC_000913/2009634‑2009829
                                                                                                  |                                                                                                  
ggATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGAcc                                                                                                  >  1:672855/1‑100 (MQ=255)
ggATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGAcc                                                                                                  >  1:1745949/1‑100 (MQ=255)
ggATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGAcc                                                                                                  >  1:3915447/1‑100 (MQ=255)
ggATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACAGTCCTACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGAcc                                                                                                  <  2:619777/100‑1 (MQ=255)
  aTGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCca                                                                                                <  1:3110137/100‑1 (MQ=255)
   tGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCcac                                                                                               <  1:451111/100‑1 (MQ=255)
     cgcgTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAg                                                                                             >  2:3247936/1‑100 (MQ=255)
     cgcgTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAg                                                                                             >  1:3299169/1‑100 (MQ=255)
       cgTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCg                                                                                           <  1:1273315/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)
         tAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGGGTCACAGCggg                                                                                         <  1:733433/100‑1 (MQ=255)
         tAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGGGTCACAGCggg                                                                                         <  2:1354151/100‑1 (MQ=255)
         tAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGGGTCACAGCggg                                                                                         <  1:1123869/100‑1 (MQ=255)
         tAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCggg                                                                                         >  1:3094325/1‑100 (MQ=255)
          aaTCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGGGTCACAGCGGGa                                                                                        <  1:1577062/100‑1 (MQ=255)
          aaTCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGGGTCACAGCGGGa                                                                                        <  2:653320/100‑1 (MQ=255)
          aaTCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGGGTCACAGCGGGa                                                                                        <  1:2106238/100‑1 (MQ=255)
          aaTCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCggg                                                                                         >  2:324169/1‑99 (MQ=255)
          aaTCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGa                                                                                        >  1:900246/1‑100 (MQ=255)
          aaTCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGa                                                                                        <  2:3463850/100‑1 (MQ=255)
           aTCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGGGTCACAGCggg                                                                                         <  2:1327194/98‑1 (MQ=255)
           aTCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGGGTCACAGCGGGaa                                                                                       <  2:1327044/100‑1 (MQ=255)
            tCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAAc                                                                                      <  1:796765/100‑1 (MQ=255)
             cACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAAc                                                                                      <  2:3807794/99‑1 (MQ=255)
             cACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACa                                                                                     >  2:796761/1‑100 (MQ=255)
                ggCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCg                                                                                  <  2:389789/100‑1 (MQ=255)
                 gCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGt                                                                                 <  1:1343376/100‑1 (MQ=255)
                 gCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGt                                                                                 <  1:2861705/100‑1 (MQ=255)
                  cTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGtt                                                                                <  1:2152285/100‑1 (MQ=255)
                   tataCGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTa                                                                               <  1:3530695/100‑1 (MQ=255)
                    ataCGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTAc                                                                              >  1:1197248/1‑100 (MQ=255)
                      gtcttGCTGAATATCCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCt                                                                            >  2:1048976/5‑100 (MQ=255)
                       cGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCtt                                                                           >  2:539725/1‑100 (MQ=255)
                       cGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCATCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCtt                                                                           >  1:952802/1‑100 (MQ=255)
                         gTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTAACttt                                                                          >  2:3636105/1‑99 (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:2442726/100‑1 (MQ=255)
                               gATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATcac                                                                   <  1:1877683/100‑1 (MQ=255)
                               gATAGTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATcac                                                                   <  2:1110290/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)
                                    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)
                                         aaGACGCGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTa                                                         <  2:3148066/100‑1 (MQ=255)
                                          aGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTAc                                                        <  1:1188744/100‑1 (MQ=255)
                                             cGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCg                                                     <  1:2383468/100‑1 (MQ=255)
                                               ggCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGcc                                                   <  1:3925098/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)
                                                            ttATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTc                                      <  2:514230/100‑1 (MQ=255)
                                                            ttATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTc                                      >  1:2184063/1‑100 (MQ=255)
                                                             tATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCt                                     <  2:3473225/100‑1 (MQ=255)
                                                               tttAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCttt                                   <  1:581068/100‑1 (MQ=255)
                                                                ttAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCtttt                                  >  1:3329704/1‑100 (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)
                                                                       aGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATg                           <  2:2860925/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)
                                                                                cccTCATTCGCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCagg                   <  2:1197239/99‑1 (MQ=255)
                                                                                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)
                                                                                   cTCATTCCAGCTGTGACACACAGCGGGAACATCGCTACCCTTGTATCACATAGATTTTACCCGCCCCCACTCCTACTCTTTTCAGTATGAtctc                      >  2:200401/1‑94 (MQ=255)
                                                                                       ttCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTc            <  1:1955895/100‑1 (MQ=255)
                                                                                        tCCCGCTGTGGGTCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCt           >  2:3063121/1‑100 (MQ=255)
                                                                                        tCCCGCTGTGACCCACAGCGGGAACATCGTTCCCTTTATATCACATAGATTTTACGCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCt           <  2:2654111/100‑1 (MQ=255)
                                                                                        tCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTc            <  1:2551977/99‑1 (MQ=255)
                                                                                        tCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCt           >  2:2450049/1‑100 (MQ=255)
                                                                                        tCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCt           >  1:1760550/1‑100 (MQ=255)
                                                                                        tCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCt           >  1:1760280/1‑100 (MQ=255)
                                                                                        tCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCt           <  2:672851/100‑1 (MQ=255)
                                                                                         cccGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTa          >  2:107069/1‑100 (MQ=255)
                                                                                         cccGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTa          >  1:638422/1‑100 (MQ=255)
                                                                                          ccGCTGTGGGTCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTAt         >  1:233115/1‑100 (MQ=255)
                                                                                          acGCTGTGACCCACAGCGGCAACATCGTTACCTTTGTATGACATAGATTTTACCCGCGCCCACTCCTATGCTGTTCAGTATGATCTCaggagggtact           >  2:90160/2‑93 (MQ=255)
                                                                                                tgACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTgg   <  2:338274/100‑1 (MQ=255)
                                                                                                tgACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCATATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTgg   <  2:1828871/100‑1 (MQ=255)
                                                                                                 gACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTAtagatctg  >  1:2693797/1‑93 (MQ=255)
                                                                                                 gACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGa  <  2:832112/100‑1 (MQ=255)
                                                                                                  |                                                                                                  
GGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGA  >  NC_000913/2009634‑2009829

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: