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

GGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTG  >  NC_000913/1369937‑1370131
                                                                                                    |                                                                                               
gggATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGt                                                                                                 >  1:423285/1‑100 (MQ=255)
  gATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGttt                                                                                               >  1:739268/1‑100 (MQ=255)
    tataCCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGttttt                                                                                             >  2:501492/1‑100 (MQ=255)
      taCCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGttttttt                                                                                           <  1:1013788/100‑1 (MQ=255)
      taCCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGttttttt                                                                                           >  1:1002106/1‑100 (MQ=255)
      taCCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGttttttt                                                                                           >  1:1847984/1‑100 (MQ=255)
      taCCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGttttttt                                                                                           >  2:32702/1‑100 (MQ=255)
      taCCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGttttttt                                                                                           >  2:1631472/1‑100 (MQ=255)
        cccACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGttttttttc                                                                                         >  1:2062310/1‑99 (MQ=255)
         ccACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGttttttttct                                                                                        <  1:851138/100‑3 (MQ=255)
          cACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGttttttttctg                                                                                       >  2:2130728/1‑97 (MQ=255)
          cACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTc                                                                                        <  1:2171198/99‑1 (MQ=255)
           aCGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTg                                                                                      <  1:439244/100‑1 (MQ=255)
            cGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTc                                                                                    <  2:2188583/100‑1 (MQ=255)
             gTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTc                                                                                    <  2:423286/100‑1 (MQ=255)
             gTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTca                                                                                   >  1:403776/1‑100 (MQ=255)
             gTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTca                                                                                   >  1:2149643/1‑100 (MQ=255)
               tGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTcac                                                                                  <  2:932311/100‑1 (MQ=255)
               tGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTcac                                                                                  <  1:1006381/100‑1 (MQ=255)
                 agaACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCacac                                                                                <  2:558272/100‑1 (MQ=255)
                 agaACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACAcc                                                                               <  1:1095192/100‑1 (MQ=255)
                    aCGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTTTGTCTCTCCtt                                                                             >  2:969079/1‑100 (MQ=255)
                    aCGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCtt                                                                             >  1:720136/1‑100 (MQ=255)
                    aCGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTa                                                                            >  1:482876/1‑100 (MQ=255)
                       ccgatctccTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTAttt                                                                         <  1:2196643/93‑1 (MQ=255)
                        cGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTAttt                                                                         <  2:568710/100‑1 (MQ=255)
                            ggCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATc                                                                    >  2:1492036/1‑100 (MQ=255)
                             gCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATcc                                                                   >  1:1203890/1‑100 (MQ=255)
                              ccTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTAGATcgg                                                                  >  2:2196628/1‑98 (MQ=255)
                               cttaaaTACATTGCAATGCATAAGGTCAAAGTTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATccc                                                                  <  1:7856/97‑1 (MQ=255)
                                tGAAATACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCa                                                                <  2:2062296/100‑1 (MQ=255)
                                tGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTTCTGTCACACCTTATTTACATccc                                                                  <  1:1941051/100‑1 (MQ=255)
                                tGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCa                                                                >  1:1361824/1‑100 (MQ=255)
                                tGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCa                                                                >  1:1031527/1‑100 (MQ=255)
                                 gAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCa                                                                >  1:523432/1‑100 (MQ=255)
                                 gAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCAt                                                               >  2:1594206/1‑100 (MQ=255)
                                 gAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCAt                                                               <  2:403777/100‑1 (MQ=255)
                                   aaGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCACCAAAg                                                             >  1:809065/1‑100 (MQ=255)
                                   aaGACATCGCAATACCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAg                                                             <  1:590348/100‑1 (MQ=255)
                                       cATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACGATGAGCGGCGCTTATGCGCCGTTTTTTTTTTTTTTTATACCTTATCTACATCCCCATAGa                                                            >  2:1375140/1‑99 (MQ=255)
                                       cATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGatta                                                         >  1:1718919/1‑100 (MQ=255)
                                       cATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGatta                                                         >  1:343340/1‑100 (MQ=255)
                                       cATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGatta                                                         >  1:1161026/1‑100 (MQ=255)
                                        aTCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGatta                                                         >  1:302318/1‑100 (MQ=255)
                                         tCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGattatt                                                       >  1:608534/1‑100 (MQ=255)
                                          cGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTAttt                                                      >  1:1855460/1‑100 (MQ=255)
                                           gCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTAttt                                                      >  2:1975048/1‑100 (MQ=255)
                                           ccgatctcGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTg                                                     <  1:1114771/93‑1 (MQ=255)
                                            cAATGTCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGc                                                    <  1:654559/100‑1 (MQ=255)
                                             aaTGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCg                                                   <  2:1161022/100‑1 (MQ=255)
                                                 ccGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAg                                               >  2:1090663/1‑100 (MQ=255)
                                                 acggaGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCTATAGATTATTTgtgt                                                  >  2:1012057/5‑97 (MQ=255)
                                                  cGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAg                                               >  1:165202/1‑100 (MQ=255)
                                                  cGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTgagatcgg                                              >  2:1114768/1‑93 (MQ=255)
                                                    aaGTTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTc                                            <  1:1792349/100‑1 (MQ=255)
                                                     aGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTca                                           <  2:1031525/100‑1 (MQ=255)
                                                      ggTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTcac                                          >  2:651366/1‑100 (MQ=255)
                                                       gTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTcac                                          <  1:1331906/100‑1 (MQ=255)
                                                         cAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACaa                                        >  2:1963327/1‑100 (MQ=255)
                                                           aaGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTCCATCCCCATAGATTATTTGCGTCAGCTCACAAATa                                     >  2:23951/1‑100 (MQ=255)
                                                           aaGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATa                                     <  2:851389/100‑1 (MQ=255)
                                                            agtttAAAAGGATTCATCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGACACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATa                                     <  1:937592/97‑1 (MQ=255)
                                                            aGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGCTTATTTGCGTCCGCTCACAAATAc                                    >  2:957086/1‑100 (MQ=255)
                                                               ttAAAAGAATTGACCATGAGCGGCGCTGATGCGCCGTTTTTTTTCTGTCACTCCTTATTTACCTCCCCATAGATTATTTGCGTCAGCTCACAAATACGCt                                 <  1:1492046/100‑1 (MQ=255)
                                                               ttAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGc                                  <  2:224807/100‑1 (MQ=255)
                                                                tAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCtt                                <  1:1142016/100‑1 (MQ=255)
                                                                tAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCtt                                >  2:1859186/1‑100 (MQ=255)
                                                                  aaaGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCt                                 >  1:585041/1‑98 (MQ=255)
                                                                  aaaGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTTTTTTTTTTTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCtttt                              >  2:1053757/1‑100 (MQ=255)
                                                                      aaTTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTTTTTTTTCCTTATTTACATCCCCATAGATTATTTGAGTCAGATCACAAATACGCttttttc                           >  2:642061/1‑97 (MQ=255)
                                                                       aTTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCt                         >  1:1795596/1‑100 (MQ=255)
                                                                        ttCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTccc                          <  2:720136/99‑1 (MQ=255)
                                                                        ttCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTg                        >  1:2076534/1‑100 (MQ=255)
                                                                          cACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTgg                       <  2:601744/100‑1 (MQ=255)
                                                                          cACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTTTTTTTTCCTTATTTAAATCGCCATAGATTATTTTCATCAGCTTACAAATACGCTTTTTccc                          >  2:43839/1‑96 (MQ=255)
                                                                              aaaaacggcgCATAAGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTaaaa                  >  1:2124452/6‑100 (MQ=255)
                                                                                gAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTaaaaaa                >  2:1786755/1‑100 (MQ=255)
                                                                                aaacggcgCATAAGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTaaa                   >  1:2112403/4‑98 (MQ=255)
                                                                                 aGCGGCGCTTATGCGCCGTTTTTTTTTTTTTACACCTTATTTACATCCCCATAGATTATTTGCATCAACTCACAAATACGCTTTTTctttt                        >  2:1410711/1‑87 (MQ=255)
                                                                                 aGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAAt               <  1:927034/100‑1 (MQ=255)
                                                                                 aGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAAt               >  1:1700570/1‑100 (MQ=255)
                                                                                    ggcgCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGa             >  1:484474/1‑99 (MQ=255)
                                                                                    ggcgCGTATCGGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGGCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGa             <  1:941669/100‑1 (MQ=255)
                                                                                     gcgcTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGAt            <  1:906335/100‑1 (MQ=255)
                                                                                     gcgcTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGAt            <  1:1101969/100‑1 (MQ=255)
                                                                                       gcTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATATATTATTTGCGTCAGCTCACTAATACTCTTTTTCACTAGTAAAAAATGAttt          <  2:1168583/99‑1 (MQ=255)
                                                                                        cTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTcc        <  1:405187/100‑1 (MQ=255)
                                                                                         ttATGCGCCGTTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGAtt           >  1:858762/1‑98 (MQ=255)
                                                                                         ttATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTcc        <  2:1447993/100‑1 (MQ=255)
                                                                                         ttATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCt       <  2:1134120/100‑1 (MQ=255)
                                                                                            tGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCg    <  2:2043518/100‑1 (MQ=255)
                                                                                             gcgcCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGc     >  2:1518399/1‑99 (MQ=255)
                                                                                             gcgcCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGt   <  2:343341/100‑1 (MQ=255)
                                                                                              cgcCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTg  <  2:344076/100‑1 (MQ=255)
                                                                                                    |                                                                                               
GGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTG  >  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: