Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 940,754 G→A 15.9% intergenic (+34/‑205) 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,7540GA15.9% 265.8 / 36.4 113intergenic (+34/‑205)serS/dmsAseryl‑tRNA synthetase/dimethyl sulfoxide reductase, anaerobic, subunit A
Reads supporting (aligned to +/- strand):  ref base G (51/43);  new base A (6/12);  total (57/56)
Fisher's exact test for biased strand distribution p-value = 1.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.82e-01

Only 100 of 118 total aligned reads displayed.

GATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTT  >  NC_000913/940655‑940851
                                                                                                   |                                                                                                 
gATGGTGGTATTGAAGTACCAGAAGGTCTGCGTCTGTATATCAACGAAATGGAAGATATTGGCTAATACCCAATTTTTCTGAATCTAAAAGGCGCCTGCg                                                                                                   <  1:11852/100‑1 (MQ=255)
gATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCg                                                                                                   <  2:1124630/100‑1 (MQ=255)
gATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCg                                                                                                   <  2:285900/100‑1 (MQ=255)
  tGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCggg                                                                                                 >  2:1547531/1‑100 (MQ=255)
  tGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCggg                                                                                                 >  1:1728051/1‑100 (MQ=255)
     tCGTATTGACGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGgcgc                                                                                              <  2:1132683/100‑1 (MQ=255)
          ttGAAGTACCATAAGTTCTGCGTCCGTTTATGTTCGGACTGTAATATAGTGGCTAATACCCAATTTTTCTGAATCTAATAAGCGCCCGCAGGCGCttttt                                                                                         <  1:2664387/100‑1 (MQ=255)
           tGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCCGCAGGCGCttttta                                                                                        <  2:916632/100‑2 (MQ=255)
           tGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCCGCAGGCGCttttta                                                                                        <  2:2204770/100‑2 (MQ=255)
           tGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCCGCAGGCGCttttta                                                                                        <  2:2838472/100‑2 (MQ=255)
            gAAGTACCAGAAGTTCTTCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGATTCTATGAATCGCCTGCGGTCTCttttttt                                                                                       <  1:142321/100‑1 (MQ=255)
            gAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCttttttt                                                                                       >  2:923623/1‑100 (MQ=255)
            gAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCttttttt                                                                                       >  1:2169043/1‑100 (MQ=255)
            gAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCCGCAGGCGCtttttag                                                                                       <  2:2039827/100‑3 (MQ=255)
             aaGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTg                                                                                      >  2:351937/1‑100 (MQ=255)
             aaGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTg                                                                                      <  1:510923/100‑1 (MQ=255)
              aGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTg                                                                                      <  1:2492825/99‑1 (MQ=255)
              aGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGt                                                                                     <  2:2728720/100‑1 (MQ=255)
               gTACCAGAAGTTCTGCGTCCGTATATGCACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGtc                                                                                    <  1:2613250/100‑1 (MQ=255)
               gTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCCGCAGGCGCtttttagatt                                                                                    <  1:2074617/100‑6 (MQ=38)
                tACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCCGCAGGCGCtttttagattc                                                                                   <  2:1327082/100‑7 (MQ=38)
                 aCCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGtctc                                                                                  <  2:2459519/100‑1 (MQ=255)
                 aCCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGtctc                                                                                  >  2:2563650/1‑100 (MQ=255)
                 aCCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGtctc                                                                                  >  2:2342131/1‑100 (MQ=255)
                 aCCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGtctc                                                                                  >  2:528412/1‑100 (MQ=255)
                 aCCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTATTTTTGtctc                                                                                  >  2:1698473/1‑100 (MQ=255)
                  ccAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTcc                                                                                 >  2:2493402/1‑100 (MQ=255)
                    aGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTTTTAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCt                                                                               >  1:1773215/1‑100 (MQ=255)
                    aGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCt                                                                               >  2:2313201/1‑100 (MQ=255)
                    aGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCCGCAGGCGCtttttagattcagaa                                                                               <  1:1052465/100‑11 (MQ=37)
                     gAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCtt                                                                              >  2:2302774/1‑100 (MQ=255)
                     gAAGTTCTGCGTCCGTATATGAACGGACTCGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCCCCTGCGGGCGCTTTTTTTGTCTCCCtt                                                                              <  1:1819605/100‑1 (MQ=255)
                      aaGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCttt                                                                             <  2:208355/100‑1 (MQ=255)
                      aaGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCttt                                                                             <  2:769705/100‑1 (MQ=255)
                            tGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATAcc                                                                       >  1:2605968/1‑100 (MQ=255)
                             gCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCg                                                                      >  2:2648371/1‑100 (MQ=255)
                              cGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGa                                                                     <  2:1106644/100‑1 (MQ=255)
                                tCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAAc                                                                   <  1:1346968/100‑1 (MQ=255)
                                      tatGAACGGACTGGAATATATTGGCTAATCCCCAATTTTTCTGACTATAACAAGCGCCTGCGGGCGCTGTTTTTGTCTCTCTTTTATACCg                                                                      >  1:757783/1‑91 (MQ=255)
                                      tatGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACaataat                                                             <  1:406753/100‑1 (MQ=255)
                                      tatGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACaataat                                                             <  2:1691341/100‑1 (MQ=255)
                                       atGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAAtt                                                            <  1:511927/100‑1 (MQ=255)
                                       atGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAAtt                                                            >  2:2321415/1‑100 (MQ=255)
                                        tGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTa                                                           >  1:2503735/1‑100 (MQ=255)
                                        tGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTa                                                           >  2:2847733/1‑100 (MQ=255)
                                         gAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTAc                                                          <  2:308712/100‑1 (MQ=255)
                                         gAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTAc                                                          <  2:1398786/100‑1 (MQ=255)
                                         gAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTGTCTCCCTTTGATACCGAACAATAATTAct                                                         >  1:75972/1‑100 (MQ=255)
                                            cGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTActcc                                                       <  2:1079899/100‑1 (MQ=255)
                                             ggACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTActcct                                                      >  1:333696/1‑100 (MQ=255)
                                             ggACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTActcc                                                       <  1:1365137/99‑1 (MQ=255)
                                              gACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTActcctc                                                     <  2:114979/100‑1 (MQ=255)
                                               aCTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTActcct                                                      <  1:1848881/98‑1 (MQ=255)
                                                cTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCAc                                                   <  2:1820676/100‑1 (MQ=255)
                                                  ggAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACtt                                                 >  2:1321243/1‑100 (MQ=255)
                                                  ggAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACtt                                                 >  2:1297618/1‑100 (MQ=255)
                                                  ggAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACtt                                                 <  2:1685235/100‑1 (MQ=255)
                                                  ggAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACtt                                                 >  1:991245/1‑100 (MQ=255)
                                                  ggAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACtt                                                 >  2:2092147/1‑100 (MQ=255)
                                                  ggAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACtt                                                 >  2:884494/1‑100 (MQ=255)
                                                     atatatTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTaca                                              <  1:261691/100‑1 (MQ=255)
                                                       atatTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACg                                            <  2:2169056/100‑1 (MQ=255)
                                                          ttGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTaa                                         <  2:865211/100‑1 (MQ=255)
                                                           tGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAAt                                        <  2:2791542/100‑1 (MQ=255)
                                                            ggCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAAta                                       >  2:2816200/1‑100 (MQ=255)
                                                              cTAATACCCAATTTTTCTGAATCTAAAATGCGCCTCCGGGCGCTTTTTTTGTCGCCCTTTGATACCTAACACTAATTACTCATCACTTACACGTAAtac                                      >  1:1642104/1‑99 (MQ=255)
                                                               tAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAAtacta                                    >  2:483814/1‑100 (MQ=255)
                                                                    cccAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTc                               >  2:1247173/1‑100 (MQ=255)
                                                                    cccAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTc                               >  2:2366465/1‑100 (MQ=255)
                                                                    cccAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTc                               >  1:1559852/1‑100 (MQ=255)
                                                                     ccAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCg                              <  2:1784959/100‑1 (MQ=255)
                                                                      cgaTTTTTCTGAAACTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTTATACCGAACAATAATTACTCCTCACTTACACGTAGTACTACTTTcta                             >  2:1977835/3‑98 (MQ=255)
                                                                      cAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGCACCATTGTTACTCCACACTTACACCAAACACTACTCTCg                              >  1:221424/1‑99 (MQ=255)
                                                                      cAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTACCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGa                             >  1:1613975/1‑100 (MQ=255)
                                                                      cAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAg                            >  2:912173/1‑100 (MQ=255)
                                                                       aaTTTTTCTGAATCTCATATGCGCCGGCGGGCGTTTTCTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAg                            <  1:815273/100‑1 (MQ=255)
                                                                       aaTTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAg                            <  1:694777/100‑1 (MQ=255)
                                                                          ttttCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGa                         <  2:1736645/100‑1 (MQ=255)
                                                                            ttCTGAATCTAAACAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCCAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGaaa                       <  2:274920/100‑1 (MQ=255)
                                                                             tCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGaaaa                      >  2:291711/1‑100 (MQ=255)
                                                                             tCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGaaaa                      <  1:1506689/100‑1 (MQ=255)
                                                                               gagacaaaaaaaGCGCCCGCAGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATc                    >  1:1717582/8‑100 (MQ=38)
                                                                                gAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATc                    <  1:1281039/99‑1 (MQ=255)
                                                                                  acaaaaaaaGCGCCCGCAGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTAc                 >  1:1767571/5‑100 (MQ=39)
                                                                                  aTCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGCGTGAAAATCTAc                 <  2:1917547/100‑1 (MQ=255)
                                                                                  aTCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAAGCTAc                 >  2:1250537/1‑100 (MQ=255)
                                                                                   tCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTTCTTTCGAGTGAAAATCTAcc                <  1:1787459/100‑1 (MQ=255)
                                                                                   tCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTAcc                <  1:2383846/100‑1 (MQ=255)
                                                                                   caaaaaaaGCGCCCGCAGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTAcc                >  1:1381432/4‑100 (MQ=255)
                                                                                    aaaaaaaGCGCCCGCAGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCt               >  1:1179054/3‑100 (MQ=255)
                                                                                    aaaaaaaGCGCCCGCAGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCt               >  2:1786490/3‑100 (MQ=255)
                                                                                     tAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATCATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTa              >  1:286449/1‑100 (MQ=255)
                                                                                     tAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTACTACTACTTTTGAGTGAAAATATACCTa              >  1:1304589/1‑100 (MQ=255)
                                                                                     aaaaaaGCGCCCGCAGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTa              >  2:2202333/2‑100 (MQ=255)
                                                                                      aaaaaGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTAt             >  1:2057989/1‑100 (MQ=255)
                                                                                          aGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATctct         >  2:162590/1‑100 (MQ=255)
                                                                                            cgcCTGCGGGCGCTTTTTTTTTTTTTTTTTAATACCAAACAATAATTACGCCTCACTTACACTTGATATTACTTTCGATTGAAAATCTACCTATCTCtt        >  1:1697929/1‑99 (MQ=255)
                                                                                              ccTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGa     <  1:804569/100‑1 (MQ=255)
                                                                                               cTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGAt    <  1:1321238/100‑1 (MQ=255)
                                                                                                 gCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGaaga  >  1:1440363/1‑97 (MQ=255)
                                                                                                   |                                                                                                 
GATGGTCGTATTGAAGTACCAGAAGTTCTGCGTCCGTATATGAACGGACTGGAATATATTGGCTAATACCCAATTTTTCTGAATCTAAAAAGCGCCTGCGGGCGCTTTTTTTGTCTCCCTTTGATACCGAACAATAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTT  >  NC_000913/940655‑940851

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: