Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,370,044 (T)9→8 100% intergenic (+41/‑172) 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,0360T.59.0% 30.3 / 165.7 83intergenic (+33/‑180)pspE/ycjMthiosulfate:cyanide sulfurtransferase (rhodanese)/alpha amylase catalytic domain family protein
Reads supporting (aligned to +/- strand):  ref base T (16/18);  new base . (28/21);  total (44/39)
Fisher's exact test for biased strand distribution p-value = 3.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.50e-01

GGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTG  >  NC_000913/1369937‑1370131
                                                                                                    |                                                                                               
gggATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑t                                                                                                 >  1:423285/1‑100 (MQ=255)
  gATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑ttt                                                                                               >  1:739268/1‑100 (MQ=255)
    tataCCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑ttttt                                                                                             >  2:501492/1‑100 (MQ=255)
      taCCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑ttttttt                                                                                           <  1:1013788/100‑1 (MQ=255)
      taCCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑ttttttt                                                                                           >  1:1002106/1‑100 (MQ=255)
      taCCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑ttttttt                                                                                           >  1:1847984/1‑100 (MQ=255)
      taCCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑ttttttt                                                                                           >  2:32702/1‑100 (MQ=255)
      taCCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑ttttttt                                                                                           >  2:1631472/1‑100 (MQ=255)
        cccACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑ttttttttc                                                                                         >  1:2062310/1‑99 (MQ=255)
         ccACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑ttttttttct                                                                                        <  1:851138/100‑3 (MQ=255)
          cACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑ttttttttctg                                                                                       >  2:2130728/1‑97 (MQ=255)
          cACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTc                                                                                        <  1:2171198/99‑1 (MQ=255)
           aCGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTg                                                                                      <  1:439244/100‑1 (MQ=255)
            cGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTCTGTc                                                                                    <  2:2188583/100‑1 (MQ=255)
             gTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTc                                                                                    <  2:423286/100‑1 (MQ=255)
             gTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTCTGTca                                                                                   >  1:403776/1‑100 (MQ=255)
             gTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTCTGTca                                                                                   >  1:2149643/1‑100 (MQ=255)
               tGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTcac                                                                                  <  2:932311/100‑1 (MQ=255)
               tGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTcac                                                                                  <  1:1006381/100‑1 (MQ=255)
                 agaACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCacac                                                                                <  2:558272/100‑1 (MQ=255)
                 agaACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTCTGTCACAcc                                                                               <  1:1095192/100‑1 (MQ=255)
                    aCGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTTTGTCTCTCCtt                                                                             >  2:969079/1‑100 (MQ=255)
                    aCGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCtt                                                                             >  1:720136/1‑100 (MQ=255)
                    aCGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTCTGTCACACCTTa                                                                            >  1:482876/1‑100 (MQ=255)
                       ccgatctccTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTCTGTCACACCTTAttt                                                                         <  1:2196643/93‑1 (MQ=255)
                        cGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTAttt                                                                         <  2:568710/100‑1 (MQ=255)
                            ggCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTCTGTCACACCTTATTTACATc                                                                    >  2:1492036/1‑100 (MQ=255)
                             gCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTCTGTCACACCTTATTTACATcc                                                                   >  1:1203890/1‑100 (MQ=255)
                              ccTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTCTGTCACACCTTATTTAGATcgg                                                                  >  2:2196628/1‑98 (MQ=255)
                               cttaaaTACATTGCAATGCATAAGGTCAAAGTTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATccc                                                                  <  1:7856/97‑1 (MQ=255)
                                tGAAATACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTCTGTCACACCTTATTTACATCCCCa                                                                <  2:2062296/100‑1 (MQ=255)
                                tGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTTCTGTCACACCTTATTTACATccc                                                                  <  1:1941051/100‑1 (MQ=255)
                                tGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTCTGTCACACCTTATTTACATCCCCa                                                                >  1:1361824/1‑100 (MQ=255)
                                tGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTCTGTCACACCTTATTTACATCCCCa                                                                >  1:1031527/1‑100 (MQ=255)
                                 gAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCa                                                                >  1:523432/1‑100 (MQ=255)
                                 gAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTCTGTCACACCTTATTTACATCCCCAt                                                               >  2:1594206/1‑100 (MQ=255)
                                 gAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTCTGTCACACCTTATTTACATCCCCAt                                                               <  2:403777/100‑1 (MQ=255)
                                   aaGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTCTGTCACACCTTATTTACATCACCAAAg                                                             >  1:809065/1‑100 (MQ=255)
                                   aaGACATCGCAATACCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAg                                                             <  1:590348/100‑1 (MQ=255)
                                       cATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACGATGAGCGGCGCTTATGCGCCGTTTTTTTTTTTTTTTATACCTTATCTACATCCCCATAGa                                                            >  2:1375140/1‑99 (MQ=255)
                                       cATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGatta                                                         >  1:1718919/1‑100 (MQ=255)
                                       cATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGatta                                                         >  1:343340/1‑100 (MQ=255)
                                       cATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGatta                                                         >  1:1161026/1‑100 (MQ=255)
                                        aTCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGatta                                                         >  1:302318/1‑100 (MQ=255)
                                         tCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGattatt                                                       >  1:608534/1‑100 (MQ=255)
                                          cGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTAttt                                                      >  1:1855460/1‑100 (MQ=255)
                                           gCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTAttt                                                      >  2:1975048/1‑100 (MQ=255)
                                           ccgatctcGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTg                                                     <  1:1114771/93‑1 (MQ=255)
                                            cAATGTCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGc                                                    <  1:654559/100‑1 (MQ=255)
                                             aaTGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCg                                                   <  2:1161022/100‑1 (MQ=255)
                                                 ccGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAg                                               >  2:1090663/1‑100 (MQ=255)
                                                 acggaGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTCTGTCACACCTTATTTACATCCCTATAGATTATTTgtgt                                                  >  2:1012057/5‑97 (MQ=255)
                                                  cGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAg                                               >  1:165202/1‑100 (MQ=255)
                                                  cGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTgagatcgg                                              >  2:1114768/1‑93 (MQ=255)
                                                    aaGTTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTc                                            <  1:1792349/100‑1 (MQ=255)
                                                     aGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTca                                           <  2:1031525/100‑1 (MQ=255)
                                                      ggTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTcac                                          >  2:651366/1‑100 (MQ=255)
                                                       gTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTcac                                          <  1:1331906/100‑1 (MQ=255)
                                                         cAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACaa                                        >  2:1963327/1‑100 (MQ=255)
                                                           aaGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTCTGTCACACCTTATTTCCATCCCCATAGATTATTTGCGTCAGCTCACAAATa                                     >  2:23951/1‑100 (MQ=255)
                                                           aaGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATa                                     <  2:851389/100‑1 (MQ=255)
                                                            agtttAAAAGGATTCATCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGACACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATa                                     <  1:937592/97‑1 (MQ=255)
                                                            aGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGCTTATTTGCGTCCGCTCACAAATAc                                    >  2:957086/1‑100 (MQ=255)
                                                               ttAAAAGAATTGACCATGAGCGGCGCTGATGCGCCG‑‑TTTTTTTTCTGTCACTCCTTATTTACCTCCCCATAGATTATTTGCGTCAGCTCACAAATACGCt                                 <  1:1492046/100‑1 (MQ=255)
                                                               ttAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGc                                  <  2:224807/100‑1 (MQ=255)
                                                                tAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCtt                                <  1:1142016/100‑1 (MQ=255)
                                                                tAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCtt                                >  2:1859186/1‑100 (MQ=255)
                                                                  aaaGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCt                                 >  1:585041/1‑98 (MQ=255)
                                                                  aaaGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTT‑TTTTTTTTTTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCtttt                              >  2:1053757/1‑100 (MQ=255)
                                                                      aaTTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTT‑TTTTTTTCCTTATTTACATCCCCATAGATTATTTGAGTCAGATCACAAATACGCttttttc                           >  2:642061/1‑97 (MQ=255)
                                                                       aTTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCt                         >  1:1795596/1‑100 (MQ=255)
                                                                        ttCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTccc                          <  2:720136/99‑1 (MQ=255)
                                                                        ttCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTg                        >  1:2076534/1‑100 (MQ=255)
                                                                          cACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTgg                       <  2:601744/100‑1 (MQ=255)
                                                                          cACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTT‑TTTTTTTCCTTATTTAAATCGCCATAGATTATTTTCATCAGCTTACAAATACGCTTTTTccc                          >  2:43839/1‑96 (MQ=255)
                                                                              aaaaacggcgCATAAGCGCCG‑‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTaaaa                  >  1:2124452/6‑100 (MQ=255)
                                                                                gAGCGGCGCTTATGCGCCG‑‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTaaaaaa                >  2:1786755/1‑100 (MQ=255)
                                                                                aaacggcgCATAAGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTaaa                   >  1:2112403/4‑98 (MQ=255)
                                                                                 aGCGGCGCTTATGCGCCG‑TTTTTTTTT‑TTTTACACCTTATTTACATCCCCATAGATTATTTGCATCAACTCACAAATACGCTTTTTctttt                        >  2:1410711/1‑87 (MQ=255)
                                                                                 aGCGGCGCTTATGCGCCG‑‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAAt               <  1:927034/100‑1 (MQ=255)
                                                                                 aGCGGCGCTTATGCGCCG‑‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAAt               >  1:1700570/1‑100 (MQ=255)
                                                                                    ggcgCTTATGCGCCG‑‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGa             >  1:484474/1‑99 (MQ=255)
                                                                                    ggcgCGTATCGGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGGCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGa             <  1:941669/100‑1 (MQ=255)
                                                                                     gcgcTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGAt            <  1:906335/100‑1 (MQ=255)
                                                                                     gcgcTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGAt            <  1:1101969/100‑1 (MQ=255)
                                                                                       gcTTATGCGCCG‑‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATATATTATTTGCGTCAGCTCACTAATACTCTTTTTCACTAGTAAAAAATGAttt          <  2:1168583/99‑1 (MQ=255)
                                                                                        cTTATGCGCCG‑‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTcc        <  1:405187/100‑1 (MQ=255)
                                                                                         ttATGCGCCGTTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGAtt           >  1:858762/1‑98 (MQ=255)
                                                                                         ttATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTcc        <  2:1447993/100‑1 (MQ=255)
                                                                                         ttATGCGCCG‑‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCt       <  2:1134120/100‑1 (MQ=255)
                                                                                            tGCGCCG‑‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCg    <  2:2043518/100‑1 (MQ=255)
                                                                                             gcgcCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGc     >  2:1518399/1‑99 (MQ=255)
                                                                                             gcgcCG‑‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGt   <  2:343341/100‑1 (MQ=255)
                                                                                              cgcCG‑‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTg  <  2:344076/100‑1 (MQ=255)
                                                                                                    |                                                                                               
GGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTG  >  NC_000913/1369937‑1370131

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: ‑