Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,531,596 (G)6→5 coding (134/762 nt) cysZ → sulfate transporter, sulfite inhibited

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,531,5910G.93.8% 276.8 / 12.7 65coding (129/762 nt)cysZsulfate transporter, sulfite inhibited
Reads supporting (aligned to +/- strand):  ref base G (2/2);  new base . (23/38);  total (25/40)
Fisher's exact test for biased strand distribution p-value = 6.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

CAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGTTATCT  >  NC_000913/2531492‑2531686
                                                                                                   |                                                                                               
cAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATg                                                                                                 <  2:1683694/100‑1 (MQ=255)
cAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATg                                                                                                 <  1:1904227/100‑1 (MQ=255)
 aGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATgg                                                                                                >  1:60768/1‑100 (MQ=255)
  gCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCCATATTTTGTTGATggg                                                                                               <  1:2085709/100‑1 (MQ=255)
   cGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATgggg                                                                                              >  2:1321173/1‑100 (MQ=255)
     gTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATgggggc                                                                                            <  2:2058217/100‑2 (MQ=255)
            cTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCt                                                                                    >  1:1117713/1‑100 (MQ=255)
            cTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCt                                                                                    >  1:654126/1‑100 (MQ=255)
            cTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCt                                                                                    >  1:1540268/1‑100 (MQ=255)
              aTTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCtgg                                                                                  <  2:674543/100‑1 (MQ=255)
              aTTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCtgg                                                                                  >  2:1025207/1‑100 (MQ=255)
               ttttGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCtggt                                                                                 >  1:357975/1‑100 (MQ=255)
                tttGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCtggtg                                                                                <  1:761750/100‑1 (MQ=255)
                tttGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCtggtg                                                                                <  2:174537/100‑1 (MQ=255)
                tttGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCtggtg                                                                                <  2:572926/100‑1 (MQ=255)
                 ttGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCtggtg                                                                                <  2:412105/100‑1 (MQ=255)
                    cgcAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGctc                                                                            >  1:73642/1‑100 (MQ=255)
                     gcAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGctct                                                                           >  1:1844064/1‑100 (MQ=255)
                      cAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCtt                                                                          <  2:1660138/100‑1 (MQ=255)
                      cAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCtt                                                                          <  2:1144176/100‑1 (MQ=255)
                        aGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTa                                                                        <  1:281917/100‑1 (MQ=255)
                         ggCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTac                                                                       <  1:418943/100‑1 (MQ=255)
                           cTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTacac                                                                     <  1:76864/100‑1 (MQ=255)
                           cTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTacac                                                                     <  1:1943739/100‑1 (MQ=255)
                            tGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTTGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTAcaca                                                                    <  1:1239230/100‑1 (MQ=255)
                            tGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTacac                                                                     <  1:1172685/100‑1 (MQ=255)
                            tGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTAcaca                                                                    <  2:1353290/100‑1 (MQ=255)
                             ggAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTAcaca                                                                    <  2:949955/99‑1 (MQ=255)
                             ggAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAg                                                                   >  2:796818/1‑100 (MQ=255)
                                aGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTc                                                                <  2:1117363/100‑1 (MQ=255)
                                 gCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCg                                                               >  1:1026870/1‑100 (MQ=255)
                                  cTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGa                                                              >  2:798511/1‑100 (MQ=255)
                                   tCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGAt                                                             >  1:585125/1‑100 (MQ=255)
                                    cGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATg                                                            <  2:156745/100‑1 (MQ=255)
                                    cGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATg                                                            <  2:156850/100‑1 (MQ=255)
                                      tctcGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTc                                                          >  1:1847315/1‑100 (MQ=255)
                                        tcGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTg                                                        <  2:107857/100‑1 (MQ=255)
                                        tcGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTg                                                        <  1:1729362/100‑1 (MQ=255)
                                        tcGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTg                                                        <  2:357976/100‑1 (MQ=255)
                                         cGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTg                                                        >  1:1438451/1‑100 (MQ=255)
                                          gCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCAGTCTGGTGGCTCTTTACACAGCTCGATGTCTGGa                                                      >  2:1928619/1‑100 (MQ=255)
                                           cAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGAt                                                     <  2:804880/100‑1 (MQ=255)
                                           cAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGAt                                                     <  1:1825500/100‑1 (MQ=255)
                                              ccTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATccc                                                  <  1:1505371/100‑1 (MQ=255)
                                              ccTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATccc                                                  <  1:514965/100‑1 (MQ=255)
                                               cTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCg                                                 >  1:1112647/1‑100 (MQ=255)
                                               cTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCg                                                 >  2:1607793/1‑100 (MQ=255)
                                                      tCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCa                                          <  2:1250326/100‑1 (MQ=255)
                                                        ggCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACTCAGCTCGATGTCTGGATCCCGACTCTCATg                                        <  2:60768/100‑1 (MQ=255)
                                                        ggCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATg                                        <  2:1695028/100‑1 (MQ=255)
                                                        ggCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATg                                        >  1:554907/1‑100 (MQ=255)
                                                            tttCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGtt                                    >  2:1457652/1‑100 (MQ=255)
                                                               cGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACg                                 >  2:287068/1‑100 (MQ=255)
                                                                attATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGt                                <  2:1529214/99‑1 (MQ=255)
                                                                    ttttACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCg                            >  1:97928/1‑100 (MQ=255)
                                                                    ttttACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCg                            <  2:588928/100‑1 (MQ=255)
                                                                       tACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGAc                         >  1:287097/1‑100 (MQ=255)
                                                                           gctgctGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGc                     <  2:899341/100‑1 (MQ=255)
                                                                            ctgctgGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCt                    >  2:1470581/1‑100 (MQ=255)
                                                                              gctgGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTAc                  <  1:288839/100‑1 (MQ=255)
                                                                              gctgGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTAc                  <  2:903519/100‑1 (MQ=255)
                                                                                    cAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGc            <  2:939292/100‑1 (MQ=255)
                                                                                     aaTATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCt           <  2:1283253/100‑1 (MQ=255)
                                                                                     aaTATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCt           <  2:301828/100‑1 (MQ=255)
                                                                                      atatTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTg          >  1:843963/1‑100 (MQ=255)
                                                                                       tatTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGa         <  1:1400300/100‑1 (MQ=255)
                                                                                         ttttGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGt       <  1:992228/100‑1 (MQ=255)
                                                                                           ttgttgAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGTTa     >  2:205310/1‑100 (MQ=255)
                                                                                             gttgAT‑GGGGGCGCGTTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGTTATc   <  1:502921/100‑1 (MQ=255)
                                                                                             gttgAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGTTATc   <  2:871140/100‑1 (MQ=255)
                                                                                               tgATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCGTGAGTTACGTTCCGGACTGGCTACAATGGCTGAGTTATCt  >  2:723286/1‑100 (MQ=255)
                                                                                                   |                                                                                               
CAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGTTATCT  >  NC_000913/2531492‑2531686

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