Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,028,159 T→A 16.9% intergenic (‑142/+29) yedQ ← / ← yodC putative membrane‑anchored diguanylate cyclase/uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,028,1590TA16.9% 194.4 / 29.5 77intergenic (‑142/+29)yedQ/yodCputative membrane‑anchored diguanylate cyclase/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (27/37);  new base A (4/9);  total (31/46)
Fisher's exact test for biased strand distribution p-value = 5.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.52e-01

CGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGA  >  NC_000913/2028064‑2028258
                                                                                               |                                                                                                   
cGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATc                                                                                                 >  1:2710130/1‑100 (MQ=255)
 gATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATcc                                                                                                <  2:1871456/100‑1 (MQ=255)
   tAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCgg                                                                                              <  1:2288658/100‑1 (MQ=255)
     gCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTA‑‑ACATCAGGCCGGATCACATCCGGCCt                                                                                           <  1:1314623/99‑1 (MQ=37)
      cGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCt                                                                                           <  2:2092825/100‑1 (MQ=255)
      cGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCt                                                                                           <  2:1465445/100‑1 (MQ=255)
      cGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATCACATCCGGCCt                                                                                           <  1:2705276/100‑1 (MQ=255)
       gTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTg                                                                                          <  2:27236/100‑1 (MQ=255)
       gTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTg                                                                                          <  1:1606336/100‑1 (MQ=255)
       gTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATCACATCCGGccn                                                                                           <  1:2581947/99‑2 (MQ=255)
       gTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATCACATCCGGCCTg                                                                                          <  2:2610131/100‑1 (MQ=255)
        tACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACATCAGGCCGGATCACATCCGGCCTgg                                                                                         <  1:2674046/100‑2 (MQ=38)
        tACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGa                                                                                         <  2:2179293/100‑1 (MQ=255)
        tACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATCACATCCGGCCTgg                                                                                         <  1:2662910/100‑2 (MQ=255)
        tACCGCAAAACAGCGATAATGTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATCACATCCGGCCTgg                                                                                         <  1:1311595/100‑2 (MQ=255)
         aCCGCTAACCCGCGATAACCTTATGATGGCATCACAATAAGACGCGACTATAACTCATTGCCCGCATTAGCACACCAGGCCCGATCACATCCGGCCtag                                                                                         <  1:1738703/99‑3 (MQ=255)
          ccGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATg                                                                                       <  1:2262302/100‑1 (MQ=255)
          ccGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATg                                                                                       >  2:1208200/1‑100 (MQ=255)
            gCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGtt                                                                                     >  2:2852745/1‑100 (MQ=255)
              aaaaCAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTaa                                                                                   <  2:566588/100‑1 (MQ=255)
               aaaCAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAAt                                                                                  <  2:1581921/100‑1 (MQ=255)
                aaCAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATg                                                                                 >  2:8027/1‑100 (MQ=255)
                aaCAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATg                                                                                 <  2:2199636/100‑1 (MQ=255)
                  cAGCGATAATTTTATAATGGTGGCATGATGAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATCACATCCGGCCTGgtgtgctaatg                                                                               <  2:740457/100‑8 (MQ=25)
                   aGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACAACAGGCCGGATGTGATCCGGCCTGATGTTAATGATa                                                                              <  1:2701554/100‑1 (MQ=255)
                    gCGAAAATTTTATAATGGTAGCATGAGAAGACGCAGCTATCGTTCATTGGAGGCAGTAGCAGACCAGGCGGGGTGTGATCCGGCCTGATGTTAATGATa                                                                              >  1:2431096/1‑99 (MQ=255)
                          aTTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTc                                                                       >  1:87738/1‑100 (MQ=255)
                            tttATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTcac                                                                     <  2:147059/100‑1 (MQ=255)
                               ataATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACAct                                                                  >  1:2048960/1‑100 (MQ=255)
                                taATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACActt                                                                 >  1:508757/1‑100 (MQ=255)
                                  aTGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCCGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACActtct                                                               <  1:1266032/100‑1 (MQ=255)
                                    ggtggCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACttcttc                                                             <  2:524279/100‑1 (MQ=255)
                                    ggtggCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACttcttc                                                             >  2:2440396/1‑100 (MQ=255)
                                     gtggCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTcc                                                            >  1:1536432/1‑100 (MQ=255)
                                     gtggCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTcc                                                            <  2:613592/100‑1 (MQ=255)
                                      tggCTTTATAAGTCGCGACTTTCATTCATTCCAGACTTTATCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATCATACCCTTTCACACttcttc                                                             <  1:474728/98‑1 (MQ=255)
                                       ggCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGc                                                          <  1:2009363/100‑1 (MQ=255)
                                       ggCAGGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGc                                                          <  1:2544962/100‑1 (MQ=255)
                                             aTAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACg                                                    <  1:1926278/100‑1 (MQ=255)
                                             aTAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACg                                                    <  2:155846/100‑1 (MQ=255)
                                                aGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGAct                                                 <  2:1317523/100‑1 (MQ=255)
                                                 gACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGActc                                                <  1:677062/100‑1 (MQ=255)
                                                  aCGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGActcc                                               <  2:1438251/100‑1 (MQ=255)
                                                           tcattcatTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCg                                      <  2:852702/100‑1 (MQ=255)
                                                           tcattcatTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCACCCCCg                                      <  1:286870/100‑1 (MQ=255)
                                                            cattcattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCgg                                     >  2:1603082/1‑100 (MQ=255)
                                                             attcattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGa                                    >  1:2485233/1‑100 (MQ=255)
                                                              ttcattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGaa                                   <  1:448949/100‑1 (MQ=255)
                                                              ttcattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGaa                                   <  1:1706686/100‑1 (MQ=255)
                                                               tcattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAAc                                  <  2:1997112/100‑1 (MQ=255)
                                                                cattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACa                                 >  1:2844083/1‑100 (MQ=255)
                                                                 nttGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCACAGAACGACTCCCCTCCCCCGGAACaa                                >  2:1278265/2‑100 (MQ=255)
                                                                  ttGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACACTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAg                               >  2:971934/1‑100 (MQ=255)
                                                                    gCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCt                             >  2:2318489/1‑100 (MQ=255)
                                                                      tggCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCg                           <  1:200200/99‑1 (MQ=255)
                                                                       ggCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCgg                          >  1:153037/1‑100 (MQ=255)
                                                                       ggCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCgg                          <  1:1603076/100‑1 (MQ=255)
                                                                       ggCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCga                          >  2:200201/1‑99 (MQ=255)
                                                                         cATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGtt                        >  1:2408077/1‑100 (MQ=255)
                                                                          aTTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGttt                       <  1:916119/100‑1 (MQ=255)
                                                                           ttAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCAGGAACAAGCTCGGTTTc                      >  1:2173563/1‑100 (MQ=255)
                                                                            tAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCGCGGACCAAGCTCGGTTTCa                     >  1:313805/1‑100 (MQ=255)
                                                                            attaacatcAGGCCGGATCACATCCGGCCTGGTGTTAATGATACCCTTTCACCCTTCTTCCGCAGAACGACTCCAATCCCCCGGGACAAGCTCGGTTTCa                     >  1:1395607/9‑100 (MQ=21)
                                                                            attaacatcAGGCCGGATCACATCCGGCCTGGTGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCa                     >  2:2236457/9‑100 (MQ=25)
                                                                             ttaacatcAGGCCGGATCACATCCGGCCTGGTGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCTCCTCCCCCGGAACAAGCTCGGTTTCAt                    >  1:2701088/8‑100 (MQ=25)
                                                                              gCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATg                   <  2:2277311/100‑1 (MQ=255)
                                                                              gCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATg                   >  1:2279952/1‑100 (MQ=255)
                                                                              gCACACCAGGCCGGATGTGATCCGGCCTGATGTTAAGGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATg                   <  2:2048971/100‑1 (MQ=255)
                                                                                acacCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATg                   >  1:2124009/1‑98 (MQ=255)
                                                                                acacCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGaa                 >  2:2475565/1‑100 (MQ=255)
                                                                                 catcAGGCCGGATCACATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGaaa                >  2:23776/4‑100 (MQ=38)
                                                                                   ccAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAg              <  2:1215490/100‑1 (MQ=255)
                                                                                     aGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCt            >  2:1023691/1‑100 (MQ=255)
                                                                                     aGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCt            >  2:1981149/1‑100 (MQ=255)
                                                                                      ggCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCtt           >  2:219023/1‑100 (MQ=255)
                                                                                       gCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTc          <  2:2385824/100‑1 (MQ=255)
                                                                                            aTGTGATCCGGCCTGATGTTAACGATACCCCTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCTGAACAAGCTCGGTTTCATGAAAAGCTTc          <  1:2719034/95‑1 (MQ=255)
                                                                                             tgtgATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGCACGACTCCCCTCCCCCGGAACATGCTCGGTTTCATGACCAGCTTCCCGCtt    >  1:2098408/1‑100 (MQ=255)
                                                                                             tgtgATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCtt    >  2:65498/1‑100 (MQ=255)
                                                                                              gtgATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTg   >  2:2437840/1‑100 (MQ=255)
                                                                                              gtgATCCGGCCTGATATGAATGATACCCTGTCACACTTCGTCCGCAGAACGGCTCCCCTCCCCCGGAACAGGCTCTGTTGCGTGAAAAGCTTCCCGCt     <  1:1304790/98‑1 (MQ=255)
                                                                                               tgATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGa  >  2:2231458/1‑100 (MQ=255)
                                                                                               |                                                                                                   
CGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGA  >  NC_000913/2028064‑2028258

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: