Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 940,751 T→C 15.8% 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.8% 389.5 / 64.0 172intergenic (+31/‑208)serS/dmsAseryl‑tRNA synthetase/dimethyl sulfoxide reductase, anaerobic, subunit A
Reads supporting (aligned to +/- strand):  ref base T (88/56);  new base C (9/18);  total (97/75)
Fisher's exact test for biased strand distribution p-value = 1.05e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.58e-01

Only 100 of 173 total aligned reads displayed.

TGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATT  >  NC_000913/940657‑940850
                                                                                              |                                                                                                   
tGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCggg                                                                                                >  1:2599483/1‑100 (MQ=255)
 ggTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGgc                                                                                               <  1:1433066/100‑1 (MQ=255)
  gTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGgcg                                                                                              >  1:4056232/1‑100 (MQ=255)
     gTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCCGCAGGCGCtt                                                                                           <  2:2315742/100‑1 (MQ=255)
      tATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCttt                                                                                          >  1:1738347/1‑100 (MQ=255)
       aTTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCCGCAGGCGCtttt                                                                                         <  1:813619/100‑1 (MQ=255)
        ttGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCttttt                                                                                        >  2:1135336/1‑100 (MQ=255)
        ttGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCttttt                                                                                        >  1:3904581/1‑100 (MQ=255)
         tGAAGTACCAGAAGTTCTGCTTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCCGCAGGCGCttttta                                                                                       <  2:2797452/100‑2 (MQ=255)
         tGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAATCTCCTGCGGGCGCtttttt                                                                                       >  2:3505808/1‑100 (MQ=255)
         tGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCtttttt                                                                                       >  2:1762205/1‑100 (MQ=255)
         tGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCtttttt                                                                                       <  1:2465823/100‑1 (MQ=255)
         tGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCCGCAGGCGCttttta                                                                                       <  2:138213/100‑2 (MQ=255)
         tGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCCGCAGGCGCttttta                                                                                       <  1:3306223/100‑2 (MQ=255)
         tGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCCGCAGGCGCttttta                                                                                       <  1:1179075/100‑2 (MQ=255)
          gAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCCGCAGGCGCtttttag                                                                                      <  2:651605/100‑3 (MQ=255)
          gAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCCGCAGGCGCtttttag                                                                                      <  2:4155001/100‑3 (MQ=255)
           aaGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTg                                                                                     <  1:2300797/100‑1 (MQ=255)
           aaGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTg                                                                                     <  1:3745624/100‑1 (MQ=255)
           aaGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTg                                                                                     <  1:2346735/100‑1 (MQ=255)
           aaGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATACATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCCGCAGGCGCtttttaga                                                                                     <  1:1013760/100‑4 (MQ=38)
              tACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCCGCAGGCGCtttttagattc                                                                                  <  1:314743/100‑7 (MQ=38)
              tACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCCGCAGGCGCtttttagattc                                                                                  <  1:1129346/100‑7 (MQ=38)
               aCCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGtctc                                                                                 >  2:2725626/1‑100 (MQ=255)
               aCCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGtctc                                                                                 <  1:1723864/100‑1 (MQ=255)
               aCCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGtctc                                                                                 >  1:1489962/1‑100 (MQ=255)
               aCCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGtctc                                                                                 >  1:1321266/1‑100 (MQ=255)
                 cAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCCGCAGGCGCtttttagattcaga                                                                               <  2:1266935/100‑10 (MQ=37)
                 cAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCCGCAGGCGCtttttagattcag                                                                                <  1:1231315/99‑9 (MQ=37)
                    aaGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCttn                                                                            <  2:156295/100‑2 (MQ=255)
                    aaGTTCTGCGTCCGTATATGAACGGACTGAAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCttt                                                                            <  2:1959229/100‑1 (MQ=255)
                     aGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTg                                                                           <  1:880346/100‑1 (MQ=255)
                      gTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGa                                                                          <  2:3042575/100‑1 (MQ=255)
                       ttCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGAt                                                                         >  2:3003078/1‑100 (MQ=255)
                            cGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGa                                                                    >  1:3353445/1‑100 (MQ=255)
                            cGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGa                                                                    >  1:3487153/1‑100 (MQ=255)
                              tCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAAc                                                                  <  1:706495/100‑1 (MQ=255)
                                 gTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACaat                                                               <  2:1711509/100‑1 (MQ=255)
                                 gTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTCTGATACCGAACaat                                                               <  1:3072199/100‑1 (MQ=255)
                                   atatGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACaataa                                                             >  2:2914090/1‑100 (MQ=255)
                                   atatGAACGGACTGGAATATATTGGCTAATACCCAAGTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACaataa                                                             >  2:830219/1‑100 (MQ=255)
                                      tGAACGGACTGGAATATATTGGCTAATACCCAATTTCTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAAGTa                                                          >  2:1993795/1‑100 (MQ=255)
                                         aCGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTActc                                                       >  1:3820166/1‑100 (MQ=255)
                                          cGGACTGGAATATAATGGCTACTACCCAATTTTTCGGGACCCAAAAAGCGCCTGCGGGCGCTTTCTGTGGCTCCCCTAAATCCCGAGCAATAATTActcc                                                      <  2:4144265/100‑1 (MQ=255)
                                          cGGACTAGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTActcc                                                      <  2:3723399/100‑1 (MQ=255)
                                            gACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTActcctc                                                    <  2:1598850/100‑1 (MQ=255)
                                                ggAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACtt                                                >  1:993179/1‑100 (MQ=255)
                                                   atataaTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTACGGGCGCTTTTTTTGTCGCCCTTTGATACCGAACAATAATTAATCCTCTCTTaca                                             <  2:3554649/100‑1 (MQ=255)
                                                    tatatTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTacac                                            <  1:568266/100‑1 (MQ=255)
                                                     atatTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACg                                           <  1:2063493/100‑1 (MQ=255)
                                                     atatTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACg                                           <  2:2354467/100‑1 (MQ=255)
                                                      tatTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGt                                          >  1:973959/1‑100 (MQ=255)
                                                         tGGCTAATACCCAATTTTTCTGAATCTACAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTa                                         >  2:4005790/1‑98 (MQ=255)
                                                         tGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAAt                                       <  1:2934655/100‑1 (MQ=255)
                                                          ggCTAATCCCCCATTTTTCTGAATCTAAAAAGCGCCTGCGGCCGCCTTTTTTGTCTCCCTTTGATGCAGAACAATAATTATTCCTCACTTACACGTAAta                                      >  2:2206836/1‑100 (MQ=255)
                                                           gCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAAtac                                     <  1:2818974/100‑1 (MQ=255)
                                                            cTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAAtact                                    >  1:643706/1‑100 (MQ=255)
                                                              aaTACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTTTTTCCCTTTGATACCGAACAATAATTACTCCTTTCTTACACGTAAtact                                    >  2:349547/1‑98 (MQ=255)
                                                              aaTACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAAtactan                                  <  1:3691596/100‑2 (MQ=255)
                                                              aaTACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAAtactac                                  >  1:2111104/1‑100 (MQ=255)
                                                              aaTACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAAtactac                                  >  1:2242764/1‑100 (MQ=255)
                                                                  cccAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACttt                               >  1:2958527/1‑99 (MQ=255)
                                                                   ccAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTAACTTACACGTAATACTACTTTCg                             >  1:1344976/1‑100 (MQ=255)
                                                                    cAATTTTTCTGAATCTAAAGAGCGCCGGCGGGCGCTTTTTTTGTCTCCCTTGGATACCGAACAATAATTACGCCTCACTTACACGTAATACTACTTTCGa                            <  2:2982235/100‑1 (MQ=255)
                                                                     aaTTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAg                           <  2:2047052/100‑1 (MQ=255)
                                                                     aaTTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAg                           <  2:2877714/100‑1 (MQ=255)
                                                                     aaTTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAg                           <  2:3862741/100‑1 (MQ=255)
                                                                     aaTTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAg                           <  1:3045036/100‑1 (MQ=255)
                                                                     aaTTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAg                           <  1:1311193/100‑1 (MQ=255)
                                                                     aaTTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAg                           >  1:1264359/1‑100 (MQ=255)
                                                                      aTTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGt                          <  2:544414/100‑1 (MQ=255)
                                                                      aTTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGt                          <  1:3332800/100‑1 (MQ=255)
                                                                       tttttCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTg                         <  2:223285/100‑1 (MQ=255)
                                                                           tCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGaaaa                     >  2:3242042/1‑100 (MQ=255)
                                                                           tCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGaaaa                     >  1:3569618/1‑100 (MQ=255)
                                                                            cTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAATGaaaa                     >  2:3358606/1‑99 (MQ=255)
                                                                             tGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATc                   >  1:442246/1‑100 (MQ=255)
                                                                              gAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCt                  >  1:2569137/1‑100 (MQ=255)
                                                                              gAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCt                  >  2:3677/1‑100 (MQ=255)
                                                                              gAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCt                  >  1:2871803/1‑100 (MQ=255)
                                                                               gacaaaaaaaGCGCCCGCAGGCGCTTTTTTTGTCTCCCTTTGATACCGAACACTAATTACTCCTCACTTACACGTAATACTACTTTCGAGTTAAAATCTa                 >  2:1034396/6‑100 (MQ=38)
                                                                               actCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTa                 <  2:2533242/98‑1 (MQ=255)
                                                                               aaTCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTa                 <  2:106261/100‑1 (MQ=255)
                                                                                aTCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTAc                <  2:3258340/100‑1 (MQ=255)
                                                                                aTCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTAc                >  2:1819730/1‑100 (MQ=255)
                                                                                 tCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTAcc               <  1:1736888/100‑1 (MQ=255)
                                                                                 tCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTAcc               <  2:2131474/100‑1 (MQ=255)
                                                                                 caaaaaaaGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTAcc               >  1:1032565/4‑100 (MQ=255)
                                                                                 caaaaaaaGCGCCCGCAGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAACTc                   >  2:553157/4‑96 (MQ=38)
                                                                                 caaaaaaaGCGCCCGCAGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTAcc               >  1:104332/4‑100 (MQ=39)
                                                                                  cTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCt              >  1:2666135/1‑100 (MQ=255)
                                                                                   aaaaaaGCGCCCGCAGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTa             >  2:144736/2‑100 (MQ=255)
                                                                                   aaaaaaGCGCCCGCAGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTa             >  2:2599985/2‑100 (MQ=255)
                                                                                    aaaaaGCGCCCGCAGGCGCTTTTTTTGTCTCCCTTTAAGACCGAACAATCCTTCCTGCTCCCTTCCACCTCATACTACTCTCGAGTGAACATCTccc               >  2:2589154/1‑97 (MQ=255)
                                                                                       aaGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTAtctc         <  1:1625228/100‑1 (MQ=255)
                                                                                        aGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATctct        >  1:3414871/1‑100 (MQ=255)
                                                                                        aGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATctct        >  2:816711/1‑100 (MQ=255)
                                                                                          cgcCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCttt      >  2:3298050/1‑100 (MQ=255)
                                                                                          cgcCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCttt      >  2:711485/1‑100 (MQ=255)
                                                                                            ccTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGa    >  2:3146255/1‑100 (MQ=255)
                                                                                              |                                                                                                   
TGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATT  >  NC_000913/940657‑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: