Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 940,751 T→C 15.7% intergenic (+31/‑208) serS → / → dmsA seryl‑tRNA synthetase/dimethyl sulfoxide reductase, anaerobic, subunit A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913940,7510TC15.7% 399.8 / 61.3 171intergenic (+31/‑208)serS/dmsAseryl‑tRNA synthetase/dimethyl sulfoxide reductase, anaerobic, subunit A
Reads supporting (aligned to +/- strand):  ref base T (76/68);  new base C (7/20);  total (83/88)
Fisher's exact test for biased strand distribution p-value = 1.19e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

Only 100 of 173 total aligned reads displayed.

CTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATT  >  NC_000913/940653‑940850
                                                                                                  |                                                                                                   
cTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTg                                                                                                    >  2:3012290/1‑100 (MQ=255)
  gATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCg                                                                                                  >  1:1003859/1‑100 (MQ=255)
  gATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCCGca                                                                                                  <  1:2948753/100‑2 (MQ=255)
   aTGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCgg                                                                                                 >  2:330620/1‑100 (MQ=255)
      gTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGgcg                                                                                              >  2:1800398/1‑100 (MQ=255)
      gTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGgcg                                                                                              <  2:2525132/100‑1 (MQ=255)
      gTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGgcg                                                                                              >  2:1100562/1‑100 (MQ=255)
      gTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGgcg                                                                                              >  2:3248329/1‑100 (MQ=255)
        cGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCt                                                                                            <  1:3053596/100‑1 (MQ=255)
        cGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCt                                                                                            <  1:2737918/100‑1 (MQ=255)
          tATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCCGCAGGCGCttt                                                                                          <  2:1957009/100‑1 (MQ=255)
            ttGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCttttt                                                                                        >  1:987139/1‑100 (MQ=255)
             tGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCtttttt                                                                                       >  2:2745954/1‑100 (MQ=255)
             tGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCCGCAGGCGCttttta                                                                                       <  2:3000675/100‑2 (MQ=255)
             tGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCCGCAGGCGCttttta                                                                                       <  1:1962375/100‑2 (MQ=255)
             tGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCCGCAGGCGCttttta                                                                                       <  1:26373/100‑2 (MQ=255)
             tGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCGAAAAAGCGCCCGCAGGCGCttttta                                                                                       <  1:1240819/100‑2 (MQ=255)
              gAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCCGCAGGCGCtttttag                                                                                      <  2:1045746/100‑3 (MQ=255)
               aaGTGCCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCCGCAGGCGCtttttag                                                                                      <  2:151264/99‑3 (MQ=39)
               aaGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTg                                                                                     <  1:3389898/100‑1 (MQ=255)
               aaGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCCGCAGGCGCttttt                                                                                        <  1:2599321/97‑1 (MQ=255)
               aaGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAAGATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCCGCAGGCGCtttttaga                                                                                     <  1:1195447/100‑4 (MQ=38)
                aGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCCGCAGGCGCtttttagat                                                                                    <  2:1162956/100‑5 (MQ=39)
                 ggaCCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTg                                                                                     <  2:1974960/96‑1 (MQ=255)
                 gTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGtc                                                                                   <  2:1997221/100‑1 (MQ=255)
                  tACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCCGCAGGCGCtttttagattc                                                                                  <  1:843012/100‑7 (MQ=38)
                   aCCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGtctc                                                                                 >  2:1410436/1‑100 (MQ=255)
                   aCCAGAAGTTCTGCGTCCGTATATGAACGAACTAGAATATATTGGCTAATACCCAATTTTTCTGAATCAAAAAAGCGCCCGCAGGCGCtttttag                                                                                      <  2:906913/95‑3 (MQ=38)
                    ccAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTcc                                                                                >  2:2398193/1‑100 (MQ=255)
                    ccAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTcc                                                                                >  2:2884119/1‑100 (MQ=255)
                     cAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTccc                                                                               <  1:1275903/100‑1 (MQ=255)
                      aGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCt                                                                              >  2:515148/1‑100 (MQ=255)
                      aGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCt                                                                              >  2:964893/1‑100 (MQ=255)
                       gAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCtt                                                                             <  2:334739/100‑1 (MQ=255)
                        aaGTTCTGCGTCCGCATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCttt                                                                            <  2:2483692/100‑1 (MQ=255)
                          gTTCTGCGTTCGTATATGAACGGACTGGAATATATCGGCCAATACCCGATTTTTTTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTTCCttt                                                                            <  2:473461/98‑1 (MQ=255)
                           ttCTGCAGCCGTAGATGAACTTACTGCAATATATTGGCTAATACCCAATTTTTCTGTGTCTAAAAAGCGCCTACGGGCGCTTTTTTTGTATTCCTTTCAt                                                                         >  2:3406585/1‑100 (MQ=255)
                            tCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATa                                                                        >  2:3468642/1‑100 (MQ=255)
                            tCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATa                                                                        <  2:2333110/100‑1 (MQ=255)
                             cTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATAc                                                                       >  1:3547992/1‑100 (MQ=255)
                              tGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATAcc                                                                      <  1:681373/100‑1 (MQ=255)
                                cGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCg                                                                     >  1:5769/1‑99 (MQ=255)
                                cGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGa                                                                    <  1:2787114/100‑1 (MQ=255)
                                  tCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAAc                                                                  <  2:1673750/100‑1 (MQ=255)
                                  tCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAAc                                                                  <  2:3596613/100‑1 (MQ=255)
                                  tCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAAc                                                                  <  1:2673274/100‑1 (MQ=255)
                                  tCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAAc                                                                  <  1:2398209/100‑1 (MQ=255)
                                    cGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACaa                                                                >  1:1388493/1‑100 (MQ=255)
                                      tataTGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACaata                                                              <  1:3286879/100‑1 (MQ=255)
                                         atGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAAtt                                                           <  1:1330620/100‑1 (MQ=255)
                                          tGAACGGACTGGAATATATTGGCTAATACCCAATTTTTTTTATTTTTTAAAGCGCCTGAGGGCACTTTTTTTGTTTCTGTTTGATACCGAACAATAATTa                                                          >  2:1089751/1‑100 (MQ=255)
                                          tGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTa                                                          >  1:2779273/1‑100 (MQ=255)
                                          tGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTa                                                          >  2:2085256/1‑100 (MQ=255)
                                             aCGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTActc                                                       >  1:690867/1‑100 (MQ=255)
                                             aCGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTActc                                                       >  1:2900579/1‑100 (MQ=255)
                                             aCGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTActc                                                       <  1:608712/100‑1 (MQ=255)
                                               ggACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTActcct                                                     >  1:2715986/1‑100 (MQ=255)
                                                gACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTTTTTTTTTTTTTTACCCAAAAATAATTActcct                                                     >  2:772749/1‑99 (MQ=255)
                                                gACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTActcctc                                                    >  1:550677/1‑100 (MQ=255)
                                                  cTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCAc                                                  <  1:2176485/100‑1 (MQ=255)
                                                   tGGAATGTATTGGCTACTACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGAGCGCTTTCTTTGTCTCCCTTTGAGACCGCACAATAATTACTCCTCACt                                                 <  2:1627975/100‑1 (MQ=255)
                                                    ggAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTCCTCCTCACtt                                                >  1:671864/1‑100 (MQ=255)
                                                    ggAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACtt                                                <  1:926287/100‑1 (MQ=255)
                                                    ggAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACtt                                                >  1:1696918/1‑100 (MQ=255)
                                                    ggAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACtt                                                >  2:3118258/1‑100 (MQ=255)
                                                     agaaaaatTGGGTATTAGCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTa                                               >  2:2374269/7‑100 (MQ=255)
                                                       atatatTGGCTAATACCCAATTTTTCTGAATCCAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTaca                                             >  1:3454214/1‑100 (MQ=255)
                                                           atTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTa                                         >  2:2073388/1‑100 (MQ=255)
                                                             tGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAAt                                       >  1:3498431/1‑100 (MQ=255)
                                                             tGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAAt                                       <  1:1009128/100‑1 (MQ=255)
                                                              ggCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAAta                                      >  1:3359518/1‑100 (MQ=255)
                                                               gCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAAtac                                     >  1:3128445/1‑100 (MQ=255)
                                                                  aaTACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAAtactac                                  <  1:1621659/100‑1 (MQ=255)
                                                                      cccAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTc                              >  1:1718568/1‑100 (MQ=255)
                                                                       ccAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCg                             <  1:611274/100‑1 (MQ=255)
                                                                        cAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGa                            <  2:748293/100‑1 (MQ=255)
                                                                                cTGAATCTAAAAAGCGGCTGCGGGCACTTTTTTTGTCTCCCTTTGATACCGCACAATAATTCCTACTCAATTAAACGTTATACTACTTTGTAGTGaaa                      >  2:385615/1‑98 (MQ=255)
                                                                                cTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAAt                    <  1:2229319/100‑1 (MQ=255)
                                                                                cTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAAt                    >  2:3464251/1‑100 (MQ=255)
                                                                                 tGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATc                   >  1:2400346/1‑100 (MQ=255)
                                                                                  gAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCt                  >  1:156953/1‑100 (MQ=255)
                                                                                     tCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTAcc               <  2:1431417/100‑1 (MQ=255)
                                                                                     tCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTAcc               <  2:1003852/100‑1 (MQ=255)
                                                                                     caaaaaaaGCGCCCGCAGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGaaaaagatcg               >  2:659790/4‑94 (MQ=38)
                                                                                     caaaaaaaGCGCCCGCAGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTAcc               >  1:3310043/4‑100 (MQ=39)
                                                                                     caaaaaaaGCGCCCGCAGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTAcc               >  1:299755/4‑100 (MQ=39)
                                                                                      cTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAAAAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGACAATCTACCt              >  2:2315634/1‑100 (MQ=255)
                                                                                      cTAAAAAGCGCCTGCGGGCGCCTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGaaagatcggaa              >  1:2322425/1‑92 (MQ=255)
                                                                                      aaaaaaaGCGCCCGCAGGCGCTTTTTTTGTCTCTCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCt              >  2:1852336/3‑100 (MQ=39)
                                                                                       tAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACTTAATACTACTTTCAAGTGAAAATCTACCTa             >  1:537550/1‑100 (MQ=255)
                                                                                       tAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTa             >  1:1394637/1‑100 (MQ=255)
                                                                                       aaaaaaGCGCCCGCAGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTa             >  2:1263048/2‑100 (MQ=255)
                                                                                        aaaaaGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTAt            >  1:1767773/1‑100 (MQ=255)
                                                                                             gcgcCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCtt       >  1:3579717/1‑100 (MQ=255)
                                                                                              cgcCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCttt      >  2:474106/1‑100 (MQ=255)
                                                                                              cgcCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCttt      >  1:48650/1‑100 (MQ=255)
                                                                                               gcCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTg     <  1:1001842/100‑1 (MQ=255)
                                                                                                ccTGCGGGCGCTTTTTTTGTTTCACTTTGATACCGAACAATAATTACTCCTCCCTTACATGTAATACTACTTTCGAGTGAAAATCTAACTATATCttt      >  2:1327235/1‑98 (MQ=255)
                                                                                                 cTGCGGGCGCTTTTTTGGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGAt   <  2:2125842/100‑1 (MQ=255)
                                                                                                  tGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGAtt  >  2:2063422/1‑100 (MQ=255)
                                                                                                  |                                                                                                   
CTGATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATT  >  NC_000913/940653‑940850

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: