Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,370,029 T→A 16.0% intergenic (+26/‑187) 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,0290TA16.0% 389.7 / 63.5 162intergenic (+26/‑187)pspE/ycjMthiosulfate:cyanide sulfurtransferase (rhodanese)/alpha amylase catalytic domain family protein
Reads supporting (aligned to +/- strand):  ref base T (60/76);  new base A (6/20);  total (66/96)
Fisher's exact test for biased strand distribution p-value = 5.17e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.17e-01

Only 100 of 163 total aligned reads displayed.

GAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTG  >  NC_000913/1369932‑1370127
                                                                                                 |                                                                                                   
tagaTGGGATATACCCTCGTTGAGACCGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCTGCGCTTAt                                                                                                     <  2:1403917/97‑1 (MQ=255)
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)
      ggATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCATAAGCGCCgct                                                                                             <  1:697652/100‑3 (MQ=255)
       gATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGttt                                                                                           <  1:3551888/100‑1 (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)
          ataCCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCATAAGCGCCgctcatg                                                                                         <  1:2986671/100‑7 (MQ=255)
           taCCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGttttttt                                                                                       >  1:198239/1‑100 (MQ=255)
           taCCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGttttttt                                                                                       <  1:238953/100‑1 (MQ=255)
           taCCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGttttttt                                                                                       <  2:181128/100‑1 (MQ=255)
           taCCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGttttttt                                                                                       <  1:1589739/100‑1 (MQ=255)
            aCCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGtttttttt                                                                                      >  2:452067/1‑100 (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)
             cccACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGttttttttt                                                                                     >  2:1976938/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)
                 cGTTGAGAACTCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGt                                                                                 <  2:2591173/100‑1 (MQ=255)
                   ttGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTTCTGTc                                                                                <  1:1446348/100‑1 (MQ=255)
                   ttGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTca                                                                               >  2:1937863/1‑100 (MQ=255)
                   ttGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTca                                                                               <  2:2325012/100‑1 (MQ=255)
                    tGAGAACGCGGGTGGCATGAAAGACATCGCAAGGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTGATGCGCCGTTTTTTTTTCTGTca                                                                               <  2:2485996/99‑1 (MQ=255)
                    tGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTcac                                                                              >  2:87742/1‑100 (MQ=255)
                     gagaACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTcaca                                                                             <  2:3702170/100‑1 (MQ=255)
                         aCGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCtt                                                                         >  2:3154355/1‑100 (MQ=255)
                          cGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTa                                                                        >  1:3426864/1‑100 (MQ=255)
                             cGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTAttt                                                                     <  1:200230/100‑1 (MQ=255)
                              ggTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCCCACCTTATTTa                                                                    >  1:551985/1‑100 (MQ=255)
                               gTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTAc                                                                   >  2:2397603/1‑100 (MQ=255)
                                tGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACa                                                                  <  2:1741008/100‑1 (MQ=255)
                                  gacTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATc                                                                <  2:59565/98‑1 (MQ=255)
                                  gCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGAGCCGTTTTTTTTTCTGTTACACCTTATTTa                                                                    >  1:941743/1‑96 (MQ=255)
                                   ccTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATcc                                                               >  1:178425/1‑100 (MQ=255)
                                      gAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCa                                                            >  1:2717158/1‑100 (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)
                                              tCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTTTTTTTCAACTTATTTACATCCGCATAGa                                                        >  2:2346192/1‑96 (MQ=255)
                                                gCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTAttt                                                  >  1:1228869/1‑100 (MQ=255)
                                                gCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTAttt                                                  <  2:1809552/100‑1 (MQ=255)
                                                   aTGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCg                                               <  1:3726893/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)
                                                            gTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTcac                                      >  1:339077/1‑100 (MQ=255)
                                                            gTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCACATAGATTATTTGCGTCAGCTcac                                      >  2:334089/1‑100 (MQ=255)
                                                             tCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTcaca                                     >  1:100972/1‑100 (MQ=255)
                                                              cAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACaa                                    <  2:1608077/100‑1 (MQ=255)
                                                               aaaGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACaaa                                   <  1:3636513/100‑1 (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)
                                                                      aaaaGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCtt                            >  1:677023/1‑100 (MQ=255)
                                                                       aaaGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCttt                           >  1:260072/1‑100 (MQ=255)
                                                                       aaaGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCttt                           <  2:2211670/100‑1 (MQ=255)
                                                                       aaaGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCttt                           <  2:1988243/100‑1 (MQ=255)
                                                                          gAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCGGCTCACAAATACGCTTTTTc                        <  2:2413613/100‑1 (MQ=255)
                                                                           aaTTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTcc                       <  1:2227532/100‑1 (MQ=255)
                                                                           aaTTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTcc                       <  1:1261420/100‑1 (MQ=255)
                                                                            aTTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTTTCCCACCATCTTTACATCCCCCTAGAATATTTGCGTCAGCTCACAAATACGCCTTTTccc                      >  2:2276268/1‑100 (MQ=255)
                                                                             ttCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCt                     <  2:926495/100‑1 (MQ=255)
                                                                              tCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCGGTAACACCTTATTTACATCCTCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTg                    <  2:1231860/100‑1 (MQ=255)
                                                                              tCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTgg                   >  2:1261536/1‑100 (MQ=255)
                                                                               cACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTg                    <  1:2530817/100‑1 (MQ=255)
                                                                               cACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTgg                   <  2:345727/100‑1 (MQ=255)
                                                                               cACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTgg                   <  2:208174/100‑1 (MQ=255)
                                                                                aaaaaaaacggcgCATAAGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGt                  >  1:451593/9‑100 (MQ=255)
                                                                                 ccATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTa                 <  2:677019/100‑1 (MQ=255)
                                                                                 aaaaaaacggcgCATAAGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTa                 >  1:2915528/8‑100 (MQ=255)
                                                                                 aaaaaaacggcgCATAAGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTa                 >  1:686984/8‑100 (MQ=255)
                                                                                 aaaaaaacggcgCATAAGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTa                 >  1:3817859/8‑100 (MQ=255)
                                                                                   aaaaacggcgCATAAGCGCCGTTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTaa                >  1:667297/6‑100 (MQ=255)
                                                                                    tGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTaaaa              >  1:2609623/1‑100 (MQ=255)
                                                                                     gAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATATATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTaaaaa             <  1:2694540/100‑1 (MQ=255)
                                                                                       gcggcgCTTATGCGCCGTTTTTTTTTTTTTTTTTTCTTATTAACATCGCCATAGATTATTTGATTTAGCTCGGAAATACGCTTTTTCCCTTGTtaaaa             >  2:527271/1‑98 (MQ=255)
                                                                                       gcggcgCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAAt           <  1:872049/100‑1 (MQ=255)
                                                                                        cggcgCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATg          >  2:2212662/1‑100 (MQ=255)
                                                                                          gcgcTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAAt           <  2:3409164/97‑1 (MQ=255)
                                                                                           cgcTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGAtt       >  2:3746789/1‑100 (MQ=255)
                                                                                           cgcTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGAtt       <  2:1183917/100‑1 (MQ=255)
                                                                                            gcTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGAttt      <  1:245665/100‑1 (MQ=255)
                                                                                            gcTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGAttt      >  1:773425/1‑100 (MQ=255)
                                                                                                aTGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTg  <  2:1606625/100‑1 (MQ=255)
                                                                                                 |                                                                                                   
GAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTG  >  NC_000913/1369932‑1370127

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: