Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,009,733 C→T 16.7% intergenic (+21/‑88) 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,7330CT16.7% 277.2 / 41.4 119intergenic (+21/‑88)yedF/yedKputative TusA family sulfurtransferase/DUF159 family protein
Reads supporting (aligned to +/- strand):  ref base C (47/51);  new base T (6/14);  total (54/65)
Fisher's exact test for biased strand distribution p-value = 2.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.78e-01

Only 100 of 125 total aligned reads displayed.

GGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGC  >  NC_000913/2009634‑2009832
                                                                                                    |                                                                                                   
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)
       cgTAATCACGGCTATACGGTGCTTGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCg                                                                                              <  1:1574717/100‑1 (MQ=255)
       cgTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCg                                                                                              <  1:1273315/100‑1 (MQ=255)
       cgTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCg                                                                                              >  1:2159274/1‑100 (MQ=255)
         tAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGGGTCACAGCggg                                                                                            <  1:733433/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)
         tAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGGGTCACAGCggg                                                                                            <  1:1244741/100‑1 (MQ=255)
         tAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGGGTCACAGCggg                                                                                            <  1:1123869/100‑1 (MQ=255)
         tAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCggg                                                                                            >  1:3094325/1‑100 (MQ=255)
          aaTCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGGGTCACAGCGGGa                                                                                           <  2:653320/100‑1 (MQ=255)
          aaTCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGGGTCACAGCGGGa                                                                                           <  1:1577062/100‑1 (MQ=255)
          aaTCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGa                                                                                           >  1:900246/1‑100 (MQ=255)
           aTCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGGGTCACAGCggg                                                                                            <  2:1327194/98‑1 (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:1343376/100‑1 (MQ=255)
                  cTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCGCTCATTCCCGCTGTGACCCACAGCGGGAACATCGtt                                                                                   >  2:2227025/1‑100 (MQ=255)
                  cTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGtt                                                                                   <  1:2152285/100‑1 (MQ=255)
                   tataCGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTa                                                                                  <  1:3530695/100‑1 (MQ=255)
                   tataCGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTa                                                                                  >  2:3235665/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)
                       cGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCATCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCtt                                                                              >  1:952802/1‑100 (MQ=255)
                         gTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTAACttt                                                                             >  2:3636105/1‑99 (MQ=255)
                           gCTGGATATTCAGCAAGACGGGCCGACCATTCGATATTTAATTAAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTa                                                                          >  2:1686502/1‑100 (MQ=255)
                             tGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTa                                                                          <  2:3697557/98‑1 (MQ=255)
                               gATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATcac                                                                      <  1:1877683/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)
                                       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)
                                             cGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCg                                                        <  1:2383468/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)
                                                       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                                         <  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)
                                                                       aGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCTCTCCTATTCTTTTCAGTATg                              <  2:1868695/100‑1 (MQ=255)
                                                                       aGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTAt                               >  1:2465044/1‑99 (MQ=255)
                                                                       aGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATg                              <  2:2860925/100‑1 (MQ=255)
                                                                        gAAGTAATCCCGCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTGCCCGCCCCCACTCCTATTCTTTTCAGTATGa                             <  2:1071184/100‑1 (MQ=255)
                                                                         aaGTAATCCCTCATTCCCGCTGTGACCCACAGCTGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGAt                            <  2:1298749/100‑1 (MQ=255)
                                                                          aGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGAtc                           >  1:1178015/1‑100 (MQ=255)
                                                                          aGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGAtc                           <  2:425968/100‑1 (MQ=255)
                                                                               taacgatgttCCCGCTGTGGGTCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCag                       >  2:142082/9‑100 (MQ=37)
                                                                                cccTCATTCGCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCagg                      <  2:1197239/99‑1 (MQ=255)
                                                                                cccTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCagga                     >  1:938865/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              >  1:1760550/1‑100 (MQ=255)
                                                                                        tCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCt              >  2:2450049/1‑100 (MQ=255)
                                                                                        tCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCt              <  2:672851/100‑1 (MQ=255)
                                                                                         cccGCTGTGGGTCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTa             >  2:3739152/1‑100 (MQ=255)
                                                                                         cccGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTa             >  1:638422/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)
                                                                                                tgACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCATATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTgg      <  2:1828871/100‑1 (MQ=255)
                                                                                                 gACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTAtagatctg     >  1:2693797/1‑93 (MQ=255)
                                                                                                 gACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGa     <  2:832112/100‑1 (MQ=255)
                                                                                                   cccACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACg   <  2:1078713/100‑1 (MQ=255)
                                                                                                    ccACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGc  >  1:2160092/1‑100 (MQ=255)
                                                                                                    |                                                                                                   
GGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGC  >  NC_000913/2009634‑2009832

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: