Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,370,044 (T)9→8 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.89.7% 202.7 / 22.7 59intergenic (+33/‑180)pspE/ycjMthiosulfate:cyanide sulfurtransferase (rhodanese)/alpha amylase catalytic domain family protein
Reads supporting (aligned to +/- strand):  ref base T (4/2);  new base . (28/24);  total (33/26)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTT‑TTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTG  >  NC_000913/1369937‑1370131
                                                                                                   |                                                                                                
gggATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGt                                                                                                  >  2:1303238/1‑100 (MQ=255)
 ggATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAAAATTCACCATGAGCGGCGCTTATGCGCCGtt                                                                                                 >  2:515838/1‑100 (MQ=255)
        cccACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGttt‑tttttc                                                                                         <  2:1326777/100‑2 (MQ=255)
         ccACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGttt‑tttttct                                                                                        <  1:920444/100‑3 (MQ=255)
           aCGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGttt‑tttttctgt                                                                                      <  2:940913/100‑5 (MQ=255)
            cGTTGAGCACTCCGGTGGCATGAAAGACATCGCAATGCCGAAAGTCAAAAGTTAAAAGAATGCACCATGAGCGGCGCTTATGTGCCG‑TT‑TTTTTTCTGTc                                                                                    <  1:444155/100‑1 (MQ=255)
            cGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGttt‑tttttct                                                                                        >  2:66798/1‑95 (MQ=255)
            cGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TT‑TTTTTTCTGTc                                                                                    >  2:915866/1‑100 (MQ=255)
            cGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TT‑TTTTTTCTGTc                                                                                    <  2:480004/100‑1 (MQ=255)
                gagaACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATTCGTCG‑TT‑TTTTTTCTGTCacac                                                                                <  1:393343/100‑1 (MQ=255)
                 agaACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TT‑TTTTTTCTGTCACAcc                                                                               <  2:279595/100‑1 (MQ=255)
                 agaACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TT‑TTTTTTCTGTCACAcc                                                                               >  2:1352790/1‑100 (MQ=255)
                     cGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TT‑TTTTTTCTGTCACACCTTAt                                                                           <  1:1303228/100‑1 (MQ=255)
                     atctggTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTT‑TTTTTTCTGTCACACCTTa                                                                            <  1:295777/96‑1 (MQ=255)
                      gCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TT‑TTTTTTCTGTCACACCTTAtt                                                                          <  1:128268/100‑1 (MQ=255)
                        cGGTGGCCCGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATTAGCGGCGCTTATGCGCCG‑TT‑TTTTTTCTGTCACACCTTAtttn                                                                        <  1:132951/100‑2 (MQ=255)
                         gtttgCCTTAAAGATATTGGAATGCCGAAGGTCAATGGTTATAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TT‑TTTTTTCTGTCACACCTTAtt                                                                          <  2:269326/93‑1 (MQ=255)
                         ggTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTT‑TTTTTTCTGTCACACCTTaaga                                                                         >  2:295779/1‑96 (MQ=255)
                          gTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TT‑TTTTTTCTGTCACACCTTATTTACa                                                                      <  2:924608/100‑1 (MQ=255)
                             gCCTGCAAGACATAGCCATGCCGAAGGTCAAAGGTTCAAAGAATTCAACATGCGCGGCGCTTATGCGCCGTTT‑TTTTTTCTGTTAAAAATTAtttttatc                                                                    >  2:360109/1‑95 (MQ=255)
                              ccTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TT‑TTTTTTCTGTCACACCTTATTTACATcc                                                                   >  1:240793/1‑99 (MQ=255)
                               cTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TT‑TTTTTTCTGTCACACCTTATTTACATcccc                                                                 >  1:1614513/1‑100 (MQ=255)
                                      aCATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TT‑TTTTTTCTGTCACACCTTATTTACATCCCCATAGatt                                                          <  1:770482/100‑1 (MQ=255)
                                        aTCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TT‑TTTTTTCTGTCACACCTTATTTACATCCCCATAGcttat                                                        >  1:702795/1‑100 (MQ=255)
                                        aTCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TT‑TTTTTTCTGTCACACCTTATTTACATCCCCATAGattat                                                        >  2:53251/1‑100 (MQ=255)
                                         tCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TT‑TTTTTTCTGTCACACCTTATTTACATCCCCATAGattatt                                                       >  1:1793519/1‑100 (MQ=255)
                                            cAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TT‑TTTTTTCTGTCACACCTTATTTACATCCCCATAGattatt                                                       >  2:440991/1‑97 (MQ=255)
                                               tGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TT‑TTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTc                                                 >  2:1462245/1‑100 (MQ=255)
                                               tGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TT‑TTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTc                                                 <  1:255307/100‑1 (MQ=255)
                                               tGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TT‑TTTTTTCTGTCAAACCTTATTTACATCCCCATAGATTATTTGCGTc                                                 >  1:63286/1‑100 (MQ=255)
                                                gCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TT‑TTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCa                                                >  2:167495/1‑100 (MQ=255)
                                                 ccGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TT‑TTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAg                                               >  2:929539/1‑100 (MQ=255)
                                                  cGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TT‑TTTTTTCTGTCACACCTGATTTACATGCCCATAGATTATTTGCGTCAGc                                              >  1:1223029/1‑100 (MQ=255)
                                                   gAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TT‑TTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCt                                             <  2:813477/100‑1 (MQ=255)
                                                     aGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TT‑TTTTTTCTGTCGCACCTTATTTACATCCCCATAGATTATTTGCGTCAGc                                              >  1:1340316/1‑97 (MQ=255)
                                                      ggTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TT‑TTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTcac                                          >  1:457679/1‑100 (MQ=255)
                                                      ggTCAAAGGTTAAAAGAATTCACCATGAGCGGCGATTATGCGCCG‑TT‑TTTTTTCTGTCACACCTTATTTCCCTCCCCATAGATAATTTGCGTCAGCTcac                                          >  1:1710971/1‑100 (MQ=255)
                                                       gTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TT‑TTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTcaca                                         <  2:80641/100‑1 (MQ=255)
                                                       gTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑AT‑TTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTcaca                                         >  1:634890/1‑100 (MQ=255)
                                                         cAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TT‑TTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACaaa                                       >  2:609484/1‑100 (MQ=255)
                                                         cAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TT‑TTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACaaa                                       >  2:1455910/1‑100 (MQ=255)
                                                           aaGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TT‑TTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATa                                     <  2:791886/100‑1 (MQ=255)
                                                             ggTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TT‑TTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACg                                   >  2:762750/1‑100 (MQ=255)
                                                               ttAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TT‑TTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCt                                 <  1:304476/100‑1 (MQ=255)
                                                                 aaaaGAATTCACCATGAGCGGCGCTTATGCGCCG‑TT‑TTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCttt                               <  1:736288/100‑1 (MQ=255)
                                                                  aaaGAATTCACCATGAGCGGCGCTTATGCGCCG‑TT‑TTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCtttt                              >  2:184148/1‑100 (MQ=255)
                                                                  aaaGAATTCACCATGAGCGGCGCTTATGCGCCG‑TT‑TTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCtttt                              >  2:967602/1‑100 (MQ=255)
                                                                   aaGAATTCACCATGAGCGGCGCTTATGCGCCG‑TT‑TTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCttttt                             <  2:11204/100‑1 (MQ=255)
                                                                    aGAATTCACCATGAGCGGCGCTTATGCGCCG‑TT‑TTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTc                            <  2:1793532/100‑1 (MQ=255)
                                                                         agaaaaaaaacggcgCATAAGCGCCG‑TT‑TTTTTTCTGTCACACCTTATTTACATCCCCATAGATCATTTGCGTCAGCTCACAAATACGCTTTTTCCCTgg                       >  1:556410/11‑100 (MQ=255)
                                                                           aCCATGAGCGGCGCTTATGCGCCG‑TT‑TTTTTTCTGTCACTCCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTa                     >  1:927476/1‑100 (MQ=255)
                                                                            ccATGAGCGGCGCTTATGCGCCG‑TT‑TTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTaa                    >  1:752572/1‑100 (MQ=255)
                                                                            ccATGAGCGGCGCTTATGCGCCG‑TT‑TTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTaa                    >  2:1506368/1‑100 (MQ=255)
                                                                             cATGAGCGGCGCTTATGCGCCGTTT‑TTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTaa                    >  2:1257516/1‑100 (MQ=255)
                                                                                 aacggcgCATAAGCGCCGCTT‑TTTTTTTTTTTTTTCCTTATTTACATCCCCATATATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTTGTaaaaa                 >  1:1775962/3‑99 (MQ=255)
                                                                                 aGCGGCGCTTATGCGCCG‑TT‑TTTCTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTaaaaaa                <  1:1064925/99‑1 (MQ=255)
                                                                                 aGCGGCGCATAAGCGCCG‑TT‑TTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTTCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAAt               >  2:1685545/1‑100 (MQ=255)
                                                                                  gcggcgCTTATGCGCCG‑TT‑TTTTTTCTGTTACACCTTATTTTCATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAAt               <  2:136957/99‑1 (MQ=255)
                                                                                   cggcgCTTATGCGCCGTTTACTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAAt               <  2:800630/100‑1 (MQ=255)
                                                                                   cggcgCATAAGCGCCGTTT‑TTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAAt               >  2:1638037/1‑99 (MQ=255)
                                                                                   ccgcgcTTATGCGCCG‑TT‑TTTTTTCTTTCACGACTTATTTACATCCCCATCGATTATTTGCGTCAGCTCACATATGCGCTTTTTCCCTGGTAAAAAAtt              <  1:677891/97‑2 (MQ=255)
                                                                                    ggcgCTTATGCGCCG‑TT‑TTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGAt            <  1:1363298/100‑1 (MQ=255)
                                                                                     gtgcTTATGCGCCG‑TT‑TTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGAt            <  1:1353654/97‑1 (MQ=255)
                                                                                       gcTTATGCGCCG‑TT‑TTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGAttt          >  2:1261518/1‑99 (MQ=255)
                                                                                            atgcgctttt‑ttttttCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGc     <  2:1804353/93‑1 (MQ=255)
                                                                                              cgcCG‑TT‑TTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTg  <  2:36630/100‑1 (MQ=255)
                                                                                                   |                                                                                                
GGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTT‑TTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTG  >  NC_000913/1369937‑1370131

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑