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.

GAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTG  >  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)
       gATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑ttt                                                                                           <  1:3551888/100‑1 (MQ=255)
         tataCCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑ttttt                                                                                         <  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)
           taCCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑ttttttt                                                                                       >  1:198239/1‑100 (MQ=255)
           taCCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑ttttttt                                                                                       <  1:238953/100‑1 (MQ=255)
           taCCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑ttttttt                                                                                       <  2:181128/100‑1 (MQ=255)
           taCCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑ttttttt                                                                                       <  1:1589739/100‑1 (MQ=255)
            aCCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑tttttttt                                                                                      >  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)
             cccACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑ttttttttt                                                                                     >  2:3248696/1‑100 (MQ=255)
             cccACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑ttttttttt                                                                                     >  2:1976938/1‑100 (MQ=255)
             cccACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTAACCATGAGCGGCGATTATGCGCCG‑ttttttttt                                                                                     >  1:1507570/1‑100 (MQ=255)
              ccACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTc                                                                                    <  2:723162/100‑1 (MQ=255)
               cACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCt                                                                                   <  1:2943357/100‑1 (MQ=255)
               cACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCt                                                                                   <  1:2056223/100‑1 (MQ=255)
                 cGTTGAGAACTCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGt                                                                                 <  2:2591173/100‑1 (MQ=255)
                   ttGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTTCTGTc                                                                                <  1:1446348/100‑1 (MQ=255)
                   ttGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTca                                                                               >  2:1937863/1‑100 (MQ=255)
                   ttGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTca                                                                               <  2:2325012/100‑1 (MQ=255)
                    tGAGAACGCGGGTGGCATGAAAGACATCGCAAGGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTGATGCGCCG‑TTTTTTTTTCTGTca                                                                               <  2:2485996/99‑1 (MQ=255)
                    tGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTcac                                                                              >  2:87742/1‑100 (MQ=255)
                     gagaACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTcaca                                                                             <  2:3702170/100‑1 (MQ=255)
                         aCGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCtt                                                                         >  2:3154355/1‑100 (MQ=255)
                          cGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTa                                                                        >  1:3426864/1‑100 (MQ=255)
                             cGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTAttt                                                                     <  1:200230/100‑1 (MQ=255)
                              ggTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCCCACCTTATTTa                                                                    >  1:551985/1‑100 (MQ=255)
                               gTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTAc                                                                   >  2:2397603/1‑100 (MQ=255)
                                tGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACa                                                                  <  2:1741008/100‑1 (MQ=255)
                                  gacTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATc                                                                <  2:59565/98‑1 (MQ=255)
                                  gCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGAGCCG‑TTTTTTTTTCTGTTACACCTTATTTa                                                                    >  1:941743/1‑96 (MQ=255)
                                   ccTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATcc                                                               >  1:178425/1‑100 (MQ=255)
                                      gAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCa                                                            >  1:2717158/1‑100 (MQ=255)
                                         aGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTTTCACACCTTATTTACATCCCCATAg                                                         >  2:2680263/1‑100 (MQ=255)
                                           aCATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGat                                                       >  1:1023273/1‑100 (MQ=255)
                                           aCATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGat                                                       <  1:3154376/100‑1 (MQ=255)
                                            cATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTTTTTTTTTTCTTATTTACATCCCTCTAGat                                                       >  2:466473/1‑99 (MQ=255)
                                              tCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTTTTTTTCAACTTATTTACATCCGCATAGa                                                        >  2:2346192/1‑96 (MQ=255)
                                                gCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTAttt                                                  >  1:1228869/1‑100 (MQ=255)
                                                gCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTAttt                                                  <  2:1809552/100‑1 (MQ=255)
                                                   aTGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCg                                               <  1:3726893/100‑1 (MQ=255)
                                                      ccGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCa                                            <  2:523580/100‑1 (MQ=255)
                                                          aGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTc                                        <  2:213815/100‑1 (MQ=255)
                                                           ggTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCt                                         >  2:822170/1‑98 (MQ=255)
                                                            gTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTcac                                      >  1:339077/1‑100 (MQ=255)
                                                            gTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCACATAGATTATTTGCGTCAGCTcac                                      >  2:334089/1‑100 (MQ=255)
                                                             tCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTcaca                                     >  1:100972/1‑100 (MQ=255)
                                                              cAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACaa                                    <  2:1608077/100‑1 (MQ=255)
                                                               aaaGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACaaa                                   <  1:3636513/100‑1 (MQ=255)
                                                                aaGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAAt                                  >  1:2413625/1‑100 (MQ=255)
                                                                 aGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCT‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGcaca                                       <  2:2316170/94‑1 (MQ=255)
                                                                  ggTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTTCTGTTTCACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAAt                                  >  1:1790100/1‑99 (MQ=255)
                                                                      aaaaGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCtt                            >  2:3568077/1‑100 (MQ=255)
                                                                      aaaaGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCtt                            >  1:677023/1‑100 (MQ=255)
                                                                       aaaGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCttt                           >  1:260072/1‑100 (MQ=255)
                                                                       aaaGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCttt                           <  2:2211670/100‑1 (MQ=255)
                                                                       aaaGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCttt                           <  2:1988243/100‑1 (MQ=255)
                                                                          gAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCGGCTCACAAATACGCTTTTTc                        <  2:2413613/100‑1 (MQ=255)
                                                                           aaTTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTcc                       <  1:2227532/100‑1 (MQ=255)
                                                                           aaTTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTcc                       <  1:1261420/100‑1 (MQ=255)
                                                                            aTTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTTTCCCACCATCTTTACATCCCCCTAGAATATTTGCGTCAGCTCACAAATACGCCTTTTccc                      >  2:2276268/1‑100 (MQ=255)
                                                                             ttCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCt                     <  2:926495/100‑1 (MQ=255)
                                                                              tCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCGGTAACACCTTATTTACATCCTCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTg                    <  2:1231860/100‑1 (MQ=255)
                                                                              tCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTgg                   >  2:1261536/1‑100 (MQ=255)
                                                                               cACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTg                    <  1:2530817/100‑1 (MQ=255)
                                                                               cACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTgg                   <  2:345727/100‑1 (MQ=255)
                                                                               cACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTgg                   <  2:208174/100‑1 (MQ=255)
                                                                                aaaaaaaacggcgCATAAGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGt                  >  1:451593/9‑100 (MQ=255)
                                                                                 ccATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTa                 <  2:677019/100‑1 (MQ=255)
                                                                                 aaaaaaacggcgCATAAGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTa                 >  1:2915528/8‑100 (MQ=255)
                                                                                 aaaaaaacggcgCATAAGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTa                 >  1:686984/8‑100 (MQ=255)
                                                                                 aaaaaaacggcgCATAAGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTa                 >  1:3817859/8‑100 (MQ=255)
                                                                                   aaaaacggcgCATAAGCGCCGTTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTaa                >  1:667297/6‑100 (MQ=255)
                                                                                    tGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTaaaa              >  1:2609623/1‑100 (MQ=255)
                                                                                     gAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATATATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTaaaaa             <  1:2694540/100‑1 (MQ=255)
                                                                                       gcggcgCTTATGCGCCG‑TTTTTTTTTTTTTTTTTTCTTATTAACATCGCCATAGATTATTTGATTTAGCTCGGAAATACGCTTTTTCCCTTGTtaaaa             >  2:527271/1‑98 (MQ=255)
                                                                                       gcggcgCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAAt           <  1:872049/100‑1 (MQ=255)
                                                                                        cggcgCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATg          >  2:2212662/1‑100 (MQ=255)
                                                                                          gcgcTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAAt           <  2:3409164/97‑1 (MQ=255)
                                                                                           cgcTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGAtt       >  2:3746789/1‑100 (MQ=255)
                                                                                           cgcTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGAtt       <  2:1183917/100‑1 (MQ=255)
                                                                                            gcTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGAttt      <  1:245665/100‑1 (MQ=255)
                                                                                            gcTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGAttt      >  1:773425/1‑100 (MQ=255)
                                                                                                aTGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTg  <  2:1606625/100‑1 (MQ=255)
                                                                                                 |                                                                                                   
GAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTG  >  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: ‑