Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,370,026 T→A 15.6% intergenic (+23/‑190) pspE → / → ycjM thiosulfate:cyanide sulfurtransferase (rhodanese)/alpha amylase catalytic domain family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,370,0260TA15.6% 394.0 / 62.0 160intergenic (+23/‑190)pspE/ycjMthiosulfate:cyanide sulfurtransferase (rhodanese)/alpha amylase catalytic domain family protein
Reads supporting (aligned to +/- strand):  ref base T (60/75);  new base A (6/19);  total (66/94)
Fisher's exact test for biased strand distribution p-value = 7.62e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.19e-01

Only 100 of 162 total aligned reads displayed.

GAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCC  >  NC_000913/1369932‑1370125
                                                                                              |                                                                                                    
gagaTGGGATATGCCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATgc                                                                                                 >  1:2591228/1‑100 (MQ=255)
 agaTGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATgcg                                                                                                >  2:282691/1‑100 (MQ=255)
  gaTGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATgcgc                                                                                               >  1:964733/1‑100 (MQ=255)
  gaTGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCATAAgcgc                                                                                               <  2:1099645/100‑1 (MQ=255)
     gggATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCATAAGCGCCgc                                                                                            <  1:3248898/100‑2 (MQ=255)
      ggATATATCCACGTTGAGAAAGCCGTTGGTTTGGCAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCATAAGCGCCgct                                                                                           <  2:1407054/100‑3 (MQ=255)
      ggATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCATAAGCGCCgct                                                                                           <  1:2808018/100‑3 (MQ=255)
      ggATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCATAAGCGCCgct                                                                                           <  2:1888277/100‑3 (MQ=255)
       gATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGttt                                                                                         <  1:3551888/100‑1 (MQ=255)
        atatACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCATAAGCGCCgctca                                                                                         <  1:2621388/100‑5 (MQ=255)
         tataCCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGttttt                                                                                       <  1:1868062/100‑1 (MQ=255)
         tataCCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCATAAGCGCCgctcat                                                                                        <  1:559032/100‑6 (MQ=255)
         tataCCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCATAAGCGCCgctcat                                                                                        <  1:387934/100‑6 (MQ=255)
          attcccACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCATAAGCGCCgctcatg                                                                                       <  1:1894401/97‑7 (MQ=255)
          ataCCCACGTTTAGAACGCCGGTGGCTTGAACGACATCCCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAACGGCTCATAAGCGCCgctcat                                                                                        <  2:2765835/99‑6 (MQ=255)
           taCCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGttttttt                                                                                     >  1:198239/1‑100 (MQ=255)
           taCCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGttttttt                                                                                     <  1:1589739/100‑1 (MQ=255)
           taCCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCATAAGCGCCgctcatgg                                                                                      <  2:1123595/100‑8 (MQ=255)
            aCCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGtttttttt                                                                                    >  2:452067/1‑100 (MQ=255)
            aCCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCATAAGCGCCgctcatggt                                                                                     <  1:3576634/100‑9 (MQ=255)
            aCCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCATAAGCGCCgctcatggt                                                                                     <  1:2477799/100‑9 (MQ=255)
            aCCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCATAAGCGCCgctcatggt                                                                                     <  1:1689405/100‑9 (MQ=255)
            aCCCACGTTGAGAACGCCGGTGGCCTGAAACACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCATAAGCGCCgctcatgg                                                                                      <  2:411984/99‑8 (MQ=255)
             cccACGTTGAGAACGCCGTTGGCCTGAAAGATATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCATAAGCGCCgctcatgg                                                                                      <  2:905070/98‑8 (MQ=255)
             cccACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGttttttttt                                                                                   >  2:3248696/1‑100 (MQ=255)
             cccACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTAACCATGAGCGGCGATTATGCGCCGttttttttt                                                                                   >  1:1507570/1‑100 (MQ=255)
              ccACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTc                                                                                  <  2:723162/100‑1 (MQ=255)
               cACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCt                                                                                 <  1:2943357/100‑1 (MQ=255)
               cACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCt                                                                                 <  1:2056223/100‑1 (MQ=255)
                aCGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGttttttttttct                                                                                >  2:1469525/1‑97 (MQ=255)
                   ttGAGAACGGCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTTCTGTc                                                                              <  2:2331642/100‑1 (MQ=255)
                   ttGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTTCTGTc                                                                              <  1:1446348/100‑1 (MQ=255)
                   ttGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTca                                                                             >  2:1937863/1‑100 (MQ=255)
                    tGAGAACGCGGGTGGCATGAAAGACATCGCAAGGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTGATGCGCCGTTTTTTTTTCTGTca                                                                             <  2:2485996/99‑1 (MQ=255)
                    tGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTcac                                                                            >  2:87742/1‑100 (MQ=255)
                    tGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTca                                                                             >  2:1231930/1‑99 (MQ=255)
                     gagaACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTcaca                                                                           <  2:3702170/100‑1 (MQ=255)
                       gaACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACAcc                                                                         <  1:667357/100‑1 (MQ=255)
                         aCGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCtt                                                                       >  1:538318/1‑100 (MQ=255)
                         aCGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCtt                                                                       >  2:3154355/1‑100 (MQ=255)
                          cGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTa                                                                      >  1:3426864/1‑100 (MQ=255)
                             cGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTAttt                                                                   <  1:200230/100‑1 (MQ=255)
                                tGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACa                                                                <  1:2259185/100‑1 (MQ=255)
                                tGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACa                                                                <  2:1741008/100‑1 (MQ=255)
                                  gacTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATc                                                              <  2:59565/98‑1 (MQ=255)
                                    cTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTTCTGTCACACCTTATTTACATcc                                                             <  2:1630182/100‑1 (MQ=255)
                                      gAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCa                                                          >  1:2717158/1‑100 (MQ=255)
                                       aaaGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCAt                                                         <  2:2127522/100‑1 (MQ=255)
                                         aGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTTTCACACCTTATTTACATCCCCATAg                                                       >  2:2680263/1‑100 (MQ=255)
                                           aCATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGat                                                     >  1:1023273/1‑100 (MQ=255)
                                           aCATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGat                                                     <  1:3154376/100‑1 (MQ=255)
                                            cATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTTTTTTTTTTCTTATTTACATCCCTCTAGat                                                     >  2:466473/1‑99 (MQ=255)
                                            cATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGatt                                                    >  1:289050/1‑100 (MQ=255)
                                            cATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGatt                                                    >  1:691012/1‑100 (MQ=255)
                                             aTCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCTTAGatt                                                    >  2:567093/1‑99 (MQ=255)
                                              tCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTTTTTTTCAACTTATTTACATCCGCATAGa                                                      >  2:2346192/1‑96 (MQ=255)
                                                gCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTAttt                                                >  1:1228869/1‑100 (MQ=255)
                                                gCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTAttt                                                <  2:1809552/100‑1 (MQ=255)
                                                 cAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTTTTTTTTTCCTTATTTACATCCCCATAGATTAttt                                                >  2:2401136/1‑99 (MQ=255)
                                                   aTGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCg                                             <  1:3726893/100‑1 (MQ=255)
                                                    tGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGt                                            <  2:1885855/100‑1 (MQ=255)
                                                      ccGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCa                                          <  2:523580/100‑1 (MQ=255)
                                                          aGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTc                                      <  2:213815/100‑1 (MQ=255)
                                                           ggTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCt                                       >  2:822170/1‑98 (MQ=255)
                                                            gTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCACATAGATTATTTGCGTCAGCTcac                                    >  2:334089/1‑100 (MQ=255)
                                                             tCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTcaca                                   <  1:3790728/100‑1 (MQ=255)
                                                             tCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTcaca                                   >  1:100972/1‑100 (MQ=255)
                                                                aaGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAAt                                >  1:2413625/1‑100 (MQ=255)
                                                                 aGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCTTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGcaca                                     <  2:2316170/94‑1 (MQ=255)
                                                                  ggTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTTCTGTTTCACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAAt                                >  1:1790100/1‑99 (MQ=255)
                                                                      aaaaGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCtt                          >  2:3568077/1‑100 (MQ=255)
                                                                       aaaGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCttt                         <  2:2211670/100‑1 (MQ=255)
                                                                       aaaGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCttt                         <  2:1988243/100‑1 (MQ=255)
                                                                       aaaGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCttt                         >  1:260072/1‑100 (MQ=255)
                                                                         aGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCttttt                       <  2:1275413/100‑1 (MQ=255)
                                                                          gAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCGGCTCACAAATACGCTTTTTc                      <  2:2413613/100‑1 (MQ=255)
                                                                           aaTTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTcc                     <  1:1261420/100‑1 (MQ=255)
                                                                           aaTTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTcc                     <  1:2227532/100‑1 (MQ=255)
                                                                            agtCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTccc                    <  2:2669676/98‑1 (MQ=255)
                                                                            aTTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTTTCCCACCATCTTTACATCCCCCTAGAATATTTGCGTCAGCTCACAAATACGCCTTTTccc                    >  2:2276268/1‑100 (MQ=255)
                                                                             ttCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCt                   <  2:926495/100‑1 (MQ=255)
                                                                             ttCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCt                   <  1:2953188/100‑1 (MQ=255)
                                                                              tcgccATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTg                  <  1:1748031/97‑1 (MQ=255)
                                                                              tCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCGGTAACACCTTATTTACATCCTCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTg                  <  2:1231860/100‑1 (MQ=255)
                                                                               cACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTg                  <  1:2530817/100‑1 (MQ=255)
                                                                               cACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTgg                 <  1:87742/100‑1 (MQ=255)
                                                                                aaaaaaaacggcgCATAAGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGt                >  1:451593/9‑100 (MQ=255)
                                                                                 ccATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTa               <  2:677019/100‑1 (MQ=255)
                                                                                 aaaaaaacggcgCATAAGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTa               >  1:686984/8‑100 (MQ=255)
                                                                                   aaaaacggcgCATAAGCGCCGTTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTaa              >  1:667297/6‑100 (MQ=255)
                                                                                    tGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTaaaa            >  1:2609623/1‑100 (MQ=255)
                                                                                       gcggcgCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAAt         <  1:872049/100‑1 (MQ=255)
                                                                                          gcgcTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAAt         <  2:3409164/97‑1 (MQ=255)
                                                                                           cgcTTATGCGGCGTTTTTTTTTCGTTCGCACCGTATTTACATCCCCATAGATTATTTGCGGCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGa       <  2:3554918/98‑1 (MQ=255)
                                                                                           cgcTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGAtt     >  2:3746789/1‑100 (MQ=255)
                                                                                            gcTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGAttt    >  1:773425/1‑100 (MQ=255)
                                                                                            gcTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGAttt    <  1:245665/100‑1 (MQ=255)
                                                                                            gcTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTACCAGGGaaaaagcgtatt    <  2:2268410/100‑8 (MQ=255)
                                                                                            gcTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACAGCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGAttt    <  1:3290281/100‑1 (MQ=255)
                                                                                              ttATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTcc  <  2:198239/100‑1 (MQ=255)
                                                                                              |                                                                                                    
GAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCC  >  NC_000913/1369932‑1370125

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: